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An  Analysis  of  the  Data  Collection  Modes  of  a 
Digital  Weather  Radar  System  with  Respect  to 
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B.S.,  University  of  Utah 

Chairman  of  Advisory  Committee:  Dr.  George  Huebner 


Dual  wavelength  digital  weather  radar  data  were  collected  at  zero 
deg  elevation  angle  in  all  four  modes  of  collection  possible  with  the 
Texas  A&M  University  weather  radar  system  (1  and  2 deg  and  1 and  2 km 
modes) . These  data  were  processed  and  displayed  as  contoured  maps  of 
reflectivity.  These  maps  were  then  analyzed,  both  qualitatively  and 
quantitatively,  in  an  effort  to  determine' an  optimum  mode  of  collection 
with  respect  to  retention  of  meteorological  features  and  data  processing 
resources  required  to  produce  the  maps. 

The  results  indicated  that  significant  differences  exist  between 
the  reflectivity  displays  obtained  from  the  different  collection  modes. 
Results  from  two  of  the  modes  are  good,  while  the  other  two  modes  suf- 
fer from  enough  loss  of  detail  and  accuracy  to  make  their  usefulness 
questionable . 

Additionally,  the  Texas  A&M  University  digital  weather  radar  system 
and  the  two  computer  programs  used  to  process  the  data  were  explained. 
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CHAPTER  I 


INTRODUCTION 


Within  only  a short  time  of  the  initial  development  of  radar,  it 
was  recognized  that  the  device  could  locate  and  track  areas  of  severe 
weather  and  precipitation.  Since  that  time,  radar  has  come  to  be 
considered  one  of  the  most  important  tools  available  to  the  meteoro- 
logist for  the  detection  of  severe  local  weather  phenomena,  such  as 
tornadoes,  hailstorms,  and  severe  thunderstorms.  Until  recently, 
however,  the  usefulness  of  radar  data  was  limited  by  the  qualitative, 
rather  than  quantitative,  nature  of  the  interpretation  to  which  it 
lent  itself.  This  interpretation  was  subject  to  many  variables, 
especially  the  experience  level  of  the  observer  and  his  ability  to 
recognize  correctly  the  type  and  severity  of  a storm  by  a series  of 
variable  radar  signatures.  The  development  of  digital  radar  systems 
represents  a significant  and  important  technological  advance.  When 
this  advance  is  coupled  with  an  efficient  and  responsive  means  of  data 
reduction  and  display,  it  will  afford  the  meteorologist  the  opportunity 
quantitatively  to  observe  and  analyze  severe  weather  phenomena  in  near 
real  time. 


The  Need  for  This  Investigation 


A great  amount  of  research  has  been  done  recently  using  digitized 
radar  data,  and  the  National  Weather  Service  (NWS)  has  a plan  to  in- 
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stall  digitizers  on  many  of  their  weather  radar  systems  in  the  United 
States  to  provide  input  to  the  new  Automation  of  Field  Operations  and 
Services  (AFOS)  system.  However,  the  operational  implementation  of 
digital  weather  radar  systems  still  faces  some  problems.  One  of  these 
is  the  problem  of  how  to  assimilate,  reduce,  and  display  in  a usable 
form,  the  massive  amount  of  data  generated  by  a digital  radar  system, 
and  be  able  to  utilize  the  data  display  in  near  real  time.  Until  now, 
one  of  the  most  widely  accepted  real  time  uses  of  these  data  has  been 
in  the  computation  of  vertical  integration  of  liquid  water  content 
(VIL) , a technique  developed  by  Greets  (1971) . While  it  is  agreed  by 
Vogel  (1973),  Can ip e (1973),  and  Pittman  (1976)  that  VIL  may  indeed  be 
an  indicator  of  the  potential  for  development  of  severe  storms,  it  is 
also  agreed  that  the  technique  does  not  adequately  represent  the  im- 
portant three-dimensional  aspects  of  a severe  storm.  Both  Canipe's 
partial  vertical  integration  of  liquid  water  (PVIL)  and  Pittman's 
partial  vertically-summed  reflectivity  (Z)  maps  (PVSZ)  better  accomp- 
lish this  task.  However,  both  of  these  techniques  are  more  difficult 
to  accomplish  in  real  time.  Sieland  (1977)  developed  a computer  pro- 
gram that  makes  it  possible  to  process  digital  radar  data  and  produce 
PVIL  and  PVSZ  displays  in  near  real  time.  The  operational  usefulness 
of  this  program  could  possibly  be  greatly  enhanced  by  several  additional 
refinements.  One  such  refinement  is  the  determination  of  an  optimum 


mode  of  data  collection 


Objectives  of  This  Investigation 


The  objectives  of  this  investigation  were  two-fold.  The  primary 
objective  was  to  examine  and  compare  dual  wavelength  digital  radar  data 
collected  near-simultaneously  in  the  four  modes  of  collection  possible 
by  the  Texas  A&M  University  weather  radar  system,  for  the  purpose  of 
determining  an  optimum  mode  of  data  collection  for  operational  use  of 
the  radar  system.  This  determination  was  based  on  such  considerations 
as  the  range  of  the  storm  from  the  radar,  wavelength  of  the  radar,  and 
data  processing  requirements,  both  hardware  and  software.  These  con- 
siderations have  been  weighed  against  those  of  resolution  of  data,  re- 
tention of  detail,  and  smoothing  of  fine  scale  meteorological  features. 
Digital  data  collected  from  the  3 and  10  cm  radars  of  the  Department 
of  Meteorology  at  TAMU  have  been  processed  by  Sieland 1 s (1977)  computer 
program,  modified  to  produce  only  the  zero-tilt  reflectivity  maps. 

These  displays  have  provided  the  basis  for  the  comparisons  and  evalua- 
tions of  the  different  data  collection  modes. 

A secondary  objective  of  this  investigation  was  to  provide  a 
detailed  description  of  the  composition  of  the  TAMU  digital  weather 
radar  system  itself. 


Present  Status  of  the  Question 

The  importance  of  weather  radar  as  a means  of  detecting,  observing 
and  forecasting  severe  weather,  and  in  providing  timely  warning  of  the 
approach  of  that  severe  weather  to  concerned  parties,  is  historically 
well  established.  Until  fairly  recently,  this  use  of  the  radar  was  re- 
stricted to  radar  scope  pattern-recognition  techniques,  many  of  which 
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were  summarized  by  Whiton  (1971).  These  patterns  included:  the  hook- 
shaped echo,  which  is  normally  shaped  like  the  figure  "6";  the  V-shaped 
notch;  pendants  or  protrusions,  normally  from  the  right  rear  portion  of 
a storm;  and  fingers,  or  scalloped  edges  of  the  echoes,  protruding  from 
the  rear  portion  of  the  storm.  Later  investigations  focused  attention 
on  the  vertical,  as  well  as  the  horizontal,  structure  of  severe  storms. 
Such  features  as  the  height  of  the  echo  tops  in  relation  to  the  tropo- 
pause  (Pautz , 1963) , echo-free  vaults  or  the  weak  echo  region  (WER) 
(Bigler,  1955),  and  high  reflectivity  spikes  protruding  from  the  tops 
of  severe  storms  (Yates,  1963).  All  of  these  techniques,  however,  rely 
entirely  upon  the  ability  of  the  radar  operator  to  observe  adequately 
these  features,  and  this  ability  varies  greatly  with  experience  and  ex- 
posure to  severe  storms.  Additionally,  the  value  of  qualitative  severe 
weather  indicators  depends  on  the  operator's  ability  to  assimilate 
mentally  information  pertaining  to  both  the  vertical  and  horizontal 
structure  and  extent  of  the  storm,  an  exceedingly  difficult  task.  The 
recent  introduction  of  electronic  digitizing  systems  for  weather  radar 
has  eased,  somewhat,  the  burden  on  radar  operators,  while  it  has 
created  a requirement  for  a method  of  reducing  and  analyzing  the  data 
in  near  real  time. 

Digital  Radar  Data 

Digitizing  systems  quantize  the  reflected  power  returned  to  the 
radar  receiver  over  a prescribed  number  of  pulses,  and  store  the  data 
in  a spherical  coordinate  system  of  range  (r) , azimuth  (a) , and 
antenna  elevation  angle  (e) . A typical  data  set  for  each  scan  sequence 
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of  a particular  storm  system  may  consist  of  as  many  as  100  radials 
(each  with  200  range  gates)  per  elevation  angle,  and  up  to  10  eleva- 
tion angles.  Since  it  is  necessary  to  sample  a storm  continuously  (one 
scan  every  5 to  10  minutes),  the  amount  of  data  collected  quickly  be- 
comes unmanageable;  consequently,  reduction  and  display  of  the  data 
requires  the  use  of  computer  methods. 

One  of  the  earliest  efforts  to  display  digital  data  was  developed 
by  Marshall  (1957),  who  processed  manually  digitized  photographic  data 
and  displayed  it  in  a constant  altitude  plan-position-indicator  (CAPPI) 
format.  Greene  (1971)  developed  a computer  technique  to  produce  con- 
stant altitude  reflectivity  (Z)  maps  (CAZM) . With  this  technique, 
digital  data  were  linearly  interpolated  from  a spherical  coordinate 
system  (r,a, e)  to  a cylindrical  coordinate  system  (s,a,h),  where  s is 
the  distance  along  a flat  earth  surface  and  h is  the  height  above 
ground.  Then  the  data  were  interpolated  in  the  vertical  plane  of  a, 
and  transformed  to  a rectangular  coordinate  system  (x,y,h)  through  a 
quadratic  interpolation  scheme.  The  result  was  a plan  view  of 
interpolated  radar  reflectivity  factors  in  constant  height  planes  with 
the  data  points  aligned  vertically  above  the  range  gates  in  the  zero- 
tilt,  or  0-deg  elevation,  level.  These  CAZMs  were  constructed  in 
intervals  of  5 kft  from  0-deg  tilt  up  to  50  kft  and  proved  to  be  a 
valuable  tool  for  analyzing  severe  storm  data.  Greene  also  originated 
the  concept  of  integrating  the  liquid  water  content  of  a storm  through 
its  depth,  and  labeled  the  technique  "vertical  integration  of  liquid 
water  content"  (VIL) . This  technique  was  an  attempt  to  compress  the 
three  dimensional  aspects  of  a storm  into  two  dimensions  (x,y). 


. v v 


6 


— 

Greene  noted  that  the  VIL  maxima  of  tornado  producing  storms  would 

. 

increase  to  high  values  approximately  one  hour  prior  to  the  occurrence 
of  a tornado,  and  then  increase  to  even  higher  values  at  the  time  of  a 
confirmed  tornado.  Greene  realized  moderate  success  with  this  tech- 
nique, and  his  early  efforts  led  to  additional  research  using  VIL  and 
variations  of  VIL. 

Morgan  and  Mueller  (1972)  studied  the  VIL  and  total  water  mass  of 
a large  Illinois  hailstorm,  using  photographically  derived  digital 
data  from  a 3 cm  radar.  They  discovered  that  hail-producing  storms 
were  characterized  by  very  large  values  of  VIL  and  rapid  increases  in 
the  total  water  mass  of  the  storm.  Clark  and  Can ip e (1972)  studied  the 
applicability  of  VIL  as  an  indicator  of  rainfall  rate  (R) , for  hydro- 
logical  purposes,  and  concluded  that  VIL  was  a better  indicator  of  R 
than  simple  Z-R  relationships.  Vogel  (1973)  used  VIL  and  CAZM  analyses 
to  study  several  severe  storms,  and  he  concluded  that:  1)  CAZMs  are 
the  best  method  with  which  to  analyze  severe  storms;  2)  analyses  of 
single-tilt  data  may  severely  limit  the  detection  of  severe  storms;  and 
3)  VIL  is  an  important  indicator  of  the  potential  for  severe  storms  to 
develop,  but  it  masks  some  of  the  three-dimensional  features,  such  as 
the  tilt  of  the  maximum  reflectivity  core  and  the  weak  echo  region 
(WER) . Canipe  (1973)  and  Canipe  and  Das  (1975)  developed  a technique 
for  partial  vertical  integration  of  liquid  water  content  (PVIL)  to  com- 
bine the  simplicity  and  speed  of  the  VIL  analysis  with  the  detail  of 
the  CAZM.  In  this  technique,  the  liquid  water  content  is  integrated 
over  three  layers  of  the  storm:  0-deg  to  15  kft  (the  source  region  of 
the  storm-no  freezing);  15  kft  to  35  kft  (the  active  region-freezing); 
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and  35  kft  to  50  kft  (the  ice  region-above  the  -40  deg  isotherm). 

Canipe  found  that  small  bounded  weak  echo  regions  (BWER)^  were  present 
through  the  lower  CAZM  analyses  of  tornadic  storms.  In  addition,  he 
concluded  that  PVIL  was  superior  to  V1L  because  it  depicted  both  the 
tilt  of  the  storm  system  and  the  loss  of  mass  in  the  upper  levels  of 
the  storm,  features  which  have  proven  to  be  valuable  indicators  of 
tornadic  activity. 

Elvander  (1975)  investigated  some  severe  storm  data  from  the 
spring  of  1972  provided  by  NSSL,  and  studied  the  relationship  between 
severe  storm  events  and  both  zero-tilt  indicators  (two-dimensional)  and 
VIL.  His  results  yielded  a correlation  coefficient  for  the  zero-tilt 
indicators  of  only  0.22,  while  that  for  the  VIL  was  0.45.  Elvander 
attributed  the  low  correlation  of  the  zero-tilt  indicators  to  the  fact 
that  they  do  not  indicate  the  intensity  of  a storm  throughout  its 
depth.  Correspondingly,  he  felt  that  the  correlation  of  VIL  was  much 
higher  due  to  the  fact  that  VIL  is  a three-dimensional  measurement. 

Pittman  (1976)  extended  Canipe's  work,  including  additional  data 
from  tornadic  storms.  Instead  of  computing  PVIL,  however,  Pittman 
summed  the  reflectivity  at  each  grid  point  in  the  layer  and  converted 
the  resultant  reflectivity  to  a dBZ  value.  This  produced  partial 
vertically-summed  reflectivity  (PVSZ)  maps.  He  found  that  the  BWERs 
associated  with  tornadoes  were  evident  in  the  lower  layer  PVSZ  maps  and 

^The  BWER  is  an  intense  updraft  within  which  the  entrained  water 
vapor  is  carried  aloft  so  rapidly  that  it  does  not  have  time  to  con- 
dense. Therefore,  the  updraft  has  very  few  water  droplets  to  back- 
scatter  radar  energy  and  appears  as  a weak  echo  region  bounded  by 
higher  reflectivity  values  (Canipe,  1973;  Sieland,  1977). 


that  the  reflectivity  contours  around  the  BWERs  exhibited  hook-shaped 
patterns.  This  feature,  together  with  the  tilt  of  the  storm  becoming 
either  more  nearly  vertical,  or  inclined  towards  the  BWER,  and  a de- 
crease of  the  upper  level  reflectivity  maximum,  correlated  well  with 
the  occurrence  of  tornadoes. 

Sieland  (1977)  developed  an  improved  computer  method  by  which 
digital  radar  data  are  reduced  into  an  optimum  number  of  PVSZ  layers. 
This  new  method  resulted  in  significant  savings  of  computer  processing 
time  and  memory,  as  compared  to  the  more  cumbersome  methods  of  Canipe 
or  Pittman,  while  at  the  same  time  it  retained  all  of  the  significant 
features  of  the  storm  complex  revealed  by  the  other  methods.  By 
improving  the  interpolation  techniques  employed,  Sieland' s analyses 
showed  greater  resolution  between  the  BWERs  and  the  reflectivity  con- 
tours surrounding  them.  His  technique  also  improved  the  detection  of 
marginal  BWERs,  where  Pittman's  technique  had  a tendency  to  smooth 
over  the  data. 

Radlein  (1977)  further  refined  Sieland 's  computer  program  to  the 
extent  that  it  would  process  data  from  the  TAMU  digital  weather  radar 
system  in  the  1 deg  by  1 km  mode.  She  then  examined  simultaneous,  dual 
wavelength  data  from  a Texas  thunderstorm,  in  an  effort  to  determine 


whether  any  of  the  radar  signatures  associated  with  tornado-producing 
thunderstorms  elsewhere  were  present  in  this  central  Texas  storm. 

The  Basis  for  the  Investigation 


With  the  initiation  of  the  Digitized  Radar  Experiment  (D/RADEX) , 
the  need  for  an  operational  technique  to  process  digital  radar  data 
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became  evident.  Since  the  NWS,  and  possibly  the  U.S.  Air  Force  Air 
Weather  Service  (AWS) , will  have  many  of  their  radar  systems  digitized 
in  the  near  future,  data  reduction  schemes  that  can  be  used  in  real 
time  are  a necessity.  Any  such  system  must  allow  for  easy  analysis  of 
the  three-dimensional  aspects  of  a storm,  since  zero-tilt  reflectivity 
patterns  alone  have  proved  to  be  inadequate  as  indicators  of  severe 
storms.  It  is  also  imperative  that  any  data  reduction,  interpolation, 
and  presentation  schemes  preserve  the  significant  details  of  the  indi- 
vidual CAZM  analyses  and  the  three-dimensional  characteristics  of 
the  storm.  It  is  very  possible  that  the  retention  of  significant 
details  and  characteristic  features  of  a storm  is  highly  dependent 
upon  the  resolution  of  the  data  and,  correspondingly,  upon  the  mode  of 
data  collection. 


Techniques  of  the  Investigation 

The  data  for  this  investigation  were  obtained  from  the  dual 
wavelength  digital  weather  radar  system  in  the  Department  of 
Meteorology  at  Texas  A&M  University.  The  reduction,  interpolation, 
and  display  of  the  data  were  accomplished  by  the  Amdahl  470V/6  com- 
puter at  TAMU  using  the  technique  developed  by  Sieland. 

Zero-tilt  data  were  collected  simultaneously  at  both  3 and  10  cm 
wavelengths,  and  near  simultaneously  in  all  four  modes  of  data  col- 
lection (1  deg  by  1 km,  1 deg  by  2 km,  2 deg  by  1 km,  and  2 deg  by 
2 km) , for  targets  at  both  near  range  (approximately  50  km)  and  far 
range  (greater  than  100  km) . The  data  were  processed  and  contoured 
reflectivity  maps  generated  for  the  zero-tilt  display.  These  maps 
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were  analyzed  with  respect  to  one  another  for  such  characteristics  as 
resolution,  retention  of  significant  echo  features,  smoothing  of  fine 
scale  detail,  and  differences  due  to  range  of  the  storm  and  radar 
wavelength.  Additionally,  these  factors  were  weighed  against  the 
amount  of  computer  resources,  both  hardware  and  software,  required  to 
process  the  data.  The  ultimate  objective  was  to  determine  an  opera- 
tionally optimum  mode  of  data  collection  based  on  all  of  these  con- 
siderations. 
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CHAPTER  II 

THE  TEXAS  A&M  UNIVERSITY  WEATHER  RADAR  SYSTEM 

Basic  Radar  Theory 

The  form  of  the  radar  equation  used  in  this  investigation  follows 
the  derivation  of  Probert-Jones  (1962) , who  assumed  a more  realistic 
beam  shape  in  order  to  reduce  the  error  in  a derivation  by  Battan 
(1959).  In  the  Probert-Jones  derivation,  we  assume  that  there  is  no 
significant  attenuation  of  microwave  energy  between  the  radar  antenna 
and  the  target,  and  that  the  target  completely  fills  the  spatial 
volume  illuminated  by  the  radar  beam.  Under  these  conditions,  the 
average  backscattered  power,  Pr  (watts),  received  from  the  target  at 
range  r(km),  is  given  by 


P 

r 


Z 

e 


(1) 


2 3 —6 

where  C (w  km  m mm  ) is  the  radar  constant,  which  is  unique  for  each 

2 

individual  radar,  |k|  (dimensionless)  is  the  dielectric  constant  used 

6 —3 

in  scattering  theory,  and  (mm  m ) is  the  equivalent  radar  reflec- 
tivity factor.  The  equation  for  C is 
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it  pt  g he<j> 
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(2) 


where  P is  the  transmitted  power  during  a radar  pulse  (w) , G is  the 
antenna  gain  (dimensionless) , h is  the  pulse  length  in  space  of  the 
transmitted  pulse  (cm),  0 is  the  horizontal  beamwidth  (radians),  (f>  is 


transmitted  microwave  energy  (cm) . The  factor  2 In  2 is  the  Probert- 
Jones  correction  factor,  which  results  from  the  fact  that  the  beam 
power  density  nearly  approximates  a Gaussian  distribution;  conse- 
quently, the  transmitting  and  receiving  gains  are  not  equal.  The 

value  of  C for  each  of  the  TAMU  radars  is  found  below,  and  the  value 
2 

of  |k|  is  chosen  as  0.93  (Bat tan,  1973). 

The  final  output  of  the  radar  system,  which  is  digitized  values 
of  the  returned  power,  P , must  be  converted  to  values  of  reflectivity 
before  any  further  data  reduction  techniques  can  be  applied.  The  usual 
practice  is  to  measure  the  returned  power  in  terms  of  decibels  with 
respect  to  a standard  reference  power  level,  normally  1 raw.  Power 
levels  are  then  expressed  in  units  of  dBm,  either  above  (+)  or  below 
(-)  1 mw.  The  equation  is 


P(dBm)  = 10  log 


P (watts) 

10  TO"3 

10  watts 


(3) 


For  example,  a returned  power  (P^_)  of  10  ^ watts  would  yield  a value 
2 

of  -80  dBm  for  P. 

To  determine  values  of  Z , we  solve  Eq  (1)  for  Z and  have 

e n e 


2 

— ^~~2  P 

c|k|  r 


(4) 


^See  Battan  (1973)  for  a more  complete  discussion  of  radar  reflec- 
tivity factor,  Z,  equivalent  radar  reflectivity  factor,  Z , and  dBm. 
Also,  in  this  investigation,  In  refers  to  base  e logarithms,  and  log 
to  base  10  logarithms. 


The  logarithm  of  Eq  (4)  is 


log  Ze  = 2 log  r + log  - log  C | k| ' 


Using  the  values  of  C for  the  TAMU  radars  given  below  and  that  of 


k|  given  earlier,  we  have 


C3  = 1.0089  x 10 
C1Q  = 8.609  x 10" 


log  C3 |k| 2 = 9.0 
log  C1Q|k|2  = 10.1 


By  substituting  Eq  (6)  into  Eq  (5)  we  get 


log  Ze3  = 2 log  r + log  + 9.0 


log  Z = 2 log  r + log  P._  + 10.1 


The  digital  value  of  P^,  is  converted  to  its  dBm  equivalent  (always 

negative)  through  the  use  of  calibration  data  for  each  of  the  two 
3 

radars.  The  received  power  in  watts  is  then  related  to  the  received 
power  in  dBm  by 


log  P^_  (watts)  = 0.1  P^  (dBm)  - 3.0 


Substitution  of  Eq  (8)  into  Eq  (7)  yields 


log  Ze3  = 0.1  Pr  + 2 log  r + 6.0 


log  ZelQ  = 0.1  Pf  + 2 log  r + 7.1 


Finally,  the  value  of  Zg  is  given  by 


The  digital  system  produces  log  P rather  than  log  P . Wilk  and 
Kessler  (1970)  developed  an  equation  which  corrects  Eq  (7)*  but  the 
correction  is  usually  less  than  1 dB. 
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Z = 10(0'1  Pr  + 2 log  r + wc) 


6 -3 

where,  by  convention,  Z£  is  in  units  of  mm  m , and  is  6.0  for  the 
3 cm  radar  and  7.1  for  the  10  cm  radar. 


Assumptions  Made  in  the  Theoretical  Development 


There  are  several  assumptions  that  are  inherent  in  the  development 


of  Eq  (10).  These  include: 

1.  The  transmitted  microwave  energy  is  not  significantly 

attenuated  between  the  radar  antenna  and  the  target.  This  has  been 

shown  to  be  valid  for  radars  with  wavelengths  near  10  cm,  but  is  not 

strictly  true  for  wavelengths  in  the  vicinity  of  3 cm  (Greene,  1964). 

The  returned  power  of  a 3 cm  radar  will  suffer  some  reduction  due  to 

attenuation  by  atmospheric  gases,  cloud  droplets,  and  precipitation, 

-2 

The  attenuation  due  to  atmospheric  gases  will  be  on  the  order  of  10 

dB  km  \ while  that  due  to  cloud  droplets,  depending  on  their  physical 

-4  -1  -1 

state  and  temperature,  will  vary  between  10  and  10  dB  km  . At- 
tenuation due  to  precipitation  is  generally  proportional  to  rainfall 
rate  and  temperature,  and  from  derived  Z-R  relationships,  can  be  ap- 
proximately expressed  as  2.9  x 10  ^ dB  km  ^ (Battan,  1973). 

2.  The  spatial  volume  illuminated  by  the  radar  beam  is  completely 
filled  by  the  target.  This  condition  is  not  always  met,  particularly 
on  the  periphery  of  a storm.  The  result  will  be  a loss  of  resolution 
of  fine  scale  detail  in  these  areas  (Greene,  1971),  but  it  will  not 
significantly  affect  the  major  features  present  in  the  digital  data. 

3.  The  Rayleigh  approximation  is  used  to  describe  the  scattering 
properties  of  spherical  liquid  water  drops  having  diameters  on  the 


■ 
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order  of  0.04  X,  where  X is  the  radar  wavelength  in  centimeters.  In 
the  case  of  severe  storms,  large,  non-spherical  water  drops  and  hail- 
stones may  be  present  that  do  not  meet  the  Rayleigh  criteria;  however, 
the  enhanced  reflectivity  that  is  characteristic  of  these  particles  may 
i be  useful  in  identifying  severe  storms. 

Two  additional  assumptions  that  are  made  concerning  the  digital 
data  in  this  investigation  are: 

1.  Each  digital  datum  value  represents  a point  in  the  center  of 
the  radar  volume. 

2.  The  equivalent  radar  reflectivity  factors,  Z fi,  obtained  from 
the  digital  data,  are  representative  of  a continuous  scalar  field. 


The  Earth  Curvature  Correction 

Microwave  radiation  propagating  in  free  space  will  follow  a 
straightline  path.  However,  the  microwave  radiation  of  a radar  beam 
propagating  through  the  atmosphere  will  follow  a curved  path,  due  to 
refraction.  The  amount  of  curvature  of  the  path  depends  upon  the 
vertical  gradient  of  the  index  of  refraction  (n) . Under  normal  at- 
mospheric conditions,  where  temperature  and  humidity  decrease  with 
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Ray  theory  may  be  applied  to  the  problem  if  the  change  with  height 

of  the  index  of  refraction,  44»  is  small.  In  that  case,  the  exact  dif- 

dh 

ferential  equation  for  a ray  in  a spherically  stratified  atmosphere  is 


d2h  _ / 2 1 dnV  dhV  / R+h\2  / 1 1 dn  \ 

ds2  y R+h  n dh/Vds/  \ R ) \R+h  n dh  / 


(ID 


where  h is  the  height  of  the  beam  above  the  earth's  surface  at  a 

distance  s from  the  transmitter,  R is  the  earth's  radius,  and  n is  the 

2 

index  of  refraction.  Since  e is  usually  very  small,  = tan2  e<<l. 

Additionally,  nsl,  and  h<<R,  so  for  all  practical  purposes,  Eq  (11) 
can  be  reduced  to 


d2h  = 1 dn 

,2  R dh 
ds 


(12) 


In  this  investigation  it  is  necessary  to  consider  the  radar  beam  axis 
as  a straight  line,  so  we  will  assume  a fictitious  earth  having  a 
radius  given  by  (Battan,  1959) 


R' 


1 + R 


dn 

dh 


(13) 


It  has  been  shown  that  44  is  small  and  nearly  linear  with  a value  of 

dh 

—8  — 1 4 

-4  x 10  m . Therefore,  an  earth  curvature  correction  of  R1  = -j  R 

will  be  used,  and  the  resultant  radar  beam  may  be  considered  a straight 

line,  as  shown  in  Figure  2. 
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The  TAMU  Weather  Radar  System 

Physical  Description 

The  TAMU  Weather  Radar  System  consists  of  two  complete  weather 
radar  systems  operating  in  parallel.  The  makeup  of  the  two  systems  is 
nearly  identical,  and  each  consists  of  the  following  major  sub-systems: 
antenna,  conventional  analog  weather  radar  set,  digital  video 
integrator  and  processor  (DVIP)  with  PPI  display,  and  a nine-track 
tape  unit  to  record  the  buffered  digital  output  of  the  DVIP.  Each  of 
these  subsystems  will  be  addressed  separately.  The  technical  charac- 
teristics of  both  radar  sets  and  the  antenna  system  are  contained  in 
Table  1. 

Antennas.  The  antenna  system  consists  of  an  array  of  three  solid 
parabolic  dish  antennas  that  are  mounted  on  a single  pedestal  and 
permanently  aligned  in  relation  to  one  another.  (The  purpose  of  the 
third  antenna  is  to  allow  the  possible  installation  of  a K band  radar 
at  some  future  date.)  The  large  S band  antenna  (10  cm  radar)  is 
mounted  in  the  center,  with  the  two  smaller  X (3  cm)  and  K band  anten- 
nas mounted  on  either  side. 

Radar  sets.  The  3-cm  radar  is  basically  a modified  CPS-9  radar, 
while  the  10-cm  radar  was  constructed  at  TAMU  using  components  taken 
from  a variety  of  different  radar  sets.  Both  units  do  have  new  re- 
ceivers, however,  and  both  are  mounted  in  CPS-9  consoles.  The  system 
has  only  a single  synchronizing  pulse  generator  that  is  used  to  trigger 
both  radars.  This,  combined  with  the  single  antenna  pedestal,  insures 
that  both  radars  scan  the  same  illuminated  volume  simultaneously.  The 
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TABLE  1.  WSR/TAM-2  Weather  Radar  Technical  Characteristics. 
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synchronizing  trigger  signal  and  the  radar  receiver  logarithmic  video 
output  comprise  the  input  into  the  next  major  subsystem,  the  DVIP. 

The  Digital  Video  Integrator  and  Processor  (DVIP) 

The  DVIP  is  a high  speed  data  acquisition,  digitizing,  and 

processing  system,  which  continuously  averages  radar  logarithmic  video 

output  in  range  (range  averaging)  and  in  direction  of  antenna  scan 

(time  averaging),  using  exponentially  weighted  digital  integrator 

4 

techniques  and  synchronized  by  the  radar  system  trigger.  The  results 
are  quantitative  estimates  of  mean  detected  precipitation  returns  and 
an  improvement  in  the  accuracy  (standard  deviation)  of  the  mean  in- 
tensity estimates,  over  that  of  a logarithmic  radar  receiver,  to  1.0  dB 
or  less,  based  on  range  increments  of  1 or  2 km.  The  DVIP  installation 
provides  for  real  time  contoured  PPI  displays  and  digital  recording  of 
the  digitized,  integrated  data.  Both  radars  are  coupled  to  identical 
DVIP  systems,  except  that  the  3 cm  system  also  has  a tape  playback 
facility  incorporated  into  it. 

DVIP  accuracy.  The  DVIP  can  improve  the  accuracy,  measured  in 
terms  of  standard  deviation,  of  precipitation  intensity  estimates  to 
1.0  dB  or  less,  at  range  increments  of  1 or  2 kilometers.  This  im- 
provement  is  accomplished  by  a combination  of  range  averaging  (an 
average  over  several  sample  volumes  in  range)  and  an  average  over 
several  pulses  returned  from  the  same  volume  (continuous  pulse-to- 
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The  accuracy,  or  standard  deviation  (a),  of  the  estimate  of  the 
mean  intensity  values  of  radar  returns  from  a precipitation  volume  is 
a funrt-i^n  of  two  things,  the  effective  sample  size,  and  the  fluctua- 
tion of  the  output  of  a logarithmic  receiver.  The  latter,  fluctuating 


characteristics  of  precipitation  returns  from  a log  receiver,  are  de- 
fined by  a standard  deviation  of  5.57  dB.  Accordingly  then,  if  a 
number  of  statistically  independent  samples  (N)  from  a particular  pre- 
cipitation volume  (the  output  of  the  log  receiver)  are  averaged  by  the 
DVIP,  the  standard  deviation  of  the  intensity  estimates  is  reduced  by 


5.57  dB 

✓S’ 


(14) 


To  determine  the  number  of  independent  samples  to  be  integrated  by  the 
DVIP,  we  assume  that  the  radar  return  signal  from  a sample  precipita- 
tion volume,  V , varies  randomly  with  time  and  has  a Gaussian  distri- 
bution. The  total  number  of  samples,  N , obtained  from  the  volume,  V , 

s s 

is  dependent  upon  the  mechanics  of  the  radar  set  itself,  specifically 
antenna  beamwidth,  PRF,  and  antenna  scan  rate.  Because  there  may  be 
some  sample-to-sample  correlation,  the  number  of  independent  samples, 

N,  in  the  total  sample,  Ng,  must  be  determined.  The  degree  of  corre- 
lation between  the  samples  is  dependent  upon  both  radar  and  meteorologi- 
cal parameters,  the  latter  being  the  standard  deviation  of  the  velocity 

distribution  of  the  scatterers,  a , in  the  samnle  volume.  The  number 

v 

of  independent  samples,  N,  is  related  to  the  total  number  of  samples. 


N , by 
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N = 


4/i 7 a T N 
v s s 

X 


(15) 


where  Tg  is  the  sampling  time  (1/PRF) , and  X is  the  transmitted  wave- 
length. If  we  set  N = Ng,  assume  ov  = 2 m/s,  and  solve  for  Ts>  we  find 
that  the  PRF  required  for  independent  sampling  with  a 3 cm  (X  band) 
radar  is 


T 

s 


X 

4/F  a 

v 


.032 
4/tF  2 


.002257  sec 


(16) 


PRF 


= 443  pps 


(17) 


Therefore,  if  the  radar  PRF  is  443  pps  or  less,  all  pulse-to-pulse 
samples,  or  time  samples,  of  the  volume,  V , will  be  independent. 

The  number  of  statistically  independent  samples,  N,  that  the  DVIP 
must  average  to  meet  the  standard  deviation  requirement  of  1.0  dB  or 
less  is 


(5.57  dB)2 
2 


(18) 


o 

so  for  a = 1.0,  the  DVIP  must  process  31  independent  samples. 

The  pulse  length  of  the  TAMU  radars  is  5 microseconds,  while  the 
specified  1 km  range  increment  is  6.67  microseconds  long,  so  each 
range  increment  contains  1.34  independent  range  samples.  If  at  least 
31  independent  samples  must  be  processed  by  the  DVIP,  then  the  number 
of  time  samples,  or  pulses,  that  must  be  processed  in  the  digital 
integrator  is  23.13  (31/1.34  = 23.13).  For  the  case  where  PRF  = 186 
with  an  antenna  scan  rate  of  1.3  rpm,  we  get  23.85  pulses  per  degree  of 
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azimuth  for  a 1 km  range  increment.  The  2 km  range  increment  contains 
2.67  independent  range  samples,  so  by  solving  Eq  (18)  for  a,  and  sub- 
stituting these  numbers  (23.85  and  2.67)  into  N,  we  find  that,  for  the 
2 km  range  increment,  a = 0.698  dB.  On  the  other  hand,  the  maximum 
allowable  antenna  scan  rate  that  will  insure  an  accuracy  of  1.0  dB  is 
increased  to  2.67  rpm. 

Integration  technique.  The  only  inputs  to  the  DVIP  are  115  VAC 
power,  the  radar  logarithmic  video,  and  the  radar  system  trigger.  The 
range  increment  integration,  over  either  1 or  2 km,  and  the  number  of 
time  samples  (pulses)  to  be  integrated,  either  15  or  31,  are  controlled 
by  switches  on  the  front  panel  of  the  unit.  The  azimuth  increment  to 
be  integrated,  either  1 or  2 deg,  is  controlled  by  the  tape  unit  and 
will  be  covered  separately  later. 

The  incoming  log  video  signal  from  the  radar  set  is  first  ampli- 
fied by  a video  scaling  amplifier,  and  then  the  analog  signal  is 
sampled  once  every  1.66  microseconds  (about  every  250  meters  in  range) 
and  the  samples  fed  into  an  analog-to-digital  converter.  The  A/D  con- 
version takes  about  one  microsecond  and  converts  the  analog  video 
signal  to  digital  bits.  From  the  A/D  converter,  the  digital  data 
samples  are  sent  to  a range  bin  circulator.  Here  the  samples  are  added 
to  one  another,  either  four  or  eight  of  them,  depending  on  the  range 
increment  selected,  and  then  the  sum  is  divided  by  either  four  or  eight, 
resulting  in  an  average  value  for  that  range  increment.  The  next  step 
is  through  a multiplexer  and  into  range  bin  storage.  It  is  during  this 
step  that  the  time  sample  integration  occurs,  controlled  by  the  time  sam- 
ple switch  on  the  front  panel.  This  switch  provides  a choice  of  15  or 31 
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steps  in  the  integration  curve  between  the  10  and  90  percent  amplitude 
points,  where  each  step  covers  one  sample  time.  The  integrator  equa- 
tion forms  a mathematical  progression  where  one-eighth  of  a number  added 
to  seven-eighths  of  one-eighth  of  that  number  will  arrive  at  90  percent 
of  the  original  number  in  15  steps.  Where  one-sixteenth  is  chosen 
instead  of  one-eighth,  the  number  of  steps  between  the  10  and  90  per- 
cent amplitude  points  is  31.  The  maximum  antenna  scan  rates  that  will 
allow  15  and  31  pulses  per  degree  of  azimuth  are  2.07  rpm  and  1.00  rpm, 
respectively. 

It  should  be  noted  that  this  time  sample  integration  step,  which 
requires  15  or  31  separate  radar  pulses,  is  not  to  be  confused  with  the 
range  increment  integration,  which  requires  31  statistically  inde- 
pendent samples,  which  can  be  obtained  from  23.13  time  samples  (pulses). 
Likewise,  neither  of  these  processes  should  be  confused  with  the  250  m 
(1.66  y sec)  sampling  and  averaging  used  to  produce  mean  values  for 
the  1 and  2 km  range  increment  received  power  values,  or  the  azimuth 
increment  integration,  which  is  controlled  by  the  tape  unit. 

DVIP  output.  The  DVIP  output  signal  is  split  into  two  channels, 
the  digital  channel  and  the  display  channel.  The  display  channel  pro- 
vides either  log  video  or  contoured  log  video  to  the  system  PPI  scope. 
The  DVIP  can  be  calibrated  to  detect  and  display  six  discrete  levels 
of  precipitation  intensity.  The  presentation  of  these  six  intensity 
levels  as  three  alternating  shades  of  grey  is  called  video  contouring. 
These  threshold  levels  can  be  set  digitally  in  0.5  dR  steps.  A radar 
return  that  falls  within  level  1 is  displayed  as  grey  (half  bright- 
ness) , level  2 is  displayed  as  white  (full  brightness) , and  level  3 is 
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black  (zero  brightness).  The  cycle  then  repeats  itself;  level  4 is 
grey,  level  5 is  white,  and  level  6,  the  highest  intensity  return,  is 
black. 

The  digital  output  consists  of  the  1 or  2 km  range  increment 
integrated  video  samples,  provided  in  a buffered  8 bit  parallel  binary 
configuration  compatible  with  standard  TTL  logic  levels.  Each  digital 
output  word  represents  a digital  value  of  the  integrated  video  in- 
tensity on  a linear  scale  over  the  input  dynamic  range  of  80  dB.  At 
the  digital  output  there  are  215  2-km  or  430  1-km  digital  samples 
representing  the  range  increments  from  21  through  450  km  of  range, 
depending  upon  the  range  increment  selected. 

Range  normalization  (STC) . The  range  normalization  function 

(sensitivity  time  control-STC)  compensates  the  integrated  video  signal 

2 

for  attenuation  due  to  range.  It  provides  a 1/R  correction  (R  = radar 
range)  in  the  range  interval  from  20  to  230  km,  with  0 dB  correction  at 
20  km  and  +21  dB  correction  at  230  km,  applied  in  increments  of  0.5  dB. 
Range  normalization  is  applied  only  to  video  signals  of  precipitation 
intensity  level  2 or  higher.  Range  normalization  is  not  applied  to  the 
digital  output  channel  at  all,  and  can  be  disabled  for  the  display 
channel  as  well,  by  means  of  a front  panel  switch. 

Range  blanking.  Two  range  blanking  options  are  available,  a 
fixed  0-20  km  interval,  and  a variable  range  blanking  interval  adjust- 
able from  0-199  km.  Within  these  intervals,  all  contoured  log  video  to 
the  radar  displays  is  inhibited.  The  type  of  range  blanking  interval 
employed  is  selected  by  a front  panel  switch,  the  operation  of  which 
does  not  affect  the  digital  output  data,  since  these  data  are  blanked 
in  the  0-20  km  interval  only. 
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Precipitation  level  selection.  Six  front  panel  mounted  switches 
are  provided  to  select  any  one  or  more  of  the  six  precipitation  in- 
tensity levels  for  display.  Selection  of  one  or  more  of  the  levels 
disables  normal  contouring,  and  only  those  levels  selected  are  dis- 
played. All  selected  levels  are  displayed  as  white  (full  brightness). 

For  example,  if  level  3 is  selected,  all  precipitation  areas  with  in- 
tensity between  the  level  3 threshold  and  the  level  4 threshold  are 
displayed  as  white,  and  no  other  precipitation  areas  or  intensities 
are  displayed. 

Digital  Output  Recording  System 

The  digital  output  channel  of  the  DVIP  is  fed  into  a 9-track 
magnetic  tape  recording  unit.  This  tape  drive  unit  serves  two  func- 
tions; it  records  the  digital  data  in  a format  that  makes  it  readily 
accessible  by  computer,  and  it  controls  the  azimuth  increment  inte- 
gration of  the  data.  The  first  function  will  be  addressed  further  in 
the  section  on  data  reduction. 

Control  of  the  azimuth  integration  increment  (data  integrated 
over  either  1 deg  or  2 deg  of  azimuth)  is  accomplished  through  the  use 
of  a switch  on  the  front  panel  of  the  tape  unit.  At  the  end  of  each 
azimuth  increment,  the  tape  unit  generates  a signal  and  calls  for  a 
digital  value  from  the  DVIP  that  is  based  upon  the  most  recent  time 
samples  from  a particular  range  bin.  The  number  of  samples  called  for 
is  determined  by  the  time  sample  switch  on  the  DVIP,  which  selects 
either  15  or  31  samples.  The  digital  value  derived  from  those  samples 
represents  the  mean  received  power  from  that  entire  azimuth  increment. 
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Consequently,  it  becomes  very  important  to  match  the  antenna  scan  rate 
with  the  azimuth  increment  selected.  For  example,  if  the  time  sample 
selected  is  31  and  the  azimuth  increment  selected  is  .1  deg,  the  value 
recorded  for  that  range  bin  is  derived  from  the  last  31  time  samples 
over  that  1 deg  of  azimuth.  At  1.0  rpm  there  will  just  be  31  samples 
in  the  1 deg  of  azimuth.  If  the  azimuth  increment  selected  is  2 deg, 
the  digital  value  recorded  is  still  derived  from  the  last  31  samples. 
However,  if  the  antenna  scan  rate  is  still  1.0  rpm,  those  31  samples 
will  all  have  come  from  the  second  half  of  the  2 deg  azimuth  increment, 
and  consequently  may  not  represent  that  range  bin  accurately.  There- 
fore, with  an  azimuth  increment  of  2 deg  it  is  necessary  to  increase 
the  antenna  scan  rate  to  2.0  rpm.  Conversely,  it  is  also  important 
that  the  antenna  scan  rate  not  be  too  fast,  since,  in  that  case,  the 
recorded  values  are  derived  from  an  integration  scheme  which  requires  a 
fixed  number  of  samples,  but  does  not  get  them.  Once  again,  the  re-, 
corded  value  may  not  be  representative  of  the  actual  value. 

Data  collection  modes.  There  are  four  possible  modes  of  data 
collection  with  this  system.  They  are: 

1.  1 deg  by  1 km.  The  data  are  integrated  over  1 km  range 
increments  by  the  DVIP  and  the  integrated  time  sample  data  called  for 
by  the  tape  unit  represent  1 deg  of  azimuth, 

2.  1 deg  by  2 km.  The  data  are  integrated  over  2 km  range 
increments  by  the  DVIP  and  the  integrated  time  sample  data  called  for 
by  the  tape  unit  represent  1 deg  of  azimuth, 

3.  2 deg  by  1 km.  The  data  are  integrated  over  1 km  range 
increments  by  the  DVIP  and  the  integrated  time  sample  data  called  for 
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by  the  tape  unit  represent  2 deg  of  azimuth. 

4.  2 deg  by  2 km.  The  data  are  integrated  over  2 km  range 
increments  by  the  DVIP  and  the  integrated  time  sample  data  called  for 
by  the  tape  unit  represent  2 deg  of  azimuth. 

In  all  cases  the  data  are  recorded  one  azimuth  increment  at  a 
time,  beginning  with  the  range  bin  nearest  the  radar,  and  proceeding 
sequentially  outwards  along  the  radial. 

Table  2 presents  the  antenna  scan  rates  required  to  maintain  an 
accuracy  of  1.0  dB  in  all  four  collection  modes,  for  both  15  and  31 
time  samples. 


TABLE  2.  Antenna  Scan  Rates 

(rpm)  Required  for 

1.0  dB  Accuracy. 

Collection  Mode 

15  Time  Samples 

31  Time  Samples 

1 deg  x 1 km 

2.07 

1.00 

1 deg  x 2 km 

2.07 

1.00 

2 deg  x 1 km 

4.13 

2.00 

2 deg  x 2 km 

4.13 

2.00 

Data  shift.  Both  output  channels,  the  coded  PPI  display  channel 
and  the  digital  channel,  suffer  from  a slight  shifting  of  the  displayed 
location  of  the  data  points,  both  in  range  and  in  azimuth.  This  is 
due  to  the  nature  of  the  integration  processes  employed  by  the  DVIP, 
and  the  magnitude  of  the  effect  is  a function  of  the  range  and  azimuth 
increments  selected.  In  both  cases,  the  shift  occurs  because  the  in- 
coming reflectivity  values  are  integrated  across  the  entire  range  or 
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azimuth  increment,  from  its  leading  edge  to  its  trailing  edge.  Then 
the  integrated  value  is  displayed  at  the  end  of  the  increment,  rather 
than  at  the  mid-point  of  the  increment  over  which  the  data  were  inte- 
grated. In  the  contoured  PPI  display,  this  data  shift  manifests  itself 
in  several  ways.  First,  the  leading  edges  of  all  targets  at  all 
threshold  levels  are  shifted  outward  and  clockwise  by  one  range  inter- 
val and  one  azimuth  interval,  respectively.  The  magnitude  of  the 
shifts  depends  upon  whether  the  1 deg,  2 deg,  1 km,  or  2 km  modes  are 
being  utilized.  Additionally,  the  trailing  edges  of  all  targets  are 
shifted  outwards  1 to  2 range  increments,  depending  upon  the  amplitude 
of  the  target's  reflectivity  at  the  range  increment  boundary. 

The  data  shift  in  the  digital  output  channel  is  partially  overcome 
by  the  use  of  the  assumption  that  the  data  represent  points  in  the 
center  of  the  illuminated  radar  volume.  This  negates  the  range  shift, 
but  only  partially  counteracts  the  azimuth  shift.  This  is  because  the 
data  points  are  assumed  to  lie  on  the  radial  at  the  end  of  the  azimuth 
increment,  rather  than  at  the  angular  mid-point  of  that  increment. 

Even  so,  this  reduces  the  shift  to  only  one  half  of  one  azimuth  incre- 
ment . 

The  reduction  and  display  of  the  digital  data  output  of  the  DVIP 
are  addressed  in  the  next  chapter. 
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CHAPTER  III 

DATA  REDUCTION  AND  DISPLAY 

The  digital  output  channel  of  the  DVIP  is  recorded  on  9-track 
tape,  and  this  tape  then  becomes  the  input  into  the  TAMU  Amdahl  470V/6 
computer.  Through  the  use  of  the  program  developed  by  Sieland  (1977) 
and  a library  subroutine  called  CONREC,  the  computer  processes  the 
digital  radar  data,  and  presents  it  to  the  user  in  a functional  format. 
Both  of  these  computer  programs  are  described  and  explained  in  the 
following  sections. 


Sieland 's  RADAR  Program 

The  computer  program  developed  by  Sieland,  called  RADAR,  begins 
with  the  digital  radar  data,  which  is  recorded  in  a spherical  coordinate 
system  (r,a,e)  and  first  transforms  it  into  a cylindrical  coordinate 
system  (r,a,h)  with  three  separate  height  divisions,  or  PVSZ  layers. 
Finally,  the  data  are  transformed  into  a two-dimensional  rectangular 
coordinate  system  (x,y)  and  five  maps  are  produced,  a zero  level  re- 
flectivity map,  the  three  PVSZ  layer  maps,  and  a VIL  map.  The  con- 
version of  coordinate  systems  requires  the  use  of  two  interpolation 
schemes,  a Lagrangian  linear  or  cubic  interpolation  scheme  along  each 
radial  of  data  (in  the  r-e  plane)  and  a quadratic  scheme  for  the  con- 
version to  rectangular  coordinates.  An  extensive  explanation  of  the 
development  and  functioning  of  these  interpolation  schemes  is  contained 
in  Sieland  (1977),  the  main  points  of  which  are  presented  here. 
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PVSZ  Layers 

This  program  creates  three  PVSZ  layers,  where  the  reflectivity 
factors  are  summed  in  a vertical  column  above  a given  (x,y)  point  at 
the  surface  into  the  lower  layer  (0-deg  to  15  kft),  the  middle  layer 
(15  to  35  kft),  and  the  upper  layer  (35  to  50  kft)  of  a storm  complex. 
The  first  step  in  this  procedure  is  to  convert  the  data  from  a spheri- 
cal coordinate  system  (r,a,e)  to  a cylindrical  one  (r,a,h).  Figure  3 
shows  the  relationship  between  these  two  coordinate  systems,  and  the 
interpolation  points  used  to  convert  from  one  to  another.  (Each  dot 
at  the  intersection  of  a column  and  a radial  is  an  interpolation 
point.)  Note  that  the  number  of  interpolation  points  in  each  layer 
varies  with  the  range,  r,  so  it  becomes  necessary  to  determine  the 
number  of  interpolation  points  (and  hence,  elevation  angles  of  data) 
contained  in  each  layer  above  each  range  bin.^  This  is  accomplished  by 
a subroutine  called  DIVIDE,  which  calculates  the  number  of  tilt  angles 
contained  in  each  of  the  three  layers  at  each  range  and  stores  them  in 
the  arrays  IDIV1,  IDIV2,  and  IDIV3,  for  the  lower,  middle,  and  upper 
layers. 

Conversion  to  Cylindrical  Coordinate  System 

The  Lagrangia  i linear  or  cubic  interpolation  scheme  incorporated 
by  Sieland  to  convert  the  data  from  a spherical  to  a cylindrical 
coordinate  system  is  applied  to  the  data  along  each  radial  of  informa- 
tion at  each  separate  elevation  angle.  The  steps  involved  in  this 
interpolation  follow  (See  Figure  4) : 
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1.  Proceed  sequentially  in  range  along  a given  radial  (a^)  until 
the  first  non-zero  reflectivity  factor  is  found. 

2.  If  the  reflectivity  factors  at  the  given  range,  r^,  and  at 
r^+^  are  equal,  no  interpolation  is  required. 

3.  If  the  reflectivity  factors  are  unequal,  then  calculate  the 
radial  distance,  r',  to  the  point  on  the  data  radial  that  intersects 
the  vertical  column  extending  above  the  applicable  range  gate,  r^,  at 
the  surface.  Here  h is  given  by 

2 
r . 

h = r±  tanei  + (19) 

where  R'  = and  r'  is  given  by 

r'  = /r2  + h2  (20) 

Equation  (19)  includes  corrections  for  both  Earth's  curvature  and 
refraction  of  the  radar  beam  under  standard  atmospheric  conditions. 

4.  Compute  the  distance  from  r^  to  r' 


fir  = r' 


- r. 


(21) 


If  fir  < 1,  the  four  data  points  used  in  the  Lagrangian  cubic  interpo- 
lation formula  are  r.  , , r . , r. and  r.,„.  If  fir  > 1,  i is  increased 

i-l  i i+l  i+2  — 

by  1 and  fir  is  decreased  by  1;  if  fir  j>2,  i is  increased  by  2 and  fir  is 
decreased  by  2;  etc. 

5.  Interpolate  the  value  of  the  reflectivity  factor  at  the  point 


TSSSf 


(r',h)  by  using  the  Lagrangian  cubic  form  of  the  interpolating  poly- 
nomial given,  in  general  terms,  by 


36 


where 


(r')  = 


E 

k=0 


^i-l+k^ 


K U’) 


(22) 


Lk  = w 

j*k 


r - r 


i-l±l 


ri-l+k  ri-l+j 


(23) 


The  form  of  the  Lagrangian  interpolating  polynomial  used  in  this  pro- 
gram is  given  by 


Ze  (r')  = 

a 

o 

(r’ 

- r.)(r'  - 

ri+l)(r'  - 

ri+2) 

+ 

ai 

(r* 

- ri-l><r’ 

- ri+l)(r' 

ri+2^ 

+ 

a2 

(r ' 

- ri-l>(r’ 

- r.)(r'  - 

ri+2) 

+ 

a3 

(r ' 

- ri-l)(r’ 

- r1)(r'  - 

ri+l) 

(24) 

where 


a = 


Ze  <ri-l> 


(ri-l  ri)(ri-l  ~ ri+l)(ri-l  ” ri+2) 


a,  = 


a„  = 


Ze  (ri> 


(ri  - ri-l} (ri  - ti+l) (ri  " ri+2) 


Ze  (ri+l> 


(ri+l  ri-l)(ri+l  ri)(ri+l  ri+2) 


a,  = 


Ze  <ri+2> 


(ri+2  “ ri-l)(ri+2  ri)(ri+2  ri+l* 


The  base  interval,  Ar,  is  always  constant  for  a particular  range  incre- 
ment mode,  so  that  the  programming  steps  are  greatly  simplified.  The 
off-diagonal  elements  of  the  matrix  used  to  calculate  the  value  of 
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Z (r')  are  constant  and  thus  only  the  diagonal  elements  of  the  matrix 
e 

need  be  computed.  If  a linear  interpolation  is  required,  then  only  the 
two  points  immediately  adjacent  to  r'  are  used. 

6.  Once  the  value  of  Z^r'.h)  is  known,  the  value  of  h is  tested 
to  determine  whether  Ze(r',h)  is  to  be  added  into  the  lower  (0-deg  to 
15  kft),  the  middle  (15  to  35  kft) , or  the  upper  (35  to  50  kft)  PVSZ 
layer. 

If  the  range  increment  is  2 km  instead  of  1 km,  then  all  the 
increments  in  this  scheme  are  simply  doubled  (Ex.  r^+^  becomes  r^+2’ 
ri+2  becomes  etc.) ; the  interpolation  scheme  itself  remains  un- 

changed. This  interpolation  along  the  data  radial  is  accomplished  by 
a subroutine  called  INTRP2,  which  also  accomplishes  the  vertical 
integration  of  liquid  water  content  (VIL) . This  procedure  uses  an 
algorithm  developed  by  Greene  (1971)  and  is  based  on  the  relationship 

M*  = 3.44  x 10-6  J t0p  Ze4/7  dh  (26) 

^bottom 

-2 

where  M*  is  the  VIL  in  units  of  kgm  and  h is  the  height  in  meters. 
Conversion  to  Rectangular  Coordinate  System 


The  other  major  interpolation  scheme  is  used  to  map  the  reflectiv- 
ity values  from  cylindrical  coordinates  into  a rectangular  coordinate 
system  (x,y,h)  where  x and  y are  the  rectangular  coordinates  on  a flat 
Earth's  surface  and  h is  held  constant.  The  technique  is  termed  a 
quadratic  interpolation  scheme  but  is  basically  a finite-difference 
form  of  the  Taylor  series  expansion  in  two  dimensions,  truncated  after 
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the  second-order  terms.  Greene  (1971)  presents  a complete  discussion 

of  the  technique.  Briefly,  however,  the  scheme  entails  selection  of 

the  nearest  six  (r,a)^  grid  points  from  the  nine  possible  grid  points 

closest  to  the  (x,y)^  point  in  question.  For  example,  in  Figure  5,  the 

it*1  point  is  the  closest  to  the  desired  point  (x.y)^.  Therefore,  the 

points  (ri>ai),  are  always  selected  along  with 

(r.,a..,)  and  (r,,a.  ,).  The  sixth  point  is  determined  by  the  sign  of 

Sx  and  6y.  The  sixth  point,  in  any  case,  will  always  be  such  that  the 

point  (x,y),  is  encompassed  by  the  nearest  four  points.  With  Sx  = a'  - 
h 

cu  and  Sy  = r'  - r , the  sixth  point  is  chosen  based  on  the  following 
conditions : 


1.  Sx  > 0 


Sy  > 0 


2 .  <Sx  > 0 


Sy  < 0 


3.  Sx  < 0 


6y  > 0 


4.  Sx  < 0 


Sy  < 0 


^ri+l,cti+l* 


(ri-l*°i+l* 


(ri+l*°i-l) 


(ri-l’“i-l) 


The  quadratic  interpolation  function  is  given  by 

2 2 

Z (x,y)  = a,  + a„5x  + a„Sy  + a. Sx  + acSy  + a SxSy  (27) 

e i z J h j b 

where  Sx  is  the  percentage  of  distance  from  the  azimuth  of  the  (x,y) 
grid  point  to  the  nearest  data  azimuth,  and  Sy  is  the  percentage  of 
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distance  from  the  grid  point  to  the  nearest  range  gate. 


The  a^  coefficients  are  defined  as 


al  = Ze  ^i’V 


a2  " [Ze(ri’ai+l)  Ze(ri,ai-1)]  1 2 


a3  ^Ze^ri+l,oti^  “ Ze^ri-l,0ti^  ^ 2 


a4  = Ze(ri’ai+1)  ‘ ai  - a2 


a5  ^e^ri_i,cti^  a!  + a3 


and  a,  is  determined  according  to  the  sign  of  & and  to  be 
o 


1. 


6x  _>  0 
6y  _>  0 


a6  Ze^ri+l,ai+l^  al  a2 


a3  “ a4  a5 


2. 


6x  ^ 0 
6y  < 0 


a6  = al  + a2  “ a3  + a4  + a5 


- Ze(ri-l’ai+l) 


3. 


6x  < 0 

<Sy  21  0 

6x  < 0 
6y  < 0 


a6  = al  ~ a2  + a3  + a4  + a5 


Ze(ri+l’“i-l^ 


a6  = Ze(ri-l>ai-l)  " al  + a2 


+ a3  - a4  - a5 


(28) 


If,  however,  the  quadratic  interpolation  results  in  a negative 


value  for  Ze(x,y),  then  a linear  interpolation  scheme  which  uses  the 
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four  data  points  immediately  surrounding  (x,y)  is  employed  to  obtain 
the  interpolated  value  of  Z£(x,y).  This  linear  interpolation  takes  the 
form 

A " Ze(ri+l’oti-l)  + ^Ze(ri+l’°i)  " Ze(ri+l’ai-l) ] (6x) 

B = Ze(ri,a1  .j^)  + fze(ri»ai)  ~ Ze^ri’ai-1^ 

Ze(x,y)  = B + (A  - B)(6y)  (29) 

If  the  grid  point  falls  on  a data  radial,  then  a linear  interpo- 
lation is  performed  along  that  radial  between  the  range  gates  on  either 
side  of  the  grid  point. 

If  either  2 deg  or  2 km  data  are  being  processed,  the  appropriate 

increments  or  ri+^,  etc-)  are  all  doubled,  as  they  were  in  the 

Lagrangian  cubic  interpolation  scheme.  The  interpolations  involved  in 

this  conversion  from  cylindrical  to  rectangular  coordinates  are  per- 

t 

formed  by  a subroutine  called  QD2. 

• • Mechanics  of  the  Computer  Program 

This  computer  program  takes  digital  radar  reflectivity  data  from 
the  TAMU  digital  weather  radar  system,  which  has  been  recorded  on  9- 
track  tape,  and  uses  these  data  to  compute  and  produce: 

1.  Zero-degree  tilt  reflectivity  maps; 

2.  Partial  vertically-summed  reflectivity  (Z)  maps  (PVSZ)  for 
three  layers  (sfc  to  15  kft,  15  to  35  kft,  and  35  to  50  kft);  and 

3.  Vertically  integrated  liquid  water  content  (VIL)  maps. 

The  program  will  process  data  recorded  in  all  four  modes  possible  with 
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this  radar  system:  1 deg  by  1 km,  1 deg  by  2 lan,  2 deg  by  1 km,  and  2 deg 
by  2 km. ^ However,  it  will  only  process  data  from  one  radar  at  a time, 
either  the  3 cm  or  the  10  cm  radar.  The  complete  program,  including 
comments  and  job  control  language  (JCL) , is  provided  in  Appendix  A for 
reference. 

Program  Inputs 

The  program  requires  three  inputs,  two  on  punch  cards  and  one  on 
magnetic  tape.  The  first  data  card  is  the  radar  calibration  data.  It 
is  used  by  the  subroutine  LOAD  to  reproduce  the  radar  calibration  curve, 
which  converts  each  of  the  256  digital  integers  of  the  raw  reflectivity 
data  (numbered  from  0 to  255)  to  equivalent  dBm  values.  The  data  card 
contains  eight  integers  and  one  real  number,  all  positive,  as  in  the 
example  below: 


32  27  60  98  139  171  201  233  107.5 

The  real  number  is  the  radar's  minimum  detectable  signal  (MDS)  in  dBm. 
The  first  integer  is  the  number  of  digital  integers  that  are  assigned 
the  MDS  for  their  equivalent  dBm  value  (32) . The  next  integer  is  the 
total  number  of  digital  integers  assigned  values  before  the.  equivalent 
dBm  value  increases  to  10  dBm  above  MDS  (from  -107.5  to  -97.5  dBm) (37). 
The  third  integer  is  the  total  number  of  digital  integers  assigned 
values  before  the  equivalent  dBm  value  increases  another  10  dBm  (total 


When  first  encountered  by  this  investigator,  this  computer  pro- 
gram would  only  process  data  recorded  in  the  1 deg  by  1 km  mode.  This 
investigator  made  all  subsequent  modifications  to  the  program  to  enable 
it  to  process  data  recorded  in  the  other  three  modes  of  collection. 


from  -107.5  to  -87.5) (60).  This  process  continues  until  all  eight 
integer  locations  on  the  data  card  have  been  filled. 

The  second  data  card  provides  the  initializing  information  for  the 
radar  scan  sequence.  One  of  these  cards  is  required  for  each  complete 
scan  sequence  that  is  to  be  processed.  An  example  is 

162-3826935106402002321008200515 

which  is  broken  up  into  its  individual  components  and  defined  below 
(see  Figure  6) . 

(162) (-38) (269) (351) (06) (40) (200232) (10) (08) (200515) 

ILEFT  is  the  distance  east  (-)  or  west  (+)  from  TAMU  to  the 
•lower  left  corner  of  the  100  x 100  km  grid,  plus  two.  For  example,  160 
km  west  = 162. 

JDOWN  is  the  distance  north  (-)  or  south  (+)  from  TAMU  to  the 
lower  left  corner  of  the  100  x 100  km  grid,  plus  two.  For  example, 

40  km  north  = -38. 

KAZ  is  the  desired  starting  azimuth  for  data  processing  to  load 
the  grid  box,  moving  clockwise.  KAZ  may  be  several  degrees  short  of 
the  grid  box,  but  may  not  err  on  the  long  side.  Example  is  269  deg. 

LASTAZ  is  the  last  azimuth  needed  to  load  the  grid  box,  and  it 
may  be  several  degrees  beyond  the  end  of  the  box,  but  not  short  of  it. 
Example  is  351  deg. 

NTLT  is  the  last  elevation  angle  incremented  by  1 deg  in  the  tilt 
sequence.  Example  is  6 deg. 
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KRANGE  is  the  radial  distance  to  the  nearest  point  in  the  grid 
box.  It  also  may  err  on  the  short  side,  but  not  on  the  long  side. 
Example  is  40  km. 

ISTRT  is  the  starting  time  of  the  scan  sequence,  local  time,  in 
hours,  minutes,  and  seconds.  Example  is  20:02:32. 

IRIT  is  the  tape  unit  onto  which  the  output  data  are  written, 
defined  in  the  JCL.  Example  is  10. 

IRR  is  the  tape  unit  from  which  the  input  data  are  read,  defined 
in  the  JCL.  Example  is  08. 

NENDT  is  the  end  time  of  the  scan  sequence,  as  above.  Example  is 
20:05:15. 
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MAGNETIC  TAPE  DATA  FORMAT 


TRACK  NUKE II  (LOGICAL  TRACKS) 


RADAR  BEING  RECORDED  (3  or 
OEGREES  (AZ) 

RANGE  BIN  SELECTION 
NUFBEK  OF  PULSES  INTEGRATED 
RATE  ERROR 
PARITY  ERROR 


AZIMUTH  INFORMATION 


ELEVATION  INFORMATION 


**  ID  code  shall  be 
hard  wired  to  Indi- 
cate either  3 or  10 
cm  radar  recorder. 
Playback  shall  be  on 
the  3 cm  radar  tape 
recorder. 


RECORD) 


1 

2 

3 

4 

5 6 

7 8 9 

1 

1 

1 

1 

DATE 

100's 

1 

1 

1 

1 

DATE 

10‘S 

1 

1 

I 

1 

DAIE 

1 'i 

1 

1 

I 

I 

HOUR 

10's 

1 

I 

1 

1 

HOUR 

1 's 

I 

1 

1 

1 

MIN 

10‘s 

1 

l 

1 

I 

KIM 

I's 

1 

1 

I 

I 

SEC 

10'S 

I 

! 

1 

1 

SEC 

r* 

I 1 1 I SPARE 

I 

1 

I 

1 

SPARE 

I 1 1 I SPARE 

1 1 I 1 ID 

1 

I 

) 

J 

DS 

1 1 1 1 RI 

1 1 1 I SA 

1 1 1 1 RE 

1 

1 

I 

1 

PE 

1 

l 

1 

1 

AZ 

100's 

1 

1 

1 

1 

AZ 

10'j 

1 

1 

1 

1 

AZ 

1*s 

1 

I 

1 

1 

EL 

10's 

1 

1 

I 

I 

EL 

Ps 

I 

1 

1 

1 

EL 

O.Ps 

8-BIT 

DVIP  OATA 

RECORDED  AT  EITHER  1*  OR 
2°  OF  AZIMUTH.  IN  1 Km 
RANGE  MODE:  Z OATA  SETS 
PER  RECORD.  IN  2 KM 
RANGE  MODE:  A DATA  SETS 
PER  RECORD 


CODE 

l 

2 

ID 

3 

10  era 

OS 

P 

2° 

RI 

1 Km 

2 Km 

SA 

15 

31 

Fig.  7.  Magnetic  tape  data  recording  format 
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Each  logical  record  contains  either  two  or  four  radials  of  data, 
depending  on  the  range  increment  selected,  and  from  two  to  eight 
degrees  of  azimuthal  area  coverage,  depending  on  the  azimuth  increment 

selected.  In  the  1 km  mode,  the  first  radial  of  data  is  contained  in 

/ 

IHWORD  (1-430)  and  the  second  in  IHWORD  (431-860).  In  the  2 1cm  mode, 
the  first  radial  is  contained  in  IHWORD  (1-215) , the  second  in  IHWORD 
(216-430),  the  third  in  IHWORD  (431-645),  and  the  fourth  in  IHWORD 
(646-860).  Likewise,  in  the  1 deg  mode,  the  first  radial  of  data 
represents  a,  the  second  a + 1,  the  third  and  fourth  (in  the  2 km  mode) 
a + 2 and  a + 3 deg.  In  the  2 deg  mode,  the  radials  represent  a,  a + 2, 
a + 4,  and  a + 6 deg,  each  increment  being  two  degrees  wide,  for  a 
total  areal  coverage  of  eight  degrees. 

Main  Program 

The  first  thing  the  main  program  does  is  assign  a value  to  WAVEC, 
which  represents  the  radar  constant,  and  is  different  for  each  radar 
(see  pg.  14).  This  value  must  match  that  of  the  radar  that  generated 
the  data  being  processed. 

The  next  thing  accomplished  is  the  reconstruction  of  the  radar 
calibration  curve  and  the  assignment  of  equivalent  dBm  values  to  the 
256  digital  integers  by  the  subroutine  LOAD. 

Then  the  initializing  data  card  is  read,  and  most  of  the  values 
contained  on  it  printed  out  below  the  radar  calibration  curve  for 
reference. 

The  program  then  calls  the  subroutine  DIVIDE,  which  computes  the 
number  of  interpolation  points  contained  in  each  PVSZ  layer  above  each 


r. TTy  • 


range  gate,  and  stores  them  in  the  arrays  IDIVl(I),  IDIV2(I),  and 
IDIV3(I) . 
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Then  the  first  logical  record  of  data  is  read  from  tape  and  the 
time  of  the  record  ITIME  is  tested  against  the  desired  start  time 
1STRT.  If  ITIME  is  less  than  ISTRT,  another  record  is  read  and  tested, 
until  one  is  found  that  is  within  the  desired  time.  This  data  record 
is  then  tested  for  range  and  azimuth  increments,  and  the  appropriate 
increments  (IADD  for  azimuth  and  IADDR  for  range)  are  assigned  for  use 
throughout  the  remainder  of  the  program.  The  next  test  is  to  see  if 
the  azimuth  of  the  first  data  radial  falls  within  the  area  of 
interest,  that  is,  between  KAZ  and  LASTAZ.  If  so,  processing  begins, 
but  if  not,  the  azimuth  of  the  second  radial  is  tested,  etc.  If  the 
last  radial  in  that  record  does  not  fall  within  the  required  bounds, 
then  another  record  is  read,  and  the  procedure  continues  until  the 
first  data  radial  is  found  that  falls  within  the  area  of  interest. 

Once  a radial  is  found  that  lies  between  KAZ  and  LASTAZ,  it  is 
converted  from  azimuth  to  coded  form  II,  where  KAZ  is  II  = 1,  and  II 
increases  sequentially  clockwise.  Then  the  elevation  angle  is  tested, 
first  by  a library  subroutine  ATOI,  and  then  against  a reference  angle, 
PTILT.  ATOI  checks  to  ensure  that  the  integer  to  binary  conversion  of 
the  elevation  angle  was  done  correctly.  A non-zero  result  returned 
from  this  test  ends  the  program.  PTILT  is  a reference  angle  that  is 
initialized  at  0.0  deg  at  the  start  of  the  program.  If  the  elevation 
angle,  TILT,  is  larger  than  PTILT,  it  means  that  the  zero-tilt  portion 
of  the  scan  is  over,  and  the  program  prints  out  the  zero-tilt  reflec- 
tivity map,  increments  PTILT  up  to  the  current  value  of  TILT,  then 
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begins  processing  the  data  radials  that  go  into  the  three  PVSZ  layers. 

If,  however,  TILT  becomes  more  than  0.5  deg  smaller  than  PTILT,  then 
the  program  assumes  that  th<  scan  sequence  is  over,  and  it  terminates 
processing  data  radials,  prints  out  all  the  remaining  maps,  and  pre- 
pares to  start  a new  scan  sequence. 

Now  the  actual  processing  of  the  reflectivity  data  begins.  The 
first  180  elements  of  the  first  data  radial,  IHWORD(1-430)  or  (1-215), 
are  loaded  into  the  array  ICLASS(T.)  in  positions  21-200,  since  the 
first  20  km  of  range  are  blanked  out  by  the  DVIP,  and  the  first  ele- 
ment recorded  is  the  21st  range  gate.  Then  the  reflectivities  for 
100  km  of  range  are  computed,  beginning  at  KRANGE  plus  one  range  incre- 
ment. If  the  digital  value  in  ICLASS  is  less  than  or  equal  to  34,  the 
reflectivity  ZW(II,I),  is  set  to  0.0,  since  values  of  34  or  less 
represent  reflectivities  below  the  MDS  of  the  radar.  The  reflectivity 
is  computed  according  to  Eq  (10)  on  pg.  13.  If  the  tilt  angle  is  less 
than  0.5  deg,  then  the  interpolated  reflectivity  data  are  added  into 
the  lower  PVSZ  layer  after  being  divided  by  the  appropriate  IDIV1 
divisor.  The  interpolated  reflectivity  factors  at  each  range  gate  in 
a layer  cannot  simply  be  summed  vertically,  because  there  are  a dif- 
ferent number  of  factors  to  be  summed  in  each  layer.  The  results  would 
be  erroneous  gradients  of  reflectivity.  Thus  the  value  of  each  re- 
flectivity factor  is  divided  by  the  number  of  tilt  angles  of  data  that 
are  actually  available  in  the  applicable  layer  for  that  range.  If  the 
tilt  angle  is  greater  than  0.5  deg,  then  the  program  waits  until  all 
100  reflectivities  have  been  computed,  and  then  sends  them  to  the 
INTRP2  subroutine.  INTRP2  interpolates  the  reflectivities  along  the 
radial,  converts  the  data  to  cylindrical  coordinates,  and  assigns  the 
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reflectivities  to  the  appropriate  PVSZ  layer,  where  they  are  summed 
and  divided  by  the  appropriate  IDIV  divisor.  INTRP2  also  accomplishes 

the  vertical  integration  of  liquid  water  content  (VIL) . This  com- 
pletes the  processing  of  the  first  data  radial.  The  remaining 
radials  in  that  record  are  processed,  and  another  record  read,  etc., 
until  all  the  data  for  that  scan  sequence  have  been  read.  When  either 
a time  is  read  that  exceeds  the  end  time  of  the  sequence,  or  the  tilt 
angle  becomes  more  than  0.5  deg  less  than  PTILT,  the  processing 
terminates  and  the  output  of  the  four  remaining  maps  begins. 

The  three  PVSZ  arrays  represent  the  three  layers  of  the  storm 
complex,  and  the  value  of  the  summed  reflectivity  factors  at  any  given 
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ments,  is  provided  in  Appendix  B.  The  values  that  are  shown  for  the 
various  parameters  are  the  ones  that  were  used  to  produce  all  of  the 
maps  contained  in  this  investigation. 

The  RDRMAP  program  uses  the  output  arrays  of  the  RADAR  program 
(Sieland's  program)  as  recorded  on  magnetic  tape,  to  produce  contoured 
reflectivity  maps.  The  program  will  plot  the  entire  51  x 51  reflec- 
tivity array,  or  any  portion  of  it  that  is  of  interest  to  the  program- 
mer. The  only  input  required  by  the  program  is  the  output  tape  from 
the  RADAR  program.  All  the  plotting  parameters  are  set  or  initialized 
inside  the  program.  The  JCL  required  for  the  program  is  contained  and 
explained  in  Appendix  B.  The  main  points  of  the  program  are  dis- 
cussed briefly  below. 

All  of  the  subroutines  in  this  program  make  extensive  use  of 
coordinate  locations.  It  is  important  to  realize  that  the  locations 
addressed  by  this  program  are  based  on  a rectangular  coordinate  system, 
with  the  origin  (0,0)  in  the  lower  left  corner,  and  not  the  con- 
ventional Fortran  two-dimensional  variable  addressing  system.  However, 
the  array  address  of  the  lower  left  corner  of  the  plot  is  (1,1). 
Coordinates  may  have  two  forms,  either  integer  or  real  numbers.  The 
integers  have  a range  of  1 to  1024,  where  each  unit  represents  0.01 
inch  on  the  output  map.  The  real  numbers,  on  the  other  hand,  refer  to 
a specific  coordinate  position  (x  = 6.),  and  the  correct  distance  down 
the  axis  in  inches  is  determined  internally  by  the  program.  The  maxi- 
mum size  of  the  output  plot,  for  these  purposes,  is  10.24  in  x 10.24  in, 
including  labels,  and  the  programmer  may  utilize  any  portion  of  that 
field. 
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The  first  subroutine  called  is  SET,  which  physically  sets  up  the 
graph.  The  form  is  CALL  SET  (XA,  XB,  YA,  YB,  XC,  XD,  YC,  YD,  LTYPE) , 
where  the  arguments  are  defined  below. 

The  first  four  arguments  (XA,  XB,  YA,  YB)  establish  the  size  of 
the  area  to  be  plotted.  XA  and  YA  are  the  position  of  the  origin  of 
the  field.  For  example,  XA  = 110,  YA  = 110  indicates  that  the  origin 
of  the  plot  is  to  be  1.10  in  to  the  right  and  1.10  in  up  from  the  lower 
left  corner  of  the  field.  XB  and  YB  indicate  the  length  of  the  x and 
y axes,  starting  from  the  point  (XA,YA) . The  second  four  arguments 
(XC,XD,YC,YD)  define  the  range  of  x and  y along  their  respective  axes. 
XC  determines  the  coordinate  assigned  to  the  lower  left  corner  of  the 
grid,  and  XD  the  lower  right  corner.  For  example,  XC  = 0.,  XD  = 25. 
means  the  scale  along  the  x axis  will  range  from  0 to  25,  for  a total 
of  26  coordinate  addresses.  The  same  is  true  for  the  y axis  and  YC 
(bottom)  and  YD  (top).  The  last  argument,  LTYPE,  determines  the  type 
of  mapping  to  be  done,  either  log  or  linear,  in  the  x and  y directions. 
LTYPE  = 1 is  linear  in  both  directions. 

The  next  subroutine  called  is  PERIM,  which  has  the  form  CALL 
PERIM  (MCRX,  MINRX,  MGRY,  MINRY).  This  subroutine  draws  a box  around 
the  area  specified  by  SET,  and  places  tick  marks  around  the  outside 
perimeter  of  the  box  as  determined  by  the  arguments  of  PERIM.  MGRX 
and  MGRY  determine  the  number  of  major  divisions  along  the  x or  y 
axis.  It  is  important  to  note  that  the  divisions  are  the  holes,  and 
not  the  tick  marks.  Likewise,  MINRX  and  MINRY  determine  the  number  of 
minor  divisions  per  major  division.  Once  again,  these  are  the  holes, 
and  not  the  tick  marks. 
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The  next  subroutine  called  is  PWRX,  which  puts  the  labels  on  the 
map.  This  subroutine  must  be  called  once  for  each  label.  The  form 
is  CALL  PWRX  (X,Y, ID, N.ISIZE, ITHETA, ICNT) . The  first  two  arguments 
(X,Y)  are  the  positioning  coordinates  for  the  characters  to  be  drawn, 
and  may  be  either  integer  or  real.  For  example,  X = 15. , Y = 90 
indicates  that  the  character  is  to  be  drawn  .90  in  up  from  the  bottom 
of  the  field  (y  coordinate)  and  at  a position  corresponding  to  the 
coordinate  X = 15  along  the  x axis  of  the  graph. 

The  argument  ID  is  the  character  string  to  be  drawn  and  takes  the 
form  of  this  example:  SH'IRU'MAP,  where  8H  reserves  eight  spaces  for 
hollerith  code  (’IRU'MAP),  IRU  specifies  the  type  of  printing,  and  MAP 
are  the  characters  printed.  The  coding  of  IRU  is: 

I = .13  in  R = Roman  U = upper  case 

P = .21  in  G = Greek  L = lower  case 

The  argument  N is  simply  the  number  of  hollerith  characters  in  ID, 

ISIZE  determines  the  size  of  the  characters,  and  may  be  either  integer 
or  real.  If  it  is  real  (floating  point),  then  it  is  used  as  a multi- 
plication factor  of  I or  U in  ID.  If  it  is  integer,  it  has  a different 
function  that  is  not  relevant  to  this  application  of  this  program. 

The  argument  ITHETA  determines  along  which  orientation  the  char- 
acters are  printed,  if  ITHETA  = 0,  the  character  strings  are  printed 
in  the  x direction,  oriented  vertically  (normal  orientation).  If 
ITHETA  = 1,  the  character  strings  are  printed  in  the  y direction, 
oriented  horizontally. 

ICNT  is  the  centering  option  and  is  defined  as  follows: 

ICNT  = 0 (X,Y)  is  the  center  of  the  first  character. 
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ICNT  = l (X,Y.)  is  the  center  of  the  entire  string. 

ICNT  = 2 (X,Y)  is  the  left  edge  of  the  first  character. 

ICNT  = 3 (X,Y)  is  the  center  of  the  last  character. 

ICNT  = 4 (X,Y)  is  the  right  edge  of  the  last  character. 

The  next  subroutine  called  by  the  program  is  the  CONREC  subroutine 
itself,  and  this  is  the  one  that  actually  does  the  plotting  and  con- 
touring of  the  data  array.  The  form  of  this  subroutine  is  CALL  CONREC 
(Z  (X, Y)  ,L ,M, N , FLO  ,HI , FT.NC , NSET , NIII  ,NDOT)  . 

The  first  argument,  Z(X,Y),  is  the  array  coordinate  of  the  lower 
left  corner  of  the  field  to  be  contoured.  For  example,  Z(l,l)  will 
fill  the  entire  grid  box,  but  Z(2,2)  will  leave  a blank  margin  one  grid 
unit  wide  around  the  inside  of  the  box.  The  argument  L is  the  first 
dimension  of  the  array  Z in  the  calling  program.  The  next  two  argu- 
ments (M,N)  are  the.  number  of  data  values  to  be  plotted  in  the  x and  y 
directions.  FLO  and  HI  arc  the  values  of  the  lowest  and  highest  contour 
levels  to  be  drawn,  respectively,  and  FINC  is  the  increment  between 
contour  levels. 

The  last  three  arguments  (NSET,  NHI , NDOT ) control  scaling,  extra 
information  on  the  plot,  and  the  line  pattern,  and  should  not  be 
changed  for  use  with  this  program. 

The  last  subroutine  called  by  this  program  is  PLOT,  and  it  pro- 
vides pen  control,  and  is  also  used  to  signal  enu-of-plot.  The  form 
is  CALL  PLOT(X,Y, IPEN)  where  (X,Y)  are  the  page  coordinates  in  inches 
of  the  point  to  which  the  pen  is  to  be  moved.  For  example,  (12.,  0.) 
moves  the  pen  12  inches  in  the  x direction.  The  magnitude  of  IPEN 
specifies  the  up  or  drT,n  condition  of  the  pen  during  its  movement  to 
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(X,Y) . IPEN  = 2 is  for  pan  down,  and  IPEN  = 3 is  for  pen  up.  If  the 
sign  of  IPEN  is  positive,  the  movement  results  in  a normal  plot 
entrance.  However,  if  the  sign  of  IPEN  is  negative,  the  movement  is  an 
end-of-plot  entrance,  and  after  the  pen  moves  to  (X,Y) , the  value 
x = 0,  y = 0 is  assigned  to  the  new  pen  location,  and  the  plot  number 
is  incremented  by  one.  Additionally,  IPEN  = 99  is  used  to  empty  the 
output  buffer  at  the  completion  of  the  plot,  and  this  must  be  done 
prior  to  terminating  the  program. 

The  remainder  of  the  RDRMAP  program  is  explained  by  means  of 
comments  in  the  progrm  as  it  appears  in  Appendix  B. 

Radar  Volume  in  Space 

It  is  easy  to  see  from  Figure  8 that  the  volume  in  space  of  the 
various  range  bins  increases  rapidly  both  as  the  range  increases,  and 
as  the  collection  mode  is  changed  from  the  1 deg  x 1 km  mode  to  one  of 
the  other  modes.  The  illuminated  radar  volume,  which  is  one  half  the 
geometric  volume,  is  given  by 


V = 


h 

2 


(31) 


where  6 is  the  beamwidth  in  radians  (assumed  to  be  circular  and  1 deg 

in  diameter),  r is  the  range  in  km,  h is  the  length  of  the  range  bin  in 

3 

km,  and  V is  the  volume  in  km  . Table  3 shows  the  approximate  illumi- 
nated radar  volumes  enclosed  by  range  bins  in  each  of  the  four  collec- 
tion modes,  at  ranges  of  50  and  150  km.  Note  that  the  volume  enclosed 
by  the  2 deg  x 2 km  mode  is  four  times  that  enclosed  by  the  1 deg  x 1 km 
mode,  and  that  the  1 deg  x 1 km  volume  at  150  km  is  nine  times  the 
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TABLE  3.  Radar  Volumes  at  Near  and  Far  Ranges. 


3 

Volume  in  km 

Collection  Mode 

50  km  range 

150  km  range 

1 deg  x 1 km 

0.30 

2.70 

2 deg  x 1 km 

0.60 

5.40 

1 deg  x 2 km 

0.60 

5.40 

2 deg  x 2 km 

1.20 

10.80 

1 deg  x 1 kin  volume  at  50  km.  Even  more  dramatic  is  that  the  2 deg  x 

2 1cm  volume  at  150  km  is  36  times  as  large  as  the  1 deg  x 1 km  volume 

at  50  km,  remembering  that  the  physical  volumes  in  space  are  twice  the 

radar  volumes  shown  in  the  table.  Based  on  these  considerations,  it 

is  only  reasonable  to  suspect  that  a single  digital  value  representing 

3 

the  mean  reflectivity  in  a volume  of  10.8  km  may  not  represent  that 
volume  as  adequately  as  four  digital  values,  each  being  the  mean  re- 
flectivity of  one  fourth  of  that  larger  volume.  The  results  of  the 
investigation  into  this  and  other  questions  are  presented  in  the  next 
chapter. 


-•-C-Vr-vy  ■ 
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CHAPTER  IV 

PRESENTATION  AND  DISCUSSION  OF  RESULTS 

Data  Acquisition  and  Reduction 

The  radar  data  used  in  this  investigation  were  produced  and 
recorded  by  the  TAMU  Digital  Weather  Radar  System  during  the  evening  of 
9 May  1977,  as  an  area  of  light  and  moderate  thunderstorm  activity 
moved  through  the  College  Station  area.  Data  were  recorded  simul- 
taneously in  both  wavelengths  (3  and  10  cm)  and  as  near  simultaneously 
as  possible  in  all  four  modes  of  collection.  Only  zero-degree  tilt 
data  were  recorded,  and  all  the  data  were  recorded  within  a period  of 
only  3.5  minutes.  The  digital  data  were  processed  by  Sieland's  RADAR 
computer  program,  and  then  the  final  output  arrays  of  zero-tilt  re- 
flectivity data  were  plotted  on  the  CALCOMP  plotter.^  The  contoured 
reflectivity  maps  presented  here  are  reductions  of  the  originals,  which 
were  approximately  ten  inches  square. 

Data  Analysis 

The  data  generated  by  this  investigation  were  analyzed  both 
qualitatively  and  quantitatively.  The  qualitative  analysis  consisted 
of  physically  superimposing  the  contoured  reflectivity  maps  on  one 
another  and  examining  them  over  a light  source. 


All  contoured  maps  were  produced  using  a library  subroutine 
named  CONREC  and  plotted  on  the  CALCOMP  plotter  at  the  TAMU  Data 
Processing  Center.  All  isopleths  on  these  maps  are  in  units  of  dBZ. 
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The  quantitative  analysis  consisted  o£  a statistical  comparison  of 
corresponding  points  in  two  different  data  sets  (reflectivity  maps)  to 
determine  whether  there  was  any  significant  difference  between  them. 

The  test  employed  was  the  Student's  t test,  based  on  the  assumptions 
discussed  below. 

The  first  assumption  is  that  the  data  value  assigned  to  each  grid 
point  in  the  51  x 51  output  arrays  represents  the  mean  reflectivity 
(sample  mean,  X^)  for  that  grid  square,  and  that  each  data  point  is 

derived  from  31  samples.  We  also  assume  that  the  sample  variance 

2 2 
(s^  ) is  the  same  for  each  data  set  and  is  equal  to  (5.57)  dB.  The 

null  hypothesis  is  Hq:  = p2>  or  that  there  is  no  significant  dif- 

ference between  the  data  values  (means)  of  two  corresponding  grid 
points  on  different  reflectivity  maps  (1  deg  x 1 1cm  map  versus  2 deg 
x 2 km  map,  etc.).  The  test  statistic  is  then  given  by 

(X.  - X ) 

C " .~2  2TT/2  (32) 

(K  • %) 

2 2 2. 

where  s is  the  pooled  variance  since  s^  = s2>  and  the  critical  region 

is  defined  by 


t > t 


(1  - a/2)  (n^  + n2  - 2) 


The  1 deg  x 1 km  map  served  as  the  standard  and  the  other  three  modes 
were  tested  against  it.  The  test  was  applied  stepwise  to  every  indi- 
vidual pair  of  corresponding  grid  points  in  the  two  maps  being  compared, 
and  the  null  hypothesis  was  either  rejected  or  not  rejected  for  each 
pair  of  points.  The  pairs  of  points  were  tested  at  three  levels  of 
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significance,  ct  = .01,  .05,  and  .10,  where  a significance  level  of 
a = .01  means  that  there  is  only  a 1%  chance  of  rejecting  a true 
hypothesis;  tliat  is,  there  is  only  a 1%  chance  of  saying  that  there  is 
a significant  difference  between  the  two  points,  when  in  truth  there  is 
not.  Finally,  the  percentages  of  pairs  of  points  that  were  signifi- 
cantly different  from  one  another  were  tabulated. 

Comparison  of  3 cm  Data  in  the  Near  Range 

The  first  analyses  made  were  of  the  3 cm  data  collected  in  all 
four  modes  within  the  near  range,  or  from  20  to  60  km  from  the  radar. 
The  zero-tilt  reflectivity  maps  for  the  four  modes  are  shown  in 
Figures  9 through  12.  Additionally,  the  results  of  the  statistical 
analyses  of  the  different  modes  are  also  displayed  in  map  form. 

The  2 deg  x 1 km  Mode 

The  first  comparison  was  between  the  1 deg  x 1 km  mode  (Figure  9), 
which  was  taken  to  be  the  standard  mode  in  all  the  analyses,  and  the 
2 deg  x 1 km  mode  (Figure  10). ^ The  2 deg  x 1 km  map  retains  most  of 
the  features  of  the  1 deg  x 1 km  map,  although  there  are  a fair  number 
of  differences  apparent.  In  the  2 deg  x 1 km  map,  there  is  some 
smoothing  of  the  details  contained  in  the  1 deg  x 1 km  map,  and  a 
tendency  to  combine  small,  adjacent  echo  features  into  single  features, 
most  likely  caused  by  the  wider  integration  interval  in  the  direction 

^The  apparent  strong  reflectivity  gradients  in  the  upper  left 
hand  corner  of  all  the  near  range  maps  should  be  disregarded.  This 
apparent  gradient  is  due  to  the  range  blanking  of  the  UVIP,  which 
blanks  out  all  data  within  20  km  of  the  radar. 
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KILOMETERS  SOUTH  OF  CLL 


KILOMETERS  EAST  OF  CLL 


Fig.  10.  Three  cm  reflectivity  map  in  the  near  range,  2 deg  x 1 km 
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of  antenna  rotation.  This  phenomena  also  results  in  the  expansion  or 
shift  of  some  reflectivity  features  in  the  direction  of  antenna  rota- 
tion (clockwise).  Additionally,  several  of  the  reflectivity  features 
with  closed  isopleths  exhibit  slightly  higher  dBZ  values  on  the 
2 deg  x 1 km  map,  although  not  more  than  5 dBZ  higher,  and  the  high- 
est value  on  both  maps  is  40  dBZ.  It  also  appears  tiiat  the  deeper  a 
reflectivity  feature  is  embedded  inside  the  precipitation  area,  the 
greater  the  amount  of  smoothing  and  reduction  in  size  of  features  that 
occur  on  the  2 deg  x 1 km  map,  as  shown  in  the  southeast  corner  of 
these  maps.  The  two  maps  do  agree  quite  well  on  the  two  areas  of 
significantly  reduced  reflectivity  (near  zero)  shown  at  A and  B.  It 
is  possible  that  these  two  areas  represent  BWERs,  but  a positive  de- 
termination can  not  be  made  without  the  maps  for  the  next  several  CAZM 
levels  above  these  areas,  or  a complete  lower  level  PVSZ  map, 
integrated  from  the  0-deg  tilt  level  to  15  left. 

The  statistical  comparisons  of  these  two  maps  at  three  levels  of 
significance  (a  = .01,  .05,  and  .10)  are  shown  in  Figures  13  through 

g 

15.  Figure  13  shows  that  at  a significance  level  of  .01,  25.4  per 
cent  of  all  the  non-zero  pairs  of  corresponding  data  points  are  sig- 
nificantly different  from  each  other.  This  percentage  increases  to 
35.4  percent  for  a = .05  (Figure  14)  and  to  43.2  percent  for  a = .10 
(Figure  15) . Examination  of  these  figures  shows  that  the  differences 


On  the  statistical  maps,  a blank  means  that  the  two  corresponding 
grid  point  values  on  the  two  maps  being  compared  were  both  zero;  a 
means  that  there  was  no  significant  difference  between  the  two  corre- 
sponding grid  points  at  that  level  of  significance;  an  'xx'  means  that 
the  two  corresponding  grid  points  were  significantly  different  from  one 
another . 
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Fig.  15.  Statistical  comparison  of  the  1 deg  x 1 km  and 
2 deg  x 1 km  modes,  3 cm,  near  range,  at  .10  level  of  significance 
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occur  primarily  in  areas  with  a pronounced  reflectivity  gradient  and  in 
the  areas  of  maximum  reflectivity  values.  This  corresponds  well  with 
the  smoothing  of  detail  and  shifting  of  reflectivity  features  observed 
on  the  contoured  maps.  The  statistical  analysis  maps  also  show  the 
two  areas  of  zero  reflectivity,  A and  B,  quite  well. 

The  1 deg  x 2 km  Mode 

The  second  comparison  was  made  between  the  1 deg  x 1 km  mode  and 
the  1 deg  x 2 km  mode  (Figure  11) . The  correlation  between  these  two 
maps  is  poor,  with  most  of  the  significant  features  of  the  1 deg  x 1 km 
map  either  lost  or  greatly  distorted  on  the  1 deg  x 2 km  map.  Addi- 
tionally, many  of  the  features  appear  to  be  shifted  outwards  in  range 
several  km,  as  might  be  expected  with  the  longer  range  integration 
increment.  The  1 deg  x 2 km  map  also  shows  a considerable  smoothing  of 
the  fine  scale  details  of  the  reflectivity  contours.  The  dBZ  values 
on  the  two  maps  agree  fairly  well , however,  and  the  highest  observed 
dBZ  value  is  once  again  40  dBZ.  Perhaps  the  best  areas  of  agreement  on 
the  two  maps  are  the  two  areas  of  zero  reflectivity  at  A and  B,  al- 
though they  are  somewhat  reduced  in  size  on  the  1 deg  x 2 km  map. 

The  statistical  analysis  at  a = .05  (Figure  16)  shows  that  66.2 
percent  of  all  the  non-zero  pairs  of  data  points  differ  from  one 
another  significantly.  The  points  that  differ  significantly  are  gen- 
erally widespread  everywhere  except  in  the  areas  of  low  reflectivity. 
The  areas  A and  B are  just  barely  discernible  in  this  map. 


• XX 


XX  XX  XX  XX  XX  XX  XX  XX  XX 


XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

• 

• 

XX 

XX 

XX 

XX 

• 

XX 

XX 

• 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

xx 

. 

• 

XX 

XX 

XX 

. 

• 

XX 

XX 

. 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

XX 

. 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

. 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

XX 

. 

XX 

. 

XX 

• 

XX 

. 

. 

XX 

. 

. 

XX 

XX 

XX 

XX 

XX 

. 

. 

. 

. 

. 

. 

. 

XX 

XX 

XX 

• 

XX 

XX 

XX 

XX 

XX 

XX 

• 

• 

XX 

XX 

. 

. 

• 

XX 

• 

• 

• 

XX 

XX 

XX 

XX 

• 

• 

XX 

XX 

• 

• 

XX 

• 

XX 

XX 

• 

XX 

• 

XX 

• 

XX 

XX 

XX 

• 

XX 

• 

XX 

XX 

. 

XX 

. . 

. 

XX 

• 

XX 

XX 

• 

• 

• 

• 

XX 

• 

. 

• 

XX 

XX 

. 

XX 

XX 

• 

. 

. 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

. 

XX 

. 

• 

XX 

xx 

XX 

• 

XX 

XX 

XX 

. 

XX 

XX 

XX 

XX 

. 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

• 

• 

XX 

XX 

XX 

XX 

XX 

. 

. 

XX 

XX 

XX 

XX 

• 

. 

XX 

XX 

XX 

XX 

XX 

XX 

. 

. 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

• 

XX 

• 

XX 

XX 

XX 

x.\ 

XX 

XX 

. 

XX 

XX 

. 

XX 

XX 

- 

XX 

. 

XX 

• 

XX 

XX 

• 

XX 

• 

• 

XX 

• 

• 

XX 

XX 

XX 

XX 

. 

• 

XX 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

• 

• 

• 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

XX 

xx 

XX 

• 

• 

XX 

XX 

XX 

XX 

XX 

• 

* 

• 

. 

. 

X< 

• 

XX 

• 

XX 

XX 

XX 

• 

XX 

XX 

XX 

xx 

XX 

XX 

• 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

xx 

XX 

. 

XX 

• 

XX 

XX 

• 

XX 

XX 

• 

XX 

XX 

• 

« 

• 

• 

. 

XX 

• 

XX 

XX 

XX 

XX 

• 

XX 

XX 

. 

XX 

XX 

XX 

XX 

. 

XX 

• 

XX 

XX 

XX 

XX 

xx 

XX 

• 

XX 

. 

XX 

XX 

XX 

XX 

XX 

• 

• 

XX 

XX 

XX 

. 

. 

XX 

XX 

XX 

XX 

XX 

XX 

. 

XX 

XX 

xx 

XX 

XX 

• 

XX 

. 

XX 

XX 

. 

• 

XX 

. 

• 

Xx 

XX 

XX 

• 

• 

XX 

XX 

• 

XX 

• 

• 

• 

• 

XX 

• 

XX 

XX 

• 

• 

XX 

. 

XX 

• 

XX 

XX 

• 

• 

XX 

X X 

XX 

XX 

XX 

. 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

• 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

. 

• 

XX 

• 

XX 

. 

XX 

• 

. 

XX 

• 

• 

XX 

XX 

• 

XX 

XX 

XX 

• 

XX 

XX 

XX 

• 

• 

XX 

XX 

XX 

XX 

XX 

• 

XX 

XX 

XX 

XX 

XX 

XX 

XX 

T 

Tc  ST 

OF 

1 

X 1 

1 V JOE 

versus  1 

X 2 

>OE 

A T 

• 05 

1 SI 

[',Ntr  it 

l.ANC 

r 

Pfc  KCf  NT  Au*.  (Jb  NUN— /u»<0  POINTS  WITH  ..  I t.N  | » 1 t.  AN  T UttllWLNl'.  I 5 66*2 


72 


The  2 deg  x 2 kin  Mode 

The  final  comparison  in  this  group  was  made  between  the  1 deg  x 
1 km  map  and  the  2 deg  x 2 km  map  (Fig.  12) . The  correlation  between 
these  two  maps  is  also  quite  poor,  and  most  of  the  small  scale  detail 
has  been  washed  out  on  the  2 deg  x 2 km  map.  Host  of  the  features 
appear  to  have  been  shifted  outwards  and  clockwise  slightly,  and  many 
of  them  have  been  stretched  in  the  direction  of  antenna  rotation. 
Although  area  B still  shows  up  well,  area  A has  lost  its  identity. 

Many  of  the  dBZ  values  on  the  2 deg  x 2 km  map  are  5 dBZ  lower  than 
the  corresponding  values  on  the  1 deg  x 1 km  map;  however,  the  maximum 
dBZ  value  observed  is  still  40  dBZ. 

The  statistical  analysis  at  a — .05  (Figure  17)  shows  that  66.9 
percent  of  the  pairs  of  non-zero  data  points  are  significantly  dif- 
ferent from  one  another.  Once  again,  the  only  areas  where  the  points 
generally  do  not  differ  are  the  areas  of  lower  reflectivity.  Although 
small,  both  areas  A and  B can  be  seen  on  this  map. 

Table  4 summarizes  the  statistical  comparisons  of  all  four  modes 
at  each  of  the  three  levels  of  significance. 

Comparison  of  10  cm  Data  in  the  Near  Range 
The  zero-tilt  reflectivity  maps  for  the  10  cm  data  in  the  near 


range  that  were  recorded  simultaneously  with  the  3 cm  data  are  shown 
in  Figures  18  through  21.  Once  again,  the  1 deg  x 1 km  mode  was  used 
as  the  standard  for  comparison  in  all  the  analyses. 
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Fig.  17.  Statistical  comparison  of  the  1 deg  x 1 km  and 
2 dog  x 2 km  modes,  3 cm,  near  range.  Level  of  significance  is 

.05. 
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TABLE  4.  Summary  of  Statistical  Analyses  of  Collection  Modes,  3 cm 
Near  Range. 


The  2 deg  x 1 km  Mode 


The  first  comparison  made  was  between  the  1 deg  x 1 km  map  (Figure 


18)  and  the  2 deg  x 1 km  map  (Figure  19).  The  most  notable  feature  at 


first  glance  is  the  fact  that  the  dBZ  values  of  many  of  the  reflectiv- 


are  10  to  20  dBZ  higher  than  the  corresponding  features 


shown  on  the  3 cm  maps.  This  is  good  tangible  evidence  of  the  in- 


creased attenuation  by  precipitation  of  microwave  energy  at  wavelengths 


shorter  than  10  cm.  The  correlation  between  these  two  maps  is  fair 


and  the  2 deg  x 1 km  map  retains  most  of  the  features  of  the  1 deg  x 


1 km  map,  although  some  of  the  detail  is  smoothed  over.  Many  of  the 


low  reflectivity  features  have  been  shifted  slightly  in  the  direction 


of  antenna  rotation,  and  in  general,  the  extreme  low  reflectivity 


points  on  the  2 deg  x 1 km  map  tend  to  be  5 to  10  dBZ  higher  than  the 


corresponding  points  on  the  1 deg  x 1 km  map.  The  two  large  zero 


Percentage  of  non-zero  pairs  of  points 
with  significant  difference  at  signifi- 
cance level  a. 

1 deg  x 1 km 
mode 

compared  to 

a = .01 

a = .05 

a = .10 

2 deg  x 1 km 

25.4 

35.4 

43.2 

1 deg  x 2 km 

56.6 

66.2 

71.3 

reflectivity  areas,  A and  B,  show  up  well  on  the  1 deg  x 1 km  map 
again,  but  only  A is  well  defined  on  the  2 deg  x 1 km  map.  Area  B 
is  almost  completely  masked  by  a single  high  reflectivity  point  that 
is  located  at  the  front  edge  of  the  area,  relative  to  the  direction  of 
antenna  rotation.  The  highest  reflectivity  value  observed  on  both  maps 
is  50  dBZ. 

The  statistical  comparison  between  these  two  maps  (Figure  22)  at 
the  a = .05  significance  level  shows  that  45.2  percent  of  all  the  pairs 
of  non-zero  points  differ  from  one  another  significantly.  Most  of  the 
differences  occur  in  areas  of  high  reflectivity  gradients  and  small 
scale  high  reflectivity  features.  This  agrees  well  with  the  smoothing 
of  features  observed  in  the  2 deg  x 1 km  map.  Neither  area  A nor  B 
show  up  on  this  map. 


1 

i 


The  next  comparison  was  between  the  1 deg  x 1 km  map  and  the  1 deg 
x 2 km  map  (Figure  20).  Correlation  between  these  two  maps  is  poor. 
Most  of  the  detail  is  either  lost  or  significantly  altered,  and  only 
the  most  prominent  reflectivity  features  retain  their  identity.  There 
is  a noticeable  shift  radially  outwards  of  many  of  the  high  and  low 
reflectivity  features,  and  many  of  the  low  values  are  5 dBZ  higher 
than  those  on  the  1 deg  x 1 km  map.  Additionally,  many  of  the  high 
reflectivity  features  are  10  dBZ  higher  than  on  the  1 deg  x 1 km  map, 
as  is  the  maximum  observed  value  on  this  map  of  60  dBZ.  Once  again, 
the  zero  reflectivity  area  at  A shows  up  quite  well,  even  though  it 


m* 


is  shifted  outwards  almost  3 km,  while  that  at  B is  not  readily 
apparent. 
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Fig.  22.  Statistical,  comparison  of  the  1 deg  x 1 km  and 
2 deg  x 1 km  modes,  10  cm,  near  range.  Level  of  significance  is 
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The  statistical  analysis  of  these  two  maps  (Figure  23)  shows  that 
70.1  percent  of  all  the  non-zero  pairs  of  points  are  significantly 
different  from  one  another  at  the  a = .05  level  of  significance.  The 
only  areas  where  almost  all  of  the  points  do  not  differ  are  several  of 
the  areas  of  generally  low  reflectivity.  Neither  area  A nor  B show  up 
on  this  map. 

The  2 deg  x 2 km  Mode 

The  final  comparison  in  this  section  was  between  the  1 deg  x 1 km 
map  and  the  2 deg  x 2 km  map  (Figure  21) . The  correlation  between 
these  two  maps  is  also  very  poor.  There  is  a great  deal  of  smoothing 
of  detail  and  altering  of  the  shape,  size,  and  position  of  the  various 
reflectivity  features.  Once  again,  both  the  high  and  ].ou  reflectivity 
points  tend  to  be  5 to  10  dBZ  higher  on  this  map  than  on  the  1 deg  x 
1 km  map;  however,  the  maximum  observed  dBZ  value  on  both  maps  is 
50  dBZ.  Many  of  the  features  are  shifted  radially  outwards  one  or 
more  km,  and  there  is  also  some  elongation  of  features  in  the  direction 
of  antenna  rotation,  due  no  doubt  to  the  longer  integration  intervals 
in  both  range  and  azimuth.  The  zero  reflectivity  area  at  A is  very 
apparent,  but  once  again,  the  area  at  B has  been  lost. 

The  statistical  comparison  of  these  two  maps  (Figure  24)  shows 
that  72.9  percent  of  all  the  non-zero  pairs  of  points  were  signifi- 
cantly different  from  one  another  at  the  a = .05  level  of  signifi- 
cance. Once  again,  only  in  the  areas  of  generally  low  reflectivity 
was  there  any  substantial  correlation  between  the  data  points,  and 
neither  area  A nor  area  B was  apparent  on  this  map. 
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Fig.  23.  Statistical  comparison  of  the  1 deg  x 1 km  and 
1 deg  x 2 km  modes,  10  cm,  near  range.  Level  of  significance 
is  .05. 
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Fig.  24.  Statistical  comparison  of  the  1 deg  x 1 km  and 
2 deg  x 2 km  modes,  10  cm,  near  range.  Level  of  signifcance  is 


Finally,  a summary  of  the  statistical  comparisons  of  all  four 
modes  of  collection  at  each  of  the  three  levels  of  significance  is 
presented  in  Table  5. 

Comparison  of  3 cm  Data  in  the  Far  Range 

The  next  set  of  analyses  was  made  of  the  3 cm  data  collected  in 
all  four  modes  in  the  far  range.  Initially,  the  far  range  was  to  be 
distances  of  approximately  150  km  from  the  radar.  However,  there  were 
no'  radar  returns  of  any  significance  at  that  range,  so  the  area 
farthest  from  the  radar  that  had  any  clearly  defined  reflectivity 
features  was  chosen  instead,  and  that  is  the  area  presented  in 
Figures  25  through  28.  These  figures  cover  the  area  which  ranges  from 

TABLE  5.  Summary  of . Statistical  Analyses  of  Collection  Modes,  10  cm, 
hear  Range. 


Percentage  of  non-zero  pairs  of  points  with 
significant  difference  at  significance 
level  a. 

1 deg  x 1 km 
mode 

compared  to 

a = . 01 

a = .05 

a = .10 

2 deg  x 1 km 

36.5 

45.2 

51.5 

1 deg  x 2 km 

63.5 

70.1 

75.7 

2 deg  x 2 km 

64.9 

72.9 

76.9 

75  to  140  km  from  the  radar.  Specifically,  the  major  reflectivity 
feature  in  the  northeast  corner  of  the  maps  is  106  km  from  the  radar. 
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Fig.  27.  Three  cm  reflectivity  map  in  the  far  range,  1 deg  x 2 km 
>de. 
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Fig.  28.  Three  cm  reflectivity  map  in  the  far  range,  2 deg  x 2 km 
>de. 
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The  2 dee  x 1 km  Mode 


The  first  comparison  of  this  series  was  between  the  1 deg  x 1 km 

map  and  the  2 deg  x 1 km  map  (Fig.  26).  The  correlation  between  these 

two  maps  is  quite  good,  and  the  2 deg  x 1 km  map  retains  almost  all  of 

the  significant  features  of  the  1 deg  x 1 km  map.  There  is  some 

smoothing  of  fine  scale  details  and  a slight  shift  of  some  features 
in  the  direction  of  antenna  rotation,  as  well  as  a slight  elongation  of 
some  features  in  that  direction.  Additionally,  three  small,  adjacent 
reflectivity  features  on  the  1 deg  x 1 km  map  (location  A)  are  combined 
into  one  larger  area  on  the  2 deg  x 1 km  map,  and  the  extent  of  that 
area's  region  of  highest  reflectivity  is  smaller.  A different  situa- 
tion prevails  in  cell  B,  however,  where  the  maximum  reflectivity 
on  the  2 deg  x 1 km  map  is  30  dBZ,  whereas  on  the  1 deg  x 1 Ion  map  it 
is  only  25  dBZ. 

The  statistical  comparison  of  these  two  maps  (Figure  29)  indi- 
cates that  only  26.2  percent  of  the  non-zero  pairs  of  points  were  sig- 
nificantly different  from  one  another  at  the  a = .05  level  of  signifi- 
cance. This  pair  of  maps  corresponds  better  statistically  than  any  of 
the  other  pairs  of  maps  analyzed,  but  this  is  no  doubt  due,  at  least 
in  part,  to  the  relatively  low  reflectivity  values  that  comprise  most 
of  the  map.  The  areas  of  greatest  disagreement  are  cell  B and  area  A, 
as  could  be  expected,  where  almost  55  percent  of  the  points  differ 
significantly. 


the  1 deg  x 2 km  map  (Figure  27).  The  correlation  between  these  two 
maps  is  poor.  Most  of  the  detail  has  either  been  smoothed  over  or 
altered,  and  the  trailing  edges  of  most  of  the  reflectivity  features 
have  been  shifted  outwards.  Area  A and  cell  B have  been  combined  into 
one  elongated  feature,  and  once  again,  the  maximum  reflectivity  of 
this  area  is  5 dBZ  higher  than  on  the  1 deg  x 1 km  map. 

The  statistical  comparison  of  these  two  maps  (Figure  30)  shows 
that  48.5  percent  of  all  the  non-zero  pairs  of  points  vary  signifi- 
cantly from  one  another  at  the  a = .05  level  of  significance.  Once 
again,  the  region  of  greatest  difference  is  in  the  areas  A and  B, 
although  differences  are  widespread  throughout  the  rest  of  the  map, 
also. 

The  2 deg  x 2 km  Mode 

The  final  comparison  in  this  section  was  between  the  1 deg  x 1 km 
map  and  the  2 deg  x 2 km  map  (Figure  28) . The  correlation  between 
these  maps  is  also  poor.  Most  of  the  detail  is  lost  and  the  shape  of 
many  of  the  features  has  been  substantially  altered.  Many  of  the 
features,  particularly  the  low  reflectivity  features,  are  shifted 
slightly  outwards  and  in  the  direction  of  antenna  rotation.  Area  A 
and  cell  B have  been  combined  into  a feature  that  bears  little  re- 
semblance to  those  on  the  1 deg  x 1 km  map. 

The  statistical  comparison  of  these  two  maps  (Figure  31)  shows 
that  45.7  percent  of  the  non-zero  pairs  of  points  differ  significantly 
from  one  another  at  the  a = .05  significance  level.  The  percentage 
is  much  higher  in  the  region  of  area  A and  cell  B (almost  66  percent) , 
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Fig.  30.  Statistical  comparison  of  the  1 deg  x 1 km  and 
1 deg  x 2 km  modes,  3 cm,  far  range.  Level  of  significance  is 

.05. 
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Fig.  31.  Statistical  comparison  of  the  1 deg  x 1 km  and 
2 deg  x 2 km  modes,  3 cm,  far  range.  Level  of  significance  is 


but  the  differences  are  still  widespread  throughout  the  entire  map. 

Lastly,  the  summary  of  the  statistical  comparisons  of  all  four 
modes  at  all  three  levels  of  significance  is  presented  in  Table  6. 


TABLE  6.  Summary  of  Statistical  Analyses  of  Collection  Modes,  3 cm. 
Far  Range. 


Percentage  of  non-zero  pairs  of  points  with 
significant  differences  at  significance 
level  a 


1 deg  x 1 km 
mode 

compared  to 

a = .01 

£ 

II 

o 

Ln 

i 

a = .10 

2 deg  x 1 km 

17.7 

26.2 

31.9 

, 

1 deg  x 2 km 

37.3 

48.5 

55.9 

2 deg  x 2 km 

34.3 

45.7 

. 53.0 

Comparison  of  10  cm  Data  in  the  Far  Range 

The  zero-tilt  reflectivity  maps  for  the  10  cm  data,  recorded 
simultaneously  with  the  3 cm  data,  are  presented  in  Figures  32  through 
35.  As  before,  the  1 deg  x 1 1cm  map  was  used  as  the  standard  for 
comparison. 

The  2 deg  x 1 km  Mode 

The  first  comparison  in  this  set  was  between  the  1 deg  x 1 km  map 
(Figure  32)  and  the  2 deg  x 1 km  map  (Figure  33).  The  correlation  be- 
tween these  two  maps  is  fairly  good,  and  most  of  tv a detail  and  shape 
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Fig.  35.  Ten  cm  reflectivity  map  in  the  far  range,  2 deg  x 2 km 

ide. 
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of  the  features  on  the  1 deg  x 1 kin  map  are  carried  over  to  the 
2 deg  x 1 kin  map.  The  correspondence  of  cell  B on  the  two  maps  is 
particularly  good,  even  though  the  point  of  maximum  reflectivity  is 
elongated  somewhat  in  the  direction  of  antenna  rotation.  The  high  and 
low  reflectivity  values  match  quite  well  in  all  regions  of  the  two 
maps,  and  the  maximum  observed  reflectivity  value  on  both  maps  is  35 
dBZ,  compared  to  30  dBZ  on  the  corresponding  3 cm  data  maps. 

The  statistical  comparison  of  these  two  maps  (Figure  36)  indi- 
cates that,  at  the  a = .05  significance  level,  29.0  percent  of  all  the 
non-zero  pairs  of  points  differ  from  one  another  significantly.  Most 
of  the  differences  occur  in  regions  of  high  reflectivity  or  reflec- 
tivity gradient,  as  is  to  be  expected. 

The  1 deg  x 2 km  Mode 

The  next  pair  of  maps  analyzed  consisted  of  the  1 deg  x 1 km  map 
and  the  1 deg  x 2 km  map  (Figure  34) . The  correlation  between  these 

two  maps  is  poor,  as  much  of  the  detail  has  been  lost  and  the  shapes  of 

many  reflectivity  features  have  been  altered  substantially.  The  trail- 
ing edges  of  most  and  the  leading  edges  of  many  reflectivity  features 
have  been  shifted  radially  outward  two  km  or  more,  and  many  features 
with  small,  closed  isopleths  have  been  expanded  in  size.  Cell  B has 
lost  some  of  its  identity,  and  the  maximum  reflectivity  value  of  that 
cell  is  only  30  dBZ,  compared  to  35  dBZ  on  the  1 deg  x 1 km  map. 

The  statistical  analysis  of  these  two  maps  (Figure  37)  shows  that 

52.2  percent  of  the  non-zero  pairs  of  points  differed  significantly 
from  one  another  at  the  a = .05  level  of  significance.  Although  the 
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Fig.  37.  Statistical  comparison  of  the  1 deg  x 1 km  and 
1 deg  x 2 km  modes,  10  cm,  far  range.  Level  of  significance  is 
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differences  are  more  concentrated  in  the  regions  of  high  reflectivity 
gradient  and  dBZ  values,  the  only  regions  where  differences  are  not 
widespread  are  those  broad  areas  of  generally  low  reflectivity. 
Elsewhere,  differences  are  the  rule,  rather  than  the  exception. 

The  2 deg  x 2 km  Mode 

The  final  comparison  was  between  the  1 deg  x 1 km  map  and  the 
2 deg  x 2 km  map  (Figure  35) . The  correlation  between  these  maps  is 
quite  poor.  Most  of  the.  detail  is  either  lost  or  significantly 
altered,  and  many  of  the  features  are  shifted  outwards.  Reflectivity 
values  over  the  entire  map  tend  to  be  5 dBZ  lower  than  on  the  1 deg 
x 1 lem  map,  including  the  maximum  observed  values.  Cell  B on  the  2 deg 
x 2 km  map  bears  little  or  no  resemblance  to  its  counterpart  on  the 
1 deg  x 1 km  map. 

The  statistical  comparison  of  these  two  maps  (Figure  38)  shows 
that  50.4  percent  of  the  non-zero  pairs  of  points  differ  from  one 
another  significantly  at  the  a = .05  level  of  significance.  The  only 
regions  where  differences  are  not  the  rule  are  the  broad  areas  of  low 
reflectivity  where  the  values  are  relatively  constant. 

Finally,  Table  7 presents  a summary  of  the  statistical  compari- 
sons of  all  four  modes  at  all  three  levels  of  significance.  Addi- 
tionally, Table  8 presents  the  same  summary,  but  averaged  over  all  the 
data,  both  3 and  10  cm  and  near  and  far  ranges.  The  analyses  con- 
sistently show  that  while  the  2 deg  x 1 km  mode  data  correspond  fairly 
well  to  the  1 deg  x 1 km  mode  data,  neither  of  the  other  two  modes  do. 
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Fig.  38.  Statistical  comparison  of  the  1 deg  x 1 km  and 
2 deg  x 2 km  modes,  10  cm,  far  range.  Level  of  significance 
is  .05. 


TABLE  7.  Summary  of  Statistical  Analyses  of  Collection  Modes,  10  cm 
Far  Range. 


Percentage  of  non-zero  pairs  of  points  with 
significant  difference  at  significance 
level  a. 


TABLE  8.  Averages  of  Statistical  Analyses  of  Collection  Modes,  Taken 
Over  Both  Wavelengths  and  Ranges. 


Percentage  of  non-zero  pairs  of  points  with 
significant  difference  at  significance 
level  a. 


1 deg  x 1 km 
mode 

compared  to 


.01  ct  = .05 


2 deg  x 1 km 

25.1 

34.0 

40.4 

1 deg  x 2 km 

49.5 

59.3 

65.0 

2 deg  x 2 km 

49.5 

59.0 

64.6 
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Summary 

The  results  of  the  analyses  of  the  various  data  collection  modes 
are  somewhat  surprising  in  that  they  contradict  the  results  that  were 
initially  anticipated.  Specifically,  it  was  expected  that  the  2 km 
modes  would  produce  better  results  than  the  2 deg  modes,  especially  at 
longer  ranges.  However,  such  has  not  proven  to  be  the  case.  The  ex- 
planation probably  lies  in  the  method  of  integration  employed  along 
the  data  radial.  The  2 km  mode  uses  the  average  of  eight  250  m incre- 
ments of  return  power,  whereas  the  1 km  mode  uses  only  four.  This 
averaging  over  2 km  may  be  enough  to  substantially  alter  or  mask  the 
shape  and  intensity  of  any  but  the  most  prominent  reflectivity  fea- 
tures. 

The  integration  across  the  azimuth  increment  is  slightly  dif- 
ferent, however,  in  that  it  integrates  the  last  31  average  values  for 
that  range  bin,  either  over  a space  of  1 or  2 deg  of  azimuth.  Con- 
sequently, in  the  2 deg  x 1 km  mode, it  uses  the  last  31  one-km  range 
bin  averages,  each  of  which  is  more  accurate  than  the  2 km  range  bin 
averages  that  are  used  by  the  1 deg  x 2 km  and  2 deg  x 2 km  modes.  The 
resultant  differences  between  the  modes  when  compared  to  the  1 deg  x 
1 km  mode  are  readily  apparent  in  Table  8. 

Data  Processing  Requirements 

For  purposes  of  this  discussion,  the  data  processing  requirements 
may  be  considered  to  be  computer  execution  time,  amount  of  magnetic 
tape  required,  and  the  time  required  to  complete  one  scan  sequence. 
These  prrameters  are  discussed  briefly  below,  and  then  summarized  for 


The  execution  time  required  by  the  computer  to  process  the  data 
varies  substantially  from  mode  to  mode.  The  execution  time  for  the 
1 deg  x 1 km  mode  is  the  lowest,  followed  closely  by  the  2 deg  x 1 km 
mode.  The  times  shown  in  Table  9 represent  the  time  required  to  pro- 
duce one  0-deg  tilt  reflectivity  map,  averaged  over  both  the  3 and  10 
cm  data  at  both  the  near  and  far  ranges. 

These  results  also  appear  to  conflict  with  the  execution  times 
that  were  initially  expected.  However,  this  is  most  likely  due  to  the 
fact  that  the  simplest  mode,  in  terms  of  computer  execution,  (1  deg  x 
1 km)  was  run  first,  and  then  the  program  was  modified  to  accommodate 
the  next  mode  to  be  run,  etc.,  until  all  the  modes  had  been  run.  Each 
modification  of  the  program  introduced  numerous  additional  executable 
steps,  many  of  which  fell  inside  computational  loops  which  required  a 
great  number  of  iterations.  Consequently,  even  though  the  number  of 
bits  of  data  to  be  processed  decreased,  the  number  of  executable  state- 
ments increased,  and  hence  the  apparent  discrepancy  in  execution  times. 
This  problem  should  not  exist  when  processing  data  with  the  RADAR  pro- 
gram  in  its  general  form  as  it  appears  in  Appendix  A. 

Magnetic  Tape  Requirements 


The  length  of  each  logical  record  recorded  on  the  tape  is  the 
same  (1744  bytes),  regardless  of  the  mode  in  use,  so  the  difference 
comes  from  the  number  of  degrees  of  azimuth  that  are  covered  by  each 
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logical  record.  Knowing  this,  and  assuming  that  each  scan  sequence 
consists  of  ten  revolutions  of  the  antenna,  it  is  easy  to  determine  the 
total  number  of  logical  records  required  to  record  one  complete  scan 
sequence.  Both  the  2 deg  x 1 km  and  the  1 deg  x 2 km  modes  require 
only  one  half  as  much  tape  as  the  1 deg  x 1 km  mode,  and  the  2 deg  x 2 
km  mode  requires  only  one  fourth  as  much. 

Time  Required  for  One  Scan  Sequence 

The  time  required  for  one  complete  scan  sequence  is  dependent 
upon  the  radar  antenna  rpm  which  should  correspond  to  the  antenna  rpm 
required  by  the  DVIP  to  maintain  the  desired  data  accuracy  of  1.0  dB. 
This  assumes  that  a scan  sequence  consists  of  ten  antenna  revolutions, 
each  at  a different  elevation  angle.  Note  that  the  time  required  for 
the  1 deg  modes  (using  31  time  samples)  is  ten  minutes,  whereas  it  is 
five  minutes  for  the  2 deg  modes.  If  only  15  time  samples  are  used, 
then  the  times  decrease  accordingly.  However,  it  should  be  noted 
that  operation  of  the  antenna  at  speeds  in  excess  of  3.00  rpm  is  un- 
desirable due  to  the  added  physical  stresses  placed  on  the  antenna 


system. 


■ — — B u * 


109 


CHAPTER  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

The  primary  objective  of  this  investigation  was  to  examine  and 
compare  dual  wavelength  digital  radar  data  collected  near-simultaneously 
in  the  four  modes  of  collection  possible  by  the  TAMU  weather  radar 
system,  in  an  effort  to  determine  whether  there  is  an  optimum  mode  of 
data  collection  for  operational  use  of  the  radar  system.  This  deter- 
mination was  based  on  such  considerations  as  resolution  of  data,  re- 
tention of  detail,  and  smoothing  of  fine  scale  meteorological  features. 
These  were  weighed  against  such  considerations  as  range  of  the  storm 
from  the  radar,  wavelength  of  the  radar,  and  data  processing  require- 
ments. Digital  data  collected  from  the  3 and  10  cm  radars  of  the 
Department  of  Meteorology  at  TAMU  were  processed  by  a modified 
version  of  Sieland's  (1977)  computer  program,  and  then  displayed  as 
contoured  zero-deg  tilt  reflectivity  maps.  These  maps  provided  the 
basis  for  the  comparisons  and  evaluations  of  the  different  data  col- 
lection modes. 

The  secondary  objective  of  this  investigation  was  to  provide  a 
detailed  description  of  the  composition  of  the  TAMU  digital  weather 
radar  system  itself. 

This  investigation  has  led  to  the  following  conclusions. 

1.  A significant  difference  does  exist  between  the  zero-tilt 
reflectivity  map  produced  by  data  collected  in  the  1 deg  x 1 km  mode 
and  the  maps  produced  by  data  collected  in  the  other  three  modes. 


•j 
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Assuming  the  1 deg  x 1 km  map  to  be  the  reference  standard,  the  2 deg 
x 1 km  map  most  closely  approximates  the  standard,  in  terms  of  reten- 
tion of  detail,  magnitude  of  reflectivity  values,  and  shifting  or 
altering  of  reflectivity  features.  The  remaining  two  modes,  1 deg  x 2 
km  and  2 deg  x 2 km,  resemble  the  1 deg  x 1 km  map  much  less  closely. 

2.  For  all  cases  wherein  the  digital  radar  data  are  to  be  used 
for  research  or  other  scientific  applications,  where  resolution  and 
fine  scale  detail  are  of  utmost  importance,  but  time  is  not,  the  1 deg 
x 1 km  mode  should  be  utilized.  This  mode  affords  the  maximum  amount 
of  detail  and  resolution,  but  it  also  requires  twice  as  much  time  and 
magnetic  tape  as  the  2 deg  x 1 km  mode. 

3.  In  those  cases  wherein  the  application  of  the  digital  radar 
data  is  to  be  an  operational  one,  and  time  is  of  much  more  importance, 
the  2 deg  x 1 km  mode  should  prove  to  be  adequate.  This  mode  retains 
most  of  the  detail  and  significant  features  of  the  1 deg  x 1 km  mode, 
while  completing  a scan  in  only  one  half  the  time  required  for  the 
latter,  and  using  only  one  half  as  much  magnetic  tape. 

4.  The  remaining  two  modes,  the  1 deg  x 2 km  and  the  2 deg  x 2 km 
modes,  appear  to  suffer  enough  loss  of  detail  and  introduce  enough 
distortion  to  be  of  limited  usefulness  at  this  point.  However,  further 
research,  perhaps  utilizing  data  from  complete  scan  sequences  taken 


through  storms  with  pronounced  severe  weather  features,  may  negate  this 


conclusion. 
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Recommenda t ions 

The  results  of  this  investigation  suggest  the  following  recom- 
mendations concerning  the  processing  of  digital  radar  data. 

1.  The  RADAR  computer  program  should  be  revised  so  that  it  packs 
data  with  return  power  values  that  are  below  the  radar's  minimum 
detectable  signal  in  order  to  substantially  reduce  the  amount  of  com- 
puter memory  utilized  and  the  amount  of  time  required  to  read  the 
various  arrays.  Additionally,  the  program  should  be  streamlined, 
eliminating  redundant  or  unnecessary  steps  that  increase  execution 
time. 

2.  The  program  should  be  revised  to  incorporate  a variable  sized 
box  for  the  area  to  be  analyzed,  instead  of  the  current  100  deg  x 100 
km  box.  The  current  memory  size  could  be  retained,  while  making  the 
lengths  of  the  sides  of  the  box  flexible. 

3.  The  program  should  be  permanently  stored  as  a library  program 
that  could  be  called  up  and  initialized  with  the  appropriate  parameters, 
instead  of  manually  reading  a deck  of  more  than  1,000  cards  for  every 
program  run,  as  is  now  the  practice. 

4.  The  revised  program  should  then  be  programmed  and  tested  on  a 
mini-computer  to  determine  the  time  required  to  produce  the  PVSZ  maps, 
and  ultimately,  the  feasibility  of  the  operational  implementation  of 
this  or  similar  type  programs. 

5.  A comparison  of  data  integrated  over  15  time  samples  and  data 
integrated  over  31  time  samples  should  be  conducted  to  determine 
whether  it  is  feasible  to  use  data  collected  in  the  15  time  sample  mode, 
and  hence  decrease  substantially  the  time  required  to  complete  an 
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entire  scan  sequence. 

The  following  additional  recommendations  are  also  made. 

1.  A radar  standard  operating  procedure  should  be  established 
to  insure  that  all  radar  data  are  recorded  at  the  antenna  rpm  and 
time  sample  rate  required  to  produce  the  same  level  of  accuracy  (1.0 
dB)  for  all  modes  of  collection. 

i ; 

2.  An  analysis  should  be  conducted  of  the  rainfall  rate  of  a 
storm  derived  from  the  VIL  values  produced  by  data  obtained  in  the 
different  modes  of  collection.  This  would  require  a relatively  steady 
state  storm,  as  the  computation  of  VIL  requires  that  a complete  scan 
sequence  be  recorded. 

3.  A continuing  effort  should  be  undertaken  to  observe  and  study 
tornadic  storms  in  this  region  of  central  Texas.  Data  from  such  a 
storm  should  be  recorded  in  all  four  modes  of  collection  to  test  the 
validity  of  the  conclusions  of  this  investigation. 
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APPENDIX  A 


This  appendix  contains  a commented  copy  of  the  computer  program 
(RADAR)  that  was  used  to  process  the  digital  data  recorded  by  the 
radar  system. 
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THIS  PROGRAM  USES  DIGITIZED  RADAR  DATA  FROM  THE  TEXAS  ACM  RAOAR 
SYSTEM.  THESE  DATA.  WRITTEN  ON  MAGNETIC  TAPE.  ARE  READ  ANO  USEO 
TO  COMPUTE  AND  PROOUCE  THE  FOLLOWING: 

1.  ZERO  DEGREE  TILT  REFLECTIVITY  MAPS 

2.  PARTIAL  VERTICALLY  SUMMEO  REFLECTI V I TY I Z I MAPS  (PVSZI  FOR 

3 LAYERS  ISFC  - 15  KFT.  15  KF T - 35  KFT.  C 35  KFT  - 50  KFT) 

3.  VERTICALLY  INTEGRATED  LIQUID  WATER  CONTENT  (VIL)  MAP 

THE  program  will  operate  on  oata  recorded  in  all  four  mooes 

POSSIBLE  BY  THIS  RADAR  SYSTEM:  I DEG  BY  1 KM.  1 DEG  BY  2 KM. 

2 DEG  BY  1 KM.  ANO  2 DEG  BY  2 KM.  HOWEVER.  IT  WILL  ONLY  PROCESS 
OATA  FOR  ONE  RADAR  13  CM  OR  10  CM)  AT  A TIME.  To  PROCESS  DATA 
FROM  THE  OTHER  RADAR.  YOU  MUST  CHANGE  THE  RADAR  CONSTANT  IWAVEC). 
THE  RADAR  CALIBRATION  DATA  CARO.  AND  THE  JCL  FOR  THE  INPUT  TAPE. 

definitions: 


ARRAYS 

ZW  - CONTAINS  RAW  REFLECTIVITY  DATA  VIA  CODED  AZIMUTH  ANO  RANGE 
Z - CONTAINS  A PVSZ  OR  VIL  ARRAY  AFTER  IT  IS  RETURNED  FROM  OD2. 

THE  ARRAY  WILL  BE  IN  X.  V COORDINATES 
PVSZI.  2.  3 - THE  THREE  PVSZ  LAYERS 
VIZ  - THE  VIL  ARRAY 

A - THE  MATRIX  USEO  FOR  THE  LAGRANCIAN  INTERPOLATION 
IDIVt.  2.  3 - The  DIVISORS  USED  FOR  OBTAINING  THE  AVERAGE 

REFLECTIVITY  IN  A LAYER.  THIS  CONSISTS  OF  THE  NUMBER  OF 

data  radials  in  a given  layer  at  a given  range,  see 
subroutine  divioe. 

ZP  - A DUMMY  ARRAY  USEO  FO«  THE  002  SUBROUTINE  WHICH  GIVES  THE 
PVSZ  ARRAYS  BACK  TO  THE  MAIN  PROGRAM  IN  X.  Y COORDINATES 
ICLASS  - CONTAINS  RAW  DIGIT  I ZED  REFLECTIVITY  DATA 
CALIB  - CONTAINS  THE  EQUIVALENT  DBM  VALUES  FOR  256  DIGITAL 

integers  numbered  from  o to  255 

CORECT  - CONTAINS  THE  EARTH  CURVATURE  CORRECTION  SASEO  ON  RANGE 
HT I • HT 2 - THE  HEIGHT  ARRAYS  TO  COMPUTE  DH  IN  THE  VIL  ROUTINE 
PZ  - CONTAINS  THE  REFLECTIVITY  FACTORS  FOR  THE  PREVIOUS 
TILT  ANGLE  USED  FOR  COMPUTATION  OF  VIL 

FOR  OATA  CAROS 

The  first  DATA  CARO  read  contains  THE  RAOAR  calibration  data. 

OIFFERENT  FOR  EACH  RAOAR. 

the  next  data  card  contains  the  initializing  information  listed 

BELOW,  IN  THE  ORDER  GIVEN,  STARTING  IN  COL.  I. 

example:  10200226900102201 9402210061 85A30 

each  lelment  listeo  is  contained  in  parentheses  below. 


ILEFT  - DISTANCE  EAST  (-)  OR  WEST  If)  OF  TAMU  TO  THE  LOWER  LEFT 
CORNER  OF  THE  100  X 100  KM  GRID.  PLUS  TWO. 

FOR  EXAMPLE.  80  WEST  = ♦ 82 « 80  EAST  = -78  EX.  1102) 

JDOWN  - DISTANCE  NORTH  (-)  OR  SOUTH  (♦)  OF  TAMU  TO  THE  LOWER 
LEFT  CORNER  OF  THE  100  X 100  KM  GRID.  PLUS  TWO 
FOR  EXAMPLE.  IS  SOUTH  = *17,  15  NORTH  = -13.  EX.  (002) 

KAZ  - DESIRED  STARTING  AZIMUTH  TO  LOAD  THE  GRID  BOX-  MOVING 
CLOCKWISE  EX.  (269). 

LASTAZ  - THE  LAST  OATA  AZIMUTH  REQUIRED  TO  LOAD  THE  GRID  BOX 
EX.  (001) 

NTLT  - THE  LAST  tilt  angle  incremented  by  I degree  IN  THE  TILT 
SEOUENCE  EX.  (02) 

KRANGE  - THE  OISTANCE  TO  THE  NEAREST  POINT  IN  THE  GRID  BOX 
EX.  (20) 
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C IRIT  - FILE  ONTO  WHICH  DATA  IS  WRITTEN.  DEFINED  IN  JCL  EX.  <101 

C 1 HR  - FILE  FROM  WHICH  DATA  IS  READ.  OEF INED  IN  JCL  EX.  (08) 

C NENOT  - END  TIME  OF  THE  SCAN  SEOUENCE  EX.  (185430) 

C 

c 

c •*»•****•  note  ••*•**••**• 

C KAZ  AND  LASTAZ  MUST  BE  ODD  NUMBERS 

C KRANGE  MUST  OE  AN  EVEN  NUMBER 

C *••**•**»♦***••*••••* 

c 

C OTHER  variables 

C WAV  EC  - LOG  OF  RADAR  CONSTANT  IN  DBM  M**3  MM**<- 6)  KM*«(-2| 

C XI l AND  XXII  - CONSTANTS  USED  FOR  CODING  THE  AZIMUTH  PARAMETER 

C OF  THE  ARRAYS 

C IAZ1  At.-D  IAZ2  - AZIMUTHS  AS  RE  AO  FROM  DATA  TAPE  (TWO  RAOIALS 

C OF  I KM  DATA  AND  FOUR  RAOIALS  OF  2 KM  DATA  ARE  CONTAINED 

C IN  EACH  RECORD) 

C LVL  - A COUNTER  USED  WHEN  PLOTTING  MAPS  TO  INDICATE  WHICH  LEVEL 

C IS  BEING  CALCULATED 

C K - AZIMUTH  PARAMETER  IN  COOED  FORM 

C MM  - COUNTER  USED  FOR  MAP  CAPTIONS 

C PTILT  - A REFERENCE  LEVEL  USED  TO  MONITOR  THE  ELEVATION  ANGLE 

C 

c 

c JCL  - JOB  CONTROL  LANGUAGE  ***** 

c with  this  jcl.  the  program  runs  in  «fortran  h extended*. 

c the  jcl  will  require  some  changes,  but  none  significant,  for 

c EACH  PROGRAM  RUN.  NOTE  THAT  THE  FIRST  CHARACTER  ON  ALL  JCL 

C CARDS  (COL.  I)  SHOULD  BE  A •/*,  REPRODUCED  HERE  AS  AN  *X*. 

C 

X/RADAR  JOB  ( X9 99. 006D, 0 09. 0 l 5 . XX ) . • YOUR  NAME  • 

C 

C 009  IS  TIME  |N  MINUTES.  015  IS  LINES  OF  PRINTED  OUTPUT  IN 

C THOUSANDS.  DON*  T CHANGE  THESE 

•c 


X4PASSWCRD  XXXX 

XASETUP 

C 

C THIS  TELLS  THE  MACHINE  NOT  TO  RUN  UNTIL  THE 

C TAPES  HAVE  BEEN  MOUNTED. 

C 

X4LEVEL  X 

X4J0BPARM  X = 0.R  = 51 2 • T— 9 • T A PE 9=  2 

X/  EXEC  FORTXCLG,GOREGN=5l2K.FXLNSPC=*3120.(lOO.lO)* 

X/FORT • SYS I N DO  * 

C 

C DO  NOT  CHANGE  THESE  LAST  THREE  CARDS 

C 

c 

c 4«44«4«*«  THE  MAIN  PROGRAM  CARO  DECK  BEGINS  HERE 
C 

c 

DIMENSION  ICLASS(200).CALIB(255) 

COMMON/ BLOCK A/ ZWI 100.IOO).PV1Z1(I00.100).PVIZ2(100.IOO).PVIZ3(IOO. 
I 100  > .VIZI  » 00. too) .KRANGE. A<  4.4) .101 VI ( ISO) .IDI V2( 1 50 ) . I DI V3( I SO ) . 
2ZP(  100.200) . PZ I I 00.100) .Z( 51 .51 ). CORE CT( ISO). ICOUNT.HTK 1001 . 

3HT2( 100) • I ADD. IADDR 
LOGICAL*!  LOGCI I 7 20  I 
INTEGER*2  I HWORO(  B50 ) 

EOul VALENCE! LOGCI  I ) , IHWORD1 I I I 
DIMENSION  IEL2I3) 

OATA  I ZRO/ • I •/, I ONE/* 2* /.NINE/* 9*/. IHW0RD/B60*0/. ICHECK/* 57* /, I ZER 


noon  orif>r»r»  oonnnn  n non  -nnnnnnnn  nnnnn 
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*o/*o*/ 

DATA  l CLASS/200* 0/. CAL IB/255* 0.0/ 

c 

VAVEC— 6 .0 

THIS  VALUE  MUST  BE  CHANCED  TO  MATCH  THE  RAOAR  THE  DATA  CAME 
FROM.  OR  WHENEVER  A NEW  RADAR  CONSTANT  IS  COMPUTED  FOR 
THE  RAO ARS 

CALL  LOAO(CALIB) 

THIS  RECONSTRUCTS  THE  CALIBRATION  CURVE  OF  THE  RADAR  AND 
CONVERTS  THE  256  DIGITAL  INTEGERS  TO  EQUIVALENT  DBM  VALUES. 

IT  ALSO  PRINTS  OUT  THE  ENTIRE  CALI  B ARRAY 

RE AO  OATA  CARO  FOR  THE  GRID  BOX.  STARTING  TIME.  OATA  FILES. 

AND  PRINT  THEM  OUT  BELOW  THE  CAL K B ARRAY. 

REA0I5. 10O.ENO=999991ILEFT, JOOWN. KAZ.LAST AZ. NTLT. KRANGE. ISTRT. IRIT 
I . IRR.NENOT 

100  FORMAT  <413.212. 16,212.16} 

WRITEI6.il II  K AZ, LAST AZ. KRANGE. ILEFT. JOOWN 
111  FORMAT  I • 0 • . • 1ST  AZ  = *.I3.*  LAST  AZ  = *.13.*  KRANGE  = ‘.I3.*  LEFT 
C = • . I 3.*  DOWN  = • . I 31 
WRITE! 6.10001 IRIT. I RR 

1000  FORMAT ( • 0 * , • WR I TE  ON  FILE  •.12,*.  READ  FROM  TAPE  *.I2) 

INITIALIZE  OR  OEFINE  COUNTERS 

NCOUNT=0 
ICOUNTsO 
K! I=KAZ-I 
KKI I*( 361-KAZ) 

PTILT=0.0 

CALL  DIVIDE  I NTLT ) 

01 V I OE  IS  A SUBROUTINE  WHICH  CALCULATES  THE  VALUES  OF  THE 
DIVISORS  NEEOEO  TO  OBTAIN  A MEAN  REFLECTIVITY  IN  A LAYER. 

READ  OATA  FROM  FILE  IRR 

15  READ 1 1 RR .177. ERR— 090. E NO  =99 99 9 I ID ATE. IT IME , ID . I DS . I R I . I SA. I RE. I PE . 
* I AZ I . I ELI • IEL2 • ILOGC 1 I I . I =2 . 1 72 0.2  I 
177  FORMAT ( A3. 1 I6.3X.6A1 . 13. A3.T22.3A1.4! 2I5AI ) I 


TEST  TIME  AGAINST  ISTRT.  IF  NOT  AT  STARTING  TIME.  REAO 
MORE  DATA.  IF  TIME  HAS  EXCEEDED  NENDT • BEGIN  PRINTING 
OUT  THE  REMAINING  MAPS. 


IFI ITIME.LT. ISTRT IGO  TO  IS 
IF! ITIME.GT. NENDT I GO  TO  99 
INOEXaO 
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IAZ2  * I AZ I «•  I A00 
I AZ3»I  AZ2J IAOO 
I AZ4«I AZ3  + I ADD 

TEST  AZIMUTH  TO  SEE  IF  I AZ  IS  BETWEEN  KAZ  AND  LASTAZ. 


IFIKAZ.LT.LASTAZI  GOT048 

IF| IAZI.LT.KAZ.AND.I AZ1 .GT. LASTAZ I G0T058 
IFiIAZl.LT.laOt  G0T05I 
GO  TO  49 

48  IFIIAZI.LT. KAZ.OR.IAZI.GT.LASTAZI  G0T058 

49  I Is  I AZ 1— KI  | 

CO  TO  S2 

51  tI~lAZI*KKII 

52  CONTINUE 

CALCULATE  TILT  AND  TEST  AGAINST  PTILT.  IF  TILT  IS  MUCH  MOPE 
THAN  PTILT.  THEN  THE  SCAN  HAS  BEEN  INCREMENTED  UPWARDS  ANO  THE 
PROGRAM  PRINTS  OUT  THE  ZERO  LEVEL  MAP.  THEN  PTILT  IS  ADJUSTED 
UPWARDS  TO  MATCH  TILT.  AND  THE  PROCESSING  CONTINUES.  IF  TILT 
IS  MUCH  LESS  THAN  PTILT.  THEN  THE  PROGRAM  ASSUMES  A NEW  SCAN 
HAS  BEEN  READ.  AND  IT  TERMINATES  PROCESSING  THE  DATA.  AND 
PRINTS  OUT  ALL  THE  MAPS. 

I EL  21 1 1 • 1 2 1 . AND  <3>  ARE  THE  TENS,  UNITS.  AND  TENTHS  DIGITS 
of  the  elevation  angle. 

ATOI  IS  A SUBROUTINE  THAT  CHECKS  TO  ENSURE  THAT  THE  DIGITS  OF 
t EL  1 WERE  CORRECTLY  CONVERTED  I NTO/FROM  BINARY  FORM.  A NON- 
ZERO RESULT  RETURNEO  FROM  THIS  TEST  ENOS  THE  PROGRAM 

IFI IEL2I I > .GT. NINE. OR. IEL2 1 2) . GT.NINE .OR. IEL2<  3 I . GT .N INE » GOTO 16 
CALL  ATOI ( IELI .3.0. IEL.NTEST I 
IF<  NTEST.NE .01 GOT 099 99 
TlLTztEL/10.0 

THIS  ENSURES  THAT  TILTdELl 


& CONTINUE 
INOEXst 

IFITILT-PT ILT .GE .0.4 1 GOT 07 6 
IF  ITILT-PTILT.LT.-0.51GO  TO  99 


PROCESSING  OF  THE  FIRST  RAOIAL  OF  DATA  BEGINS  HERE.  TEST  FOR 
RANGE  INTERVAL  IIRIi  ANO  LOAD  THE  ARRAYS  ICLASS  WITH  ISO  VALUES 
OF  REFLECTIVITY  DATA  FOR  THE  1 KM  INTERVAL  ANO  90  VALUES  FOR 
THE  2 KM  INTERVAL. 

THERE  ARE  860  REFLECTIVITY  VALUES  STOREO  IN  ONE  RECORD.  THIS 
REPRESENTS  2 RADIALS  OF  DATA  IN  THE  I KM  MODE  AND  4 RADI ALS 
IN  THE  2 KM  MODE.  IN  THE  I KM  MODE.  IHWOMOl I-430)=1ST  RADIAL. 

I 431—8601 *2N0  RADIAL.  IN  THE  2 KM  MODE.  I HWOROI 1 -215 1 = I ST. 

121 6—4301 —2ND.  (431-645 I-3RD.  AND  ( 646-860 1 *4TM. 


86 

C 

C 


600 


c 


CONTINUE 

IFI IRI • EO. I ONE  I GO  TO  5 

DO  600  1=1 .180 
I CLASS!  I»20I-IHWORD( I I 
GO  TO  6 


i 


9 DO  601  I =2 . 1 80. 2 


uuuug  uuuuuu  - ouou u u o u 
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601  {CLASS!  I T2  0 ) = I HWORDI 1/2) 

THIS  OO  LOOP  COMPUTES  THE  REFLECTIVITIES  FOR  100  RANGES  FROM 
THE  DIGITIZED  DATA  IN  ICLASS.  THE  RANGE  IS  I ♦ KRANGE.  FOR 
AN  EXPLANATION  OF  THE  EQUATION  FOR  ZWtII.II.  SEE  CHAPTER  II* 

6 OO  10  I=IAOOR. 100. IAOOR 
Kc I* KRANGE 
ZWIII, I>=0.0 
ARNG=FLOAT(  K ) 

IF  I ICLASSIKI.LE.3AIGOT010 

Z WI 11,11=1 0.**I-0. 1*CALI B( 1 CLASS1K1 > «■  2 . » ALOG  I 0< ARNG I ♦ WAVEC) 

IFI TILT.LT.0.5IGOTO31 t 
GOTO  10 

311  ZINTRP*ZWI  II  .1) 

IF  THE  TILT  ANGLE  IS  LESS  THAN  0.5  DEGREE.  THE  ZW  DATA  ARE 
AOOEO  INTO  THE  LOWER  LATER  PVSZ  ARRAT,  OTHERWISE  THEY  ARE  NOT. 
ZINTRP  IS  THE  INTERPOLATED  VALUE  OF  THE  REFLECTIVITY  AND  IT  IS 
ALWAYS  DIVIDED  BY  THE  APPROPRIATE  DIVISOR  BEFORE  SUMMATION. 

IFIZINTRP.LE.t .0) GOTO 10 

PVIZI ( II. I l=PV III ( II . I »♦  ZINTRP/IOI VI  I I) 

PZ( II . I )=Z INTRP 
0 CONTINUE 

IF  I TILT. LE. 0.51  G0T058 
CALL  INTRP2ITlLT.il  1 

THE  I NT RP 2 SUBROUTINE  USES  THE  TILT  ANGLE  AND  THE  CODED  AZIMUTH 
TO  ACCOMPLISH  THE  LAGRANGIAN  INTERPOLATION  ALONG  THE  RAOIAL  OF 
DATA  FOR  100  KM  OF  RANGE. 

58  CONTINUE 

THE  SECOND  RAOIAL  OF  DATA  IS  CONVERTED  AND  PROCESSED  AS  ABOVE. 
IF  IKAZ.LT.LASTAZ)  GO  TO  68 

IFIIAZ2.LT. KAZ. AND. IAZ2.Gr.LASTAZ)G0  TO  77 
IF  I IAZ2.LT. 180)  G0T07I 
GO  TO  69 

68  IFI IAZ2.LT .KAZ.OR.1AZ2.GT.LASTAZ) GO  TO  77 
69  I 1 = 1 AZ2-KI I 
GO  TO  72 

71  1 1 = 1 AZ2  + KK 1 1 

72  CONTINUE 

IFI INOEX.EQ. I 1G0T0166 

IFI IEL2I 1 I .GT. NINE. OR. IEL2I 2) . GT .NI NE .OR . I EL2I 3 ) . GT .N INE » GOTO 166 
CALL  ATOI I IELI .3.0. IEL.NTEST) 

IFI NTEST.NE .01 GOT  09999 
T !LT  = I EL/1  0.0 
166  CONTINUE 
I NOEX=  2 

IFITILT-PTILT.GE. 0.4 IG0T078 
IF  ITILT-PTILT.LT.-0.5IGO  TO  99 
87  CONTINUE 

C • 

IFI IRI.EQ. IONE1GO  TO  70S 
C 

OO  700  1=1 .180 

700  ICLASSI I+20)=IHWOROI 1.430) 

GO  TO  706 


'7T 
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705  DO  701  1 *2 *1  BO  .2 

701  ICLASSt I*20I=lHwOROf I/2F215I 

706  DO  II  I=IAODR.IOO.IAOOR 

k=i+krange 
ZWt II. I 1=0.0 
ARNG=FLOATtKI 

IF  I 1 CLASS IK I.LE.34) GOTO 1 1 

ZWt  I I . I > = 1 0.*«<-0. l*CALI B(  ICLASSt K I l*2.*ALOGl 0< ARNG ) ♦ HAVECI 

IFITILT.LT.0.5»G0TO321 

GOTO  It 

321  ZINTRp=Z*t  II  ,1 ) 

IFf ZlNTRP.LE.l .OlGOTOl 1 

PVIZl t I I • I I=PVIZ| ( It . I |»ZINTRP/IDtVlt l> 

PZI II. I ) = Z INT3P 
II  CONTINUE 

IFITILT.LE.O. 5 1G0T077 
CALL  I NTRP2 ( TILT .III 
77  CONTINUE 

IFt IRI .EO. IZROIGO  TO  IS 

IF  THE  RANGE  INTERVAL  IS  2 KM.  THE  THIRD  AND  FOURTH  R40IALS 
ARE  PROCESSED  AS  ABOVE.  IF  IN  I KM  MODE.  ANOTHER  DATA 
RECORD  IS  read. 

IFIKAZ.LT. LASTAZ IGO  TO  368 

IF| I AZ3.LT. KAZ. AND. I AZ3.GT.LASTAZIGO  TO  377 
IFI IAZ3.LT.160IG0  TO  371 
GO  TO  369 

368  IFI IAZ3.LT. KAZ. OR. I AZ3.GT .LAST AZ IGO  TO  377 

369  I 1 = 1 AZ  3-KI I 
GO  TO  372 

371  11  = 1 AZ  3AKK 1 1 

372  CONTINUE 
IFt IN0EX.E0.2IG0  TO  366 

IFt I El 2( 1 1 .GT. NINE. OR. IEL2I 2 I . GT. N I NE . OR. IEL2I 3 I .GT.N INE I GO  TO  366 
CALL  ATOI t IELI .3.0. IEL.NTESTI 
IFtNTEST.NE.OlGO  TO  9999 
TILT=IEL/1 0. 0 
366  CONTINUE 
I NOEX=  3 

IFtTILT-PTlLT.GE.0.4!G0  TO  78 
IFtTlLT-PTILT.LT.-0.51G0  TO  99 
68  CONTINUE 

DO  300  1*2 • 1 80*2 
600  ICLASSt lF20)=IHWORDt 1/2F430I 
DO  410  1=2.100.2 
K=|FKR*NGE 
ZWt 11.11=0.0 
ARNG=FLOAT  tK) 

IFf I CLASS! K) .LE.34IGO  TO  410 

ZVt 1 1 . I 1 = 1 0.**f-0. !*CALIBt ICLASSt K) • *2.4 ALOG 1 Ot ARNG > F VA VECJ 
IFtTiLT.LT.O.SIGO  TO  331 
GO  TO  410 

331  ZINTRP=ZWt II  .1  1 

IFf ZINTRP.LE. I.OIGO  TO  410 
PVI Z 1 1 1 1 , I >=PVlZ 1 ( 1 1 . I l+ZINTRP/IDIVlt I I 
PZf I I . I l = Z INTRP 
410  CONTINUE 

IFITILT.LE.O. 51G0  TO  377 
CALL  I NTRP2t  TILT. Ill 


uuu  u uuouu  uuuu 
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377  CONTINUE 

PROCESSING  OF  THE  FOURTH  RADIAL  BEGINS  HERE 
IF(KAZ.LT.LASTAZ)GO  TO  463 

IF!1AZ4.LT.KAZ.AN0.IAZ4.GT.LASTAZ)G0  TO  IS 
IF f I AZ4  «LT • 1 80)00  TO  471 
GO  TO  469 

463  IF( I AZ4.LT.KAZ.0R. I AZ4.GT.LASTAZIGO  TO  IS 
469  II*IAZ4-KII 
GO  TO  472 

471  II=IAZ4.KKII 

472  CONTINUE 

IFI INDEX. E0.3IG0  TO  466 

IFl IEL2I I I .GT. NINE. OR. IEL2I2) . GT.NINE .OR * I EL 21 3 1 . GT .NINE (GO  TO  466 
CALL  ATOI ( IELI .3.0. IEL.NTEST) 

IFINTEST.NE. OIGO  TO  9999 
TILT=I EL/t 0.0 
466  CONTINUE 
In0Ex=4 

IFITILT-PTILT.GE.0.4IGO  To  78 
IF(TILT-PTILT.LT.-0.5)GO  TO  99 
89  CONTINUE 

DO  900  1=2.180.2 

900  ICLASSI I+20)=1HW0R01 1/2*6451 
DO  411  1=2.100.2 
K=I«KRANGE 

zwl it. i 1=0.0 

ARNG=FLOAT(K| 

IF( ICLASSI K> .LE. 34 IGO  TO  411 

ZWI 11.11  = 10. **(-0. 1*CALIBI I CLASS! Kl ) 42. 4AL0G I 0 ( ARNG > ♦ WAVEC t 
IFItJlT.LT. 0.5)00  TO  431 
GO  TO  411 

431  ZINTRP=Zwt II .1  I 

IFI ZINTRP.LE. 1 .0 »GO  TO  4t! 

PVIZII II .1  )=PV IZ 1 1 II . I >»ZlNTRP/IOIVll I) 

PZf 1 1 • I l=Z  JNTRP 
411  CONTINUE 

IFITILT.LE.O.SIGO  TO  477 
CALL  I NTRP2I TILT. I!) 

477  CONTINUE 

GO  TO  IS 

THIS  PROCEDURE  CONTINUES  UNTIL  ALL  OF  THE  DATA  ARE  READ  FOR  A 
TILT  SEQUENCE. 


78  CONTINUE 

I COUNT  = l COUNT  4 l 
JFIICOUNT.GT.il  G0T083 

the  following  cooe  oevelops  and  plots  the  zero  degree 

REFLECTIVITY  MAP.  USING  SUBROUTINE  PLOTZ. 

DO  81  1 1 = 1 . I 00. I ADO 
DO  81  1=1 AOOR . 100.IADDR 
ZPI It. I+KRANGE>=0.0 
81  ZPIII.I»KRANGEI=ZWl II. I) 

CALL  QD2I 1 LEFT • JO OWN , KI I . KK III 

KM=I 

LVL=0 


123 


I \\ 

I n 


c 

c the  Zlx.vl  array  is  converted  to  dbz  and  plott.ed. 

c 

OO  82  1*1. SI 
DO  82  J*l.51 

IF  (Z(I.J).LE.l. 0)2(1. J)=l.O 
ZII  . Jl =I0.4AL0G1 01 Z( I. J) I 

82  CONTINUE 

CALL  PLOTZIIDATE.NENDT.LVL. ILEFT. JDOWN. MM I 

c 

C the  ARRAY  Z IS  WRITTEN  TO  TAPE  FILE  IRIT  TO  DE  USED  LATER  WITH 

C A CONREC  SUBROUTINE  WHICH  CONTOURS  THE  DATA  ANO  DRAWS  MAPS  ON 

C THE  CALCOMP  PLOTTER. 

c 

WRITE! IRIT .6006) Z 
6006  FORMAT ( 2 OOF 5. 1 , 69F5. I I 

83  CONTINUE 

DO  84  1=1 . 100 

84  HT 1(1) = HT2  < I > 

PTILT=T ILT 

C 

CO  T0I86.87. 88. 891. INDEX 
C 

C GOES  BACK  TO  WHERE  IT  LEFT  OFF  PROCESSING  THE  RADI ALS 

C 

C WHEN  TILT  BECOMES  MORE  THAN  .5  DEG  LESS  THAN  PTILT,  THE 

C PROGRAM  ASSUMES  THAT  THE  SCAN  SEGUENCE  IS  OVER  ANO  TERMINATES 

C PROCESSING.  IT  THEN  BEGINS  TO  GENERATE  THE  OUTPUT  MAPS. 

C 

99  CONTINUE 
C 

C THE  LOWER  LAYER  PVSZ  MAP  <0-15  KFT I IS  GENERATED  FIRST 

C AND  WRITTEN  ON  THE  OUTPUT  TAPE  (IRIT) 

• C 

MM=2 

LVL-MM45000 
DO  TT7  11*1.1 00. I ADD 
DO  777  1=1 ADDR.100. IAODr 
ZPIII  . I 4 KRANGE)  *0.0 
C 

C THE  dummy  ARRAY  ZP  IS  LOADEO  WITH  THE  APPROPRIATE  PVSZ  ARRAY 

C ' AND  SENT  TO  SUBROUTINE  002. 

C 

ZPIII  .14  KRANGE)  = 3 * PVIZlfll.I) 

777  CONTINUE 

CALL  0D2! ILEFT. JOOWN.KU ,KK1I ) 

Do  150  1=1,51 
DO  150  J=1.51 

!F(Z< t .JI.LE. 1.01 ZIt.J 1=1.0 
Z< I.J)=l0*ALOG!0!Z!l.J)) 

150  CONTINUE 

CALL  PLOTZIIOATE.NENDT.LVL. ILEFT. JDOWN. MM  I 
WRITE! IRIT . 6006) Z 
C 

C NEXT  IS  THE  MIDDLE  LAYER  PVSZ  MAP  (15-35  KFT) 

c 

MM=3 

LVL-MMF5000 
DO  778  11=1.100.1 AOD 
DO  778  1= IAODR. 100. IAODR 
ZPIII  .14  KRANGE)  =0,0 
ZP(II  ,14  KRANGE)  * 4 * PV!Z2<II.I> 
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778  CONTINUE 

CALL  002! (LEFT. JOOWN. KII.KKII I 
00  160  1*1.5! 

DO  160  Jsl.SI 

IF(Z( t.JI.LE.t.O)Z(|,JI=I.O 
Z( I. J) =10*ALOG10( 2! I.JU 
160  CONTINUE 

CALL  PLOTZI IOATE.NENOT.LVL. (LEFT. JOOWN. MM 1 
WRITE! IRIT .60061 Z 
C 

C NOW  THE  UPPER  LAVER  PVSZ  MAP  (35-50  KFT)  IS  CENERATED 

C 

MM*4 

LVL*MM*5000 

00  779  11=1.100. IAOO 

00779  1*1 AOOR.IOO. IAOOR 

ZP(  I I .14  KRANG6)  >0.0 

ZP(II  . I ♦ KRAnGE)  = 39  PVt  Z3( 1 1 • 1 1 

779  CONTINUE 

CALL  002(1 LEFT .JDOwN.KII , KK II) 

DO  170  1*1.51 
00  170  J=1.51 

IF(Z( I . J I .LE« |,0IZ(t.J)=l.0 
Z( I » J ) = I 0 * AL  OG 1 0 ( Z! I . J) > 

170  CONTINUE 

CALL  PLOTZdOATE.NENDT.LVL.  ILEFT  . JOOWN.  MM  | 

WRITE! IRIT .6006 I Z 
C 

C VIL  IS  CALCULATED  IN  A SIMILAR  MANNER.  vlL  IS  MEASURED 

C IN  KG/M**2 

C 

MH=5 

LVL=MM*5000 
00  7 00  11  = 1. 100, IAOO 
OO  700  I=IAOOR. 100. IAOOR 
ZP( I I • I TKR  ANGE 1 = 0.0 
700  ZP! I I . I*KRANGE)=VIZ!  I I .1 1 

CALL  002! ILEFT. JDOWN.KI I .KKI I I 
00  100  1=1.51 
DO  100  J=1.5l 
Z!I.JI=Z( I .Jl/1000.0 
100  CONTINUE 

CALL  PLOTZIIOATE.NENOT.LVL. ILEFT . JOOWN, MM ) 

WRITE! IRIT .6006 1 Z 
ENOFILE  IRIT 
C 

C THIS  PUTS  AN  ENO-OF-FILE  MARK  ON  THE  TAPE.  SO  THESE  5 MAPS 

C COMPRISE  ONE  FILE 

C 

WRITE! 6.2000) IRIT 

2000  FORMAT! • 0* , • ENO  OF  FILE  WRITTEN  ON  TAPE  ».I21 
IF  ! I OS.EO. I ONE  1 WR I TE! 6 • 72 l 7 1 

7217  format  do*.*  oegree  select  * 2*  1 

IF! IRI.EO. IONEI  WRITEI6.7317I 
7317  FORM AT I* 0’.’  RANGE  SELECT  = 2 KM. • ) 

C 

C THIS  INOICATES  WHETHER  EITHER  2 DEG  OR  2 KM  MODE  DATA  WAS 

C PROCESSED.  DEFAULT  VALUES  ARE  1 OEG  AND  1 KM, 

C 

00  781  11=1.100 
00  781  1=1 ,100 
PVIZKII, 11  = 0.0 


I’sy.i T 


PVIZ2< II, 1 )=0.0 
PVIZ3( 1 1 . I > = 0.0 
T81  VIZ(  II  .0  = 0.0 
00  782  1 = 1 .100 
MT2< I >=0.0 
782  HTI( !)=0.0 
C 

C REINITIALIZES  THE  ARRAY  VALUES  TO  ZERO. 

C THE  PROGRAM  IS  NOW  READY  TO  READ  ANOTHER  INITIAL  DATA  CARO 

C AND  BEGIN  PROCESSING  ANOTHER  SCAN  SEOUENCE. 

C 

GO  TO  I 

999  WR1 TEt 6. 300) 

300  FORMAT! *0* «• ********  READ  ERROR  ENCOUNTERED  ***»*»*»»******•> 
GO  TO  99999 

9999  WRI T E ( 6 • 30 1 ) 

301  FORMAT! •O’ .* *******  ATOI  NTEST  ERROR  SS»$**SSSSS*S»****SS*S**« ) 
99999  CONTINUE 

STOP 

ENO 

C 

C 

c 

c 

BLOCK  DATA 

COMMON/BLOCK A/ZV!  1 00 . I 00 >. P VI Z 1 1 1 00 . I 00  I , PVl Z2<  10  0. 100  1 .PVIZ3< 100. 
1 1001 .VIZ! I 00. I 00> .KRANGE .A! 4,4  I .IDIV1!I50>.IOIV2I150>.IDIV3!150>. 
2ZP1 100.200) .PZ I 100.1001.2(51.51 ) .CORECTI  I 50) . ICOUNT.hU ( 100) . 

3HT2!  100) . I ADO. IADDR 

DATA  A!  I .2 ) . A!  2.3) . A! 3.4>/3»-l. 0/ . A! I ,3) . At  2 . 4 > /2*-2. 0/. A! 1 ,4>/-3. 
70/.A!2.1).A! 3.2) .At4.J)/3*  1.0/.A(3.1),A(4:2)/2*2.0/.A< 4, 11/3.0/ 
DATA  V IZ.PVl Zl .PVIZ2.PVI Z3/40  0 0 0*0.0/ . 10  I VI . I D I V2 . ID  I V3/450  4 0/ 

DATA  ZP/2000040.  0/. Z /260 I *0 . O/.CORECT/I 50*0. 0/. PZ/ 1 0000*0. 0/  . 

END 

C 
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RSQ  » R**2 

CORECTI  I 1 = 10.375  * RSQ)/  6371.315 

C 

C 6371.315  IS  THE  RADIUS  OF  THE  EARTH  IN  KM.  IF  SOME  OTHER  UNIT 

c of  oistance  is  useo  in  the  digitazation  process,  then  the 

c RADIUS  OF  THE  EARTH  IN  THOSE  UNITS  MUST  BE  USEO  INSTEAD  OF  THIS 

C NUMBER. 

C 

HTn=4.572»C0RECT< I I 
HT22=10.66a-CORECTI I 1 
HT33=!5.240-C0RECT(! ) 

ELEV 1 = ATAN2(HT 1 1 .R 1*57.29578 
ELEV2=ATAN2IHT22.R) *57.29578 
ELEV3=ATAN2I HT33.RI*57. 29578 
C 

C IF  THE  ELEVATION  ANGLE  COMPUTED  IS  GREATER  THAN  THE  LAST  TILT 

C ANGLE  INCREMENTED  BY  1 DEGREE.  THEN  THE  01 VISOR  MUST  BE 

C CALCULATED  DIFFERENTLY. 

C 

I F( ELE VI *GT • SU0TLT1 G0T01  O 
I0IV1I  I ).  = I FI  XI  ELEV  1 1 
It  CONTINUE 

IF! ELE V2.GT • SUBTLT IGOT020 
IOIV2III  = I F I XI  ELEV2 1 - lOIVltll 

12  CONTINUE 

IFIELEV3.GT. SUBTLT 1 GO TO 30 

IDIV3II)  * IFIXIELEV3)  - IIDIVltll  * IDIV2II1I 
GOTO  1 • 

10  IDIVlt  I )•  = IFIXIIELEV1  - SUBTLT  1/2.0)  ♦ NTLT 

GOTO  1 1 

20  IDIV2III  = I F I XI 1 ELE V2  - SUBTLTI/2.0)  ♦ NTLT  - IDIVltl) 

GOTO  12 

JO  IOI V3I I ) = IFIXt  (ELEV3-SUQTLT1/2.0DNTLT  -I  1 D1  Vl  I I > » I DI  V2I  I I) 

C 

IFI1DIVII 1 I.EO.O) IDIVlt I »=! 

C 

1 CONTINUE 

RETURN 
ENO 
C 

c 

c 

c 

subroutine  intrp2<  tilt. id 

DIMENSION  RM I A 1 

COMMON/0LOCK  A/ ZW ( 1 00 . 1 00  » . PVI Z 1 1 1 00 . 1 00 ) . PVI Z2I  100 . 1 001 . PV I Z3I 1 00 . 
1 1001 .VIZI  I 00. 1 00 J .KRANGE.AI4.4 1 .101 VI ( 1501, 101 V2< 1 50 ) . I DI V3 1 I 50 > . 
2ZPI I00.200).PZ(100.100>.ZISI.S1 > .CORECTI  1 50 ) . I COUNT . HT 1 1 100) • 

3HT2I 100) .1  ADD. IADDR 

c 

C THIS  SUBROUTINE  USES  THE  CODED  AZIMUTH  III)  AND  THE  ELEVATION 

C ANGLE  TO  PERFORM  THE  LAGRANGIAN  INTERPOLATION  ALONG  THE  RADIAL. 

C 

c 

c 

C DEFINITIONS 

C 

c PHI-  THE  TILT  ANGLE  IN  RAO IANS 

C R - RANGE  OF  THE  SUMMATION  COLUMN  AT  THE  ZERO  DEGREE  TILT  ANGLE 

C H - HEIGHT  OF  THE  DATA  POINT  IN  KM. 

C HT2II)  - THE  HEIGHT  IN  METERS 

C IHT  - HEIGHT  OF  THE  DATA  POINT  IN  KFT 
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*PRN  - RADIAL  DISTANCE  TO  THE  INTERPOLATION  DATA  POINT 
IFI IAOOR.EO. I >GO  TO  77 

INITIALIZE  OFF»OI AGONAL  ELEMENTS  OF  A FOR  2 KM  MODE. 

A( 1. 2)3-2. 0 
Ait. 3)^.4. 0 
A< 1. 4)3.6. O 
A< 2. I > = 1.0 
A<2. 3)3.1. o 
AC  2. A) 3.2. 0 
AC  3. 1 ) =2.0 
A<3. 2)31.0 
AC3.4) 3.) • 0 
A<4. 1)33.0 
A(4. 2)32.0 
A(4. 3)31.0 

77  PHI=TILT*0. 01 74532 

00  I I 3 | . 1 00 . I AOOR 

THE  INDICES  FOR  THE  FOUR  DATA  POINTS  USED  IN  THE  LAGRANGI AN 
INTERPOLATION  ARE  I.  1*1.  1*2.  AND  I-I  FOR  THE  I KM  MODE  ANO 
I.  142.  144.  AND  1*2  FOR  THE  2 KM  MODE. 

IP 13) 4 I AODR 
IP23I424I ADOR 

1 M1*I**I  ADDR 

IFC IMI .LE.O.OR.IP2.GE.IOO)GOTOI 

ZINTRPxO.O 

R3  FLOATi I4KRANGE) 

M=  R4TANI  PHI >4CORECTII ) 

HT2I I )=H*I 000.0 

IFC ICOUNT.LT. I >HT1 1 1 )3HT2( I I 

DH=MT2I I l-HT I < I I 

IFCZWC  1 1 .1  ).  170. 0.0.  ANO.  ZWC  I I.  IP  I ).EQ.  O.OIGOTOl 
IHTsIFIXf H43.2BO044I.O) 

RPRM=SQRTCR**24M**2> 

OELR-RPRM-R 

IFC I ADOR • EQ. 2 > GO  TO  20 

IFCDELR.LT. 1 .0)XCNT30.0 

IN  THE  I km  MODE.  IF  DELR  IS  GREATER  THAN  I BUT  LESS  THAN  2. 
THEN  THE  INTERPOLATION  POINT  LIES  BETWEEN  141  AND  142 

IF! DELR .GE .1 • 0. ANO.DELR.LT*  2.0) XCNT  =1.0 
IFC  DELR. GE .2 • 0 . ANO.DELR.LT .3.0) XCNT=2. 0 
IF  C DELR.GE .3. 0 . ANO.DELR.LT . 4. 0 1 XCNT=3 . 0 
IF C DELR .GE. 4.0. AND. DELR. LT. 5. 0 I XCNT=4.0 
CO  TO  22 

IN  THE  2 KM  MOOE.  THE  INTERVALS  ARE  DIFFERENT. 

20  IFCDELR.LT. 2. 0>XCNT=0.0 

IFC DELR.GE. 2.0. ANO.DELR.LT. 4. 0>XCNT=2.0 
IFC OELR.GE. 4.0. ANO.DELR.LT. 6.0) XCNT=4.0 
IFC DELR. GE.b.O. AND.OELR.LT. 8.01 XCNT=6.0 
I F< DELR.GE. fl. 0. ANO.DELR.LT. I 0.01 XCNT= 8.0 
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OEPENOS  ON  THE  VALUE  OF  OELR, 

22  IP l = IP  1 ♦XCNT 
IP2=IP2*XCNT 
IMlsIMI AXCNT 

IF< IMI. LE. 0. OR .IP2.GE <100) GOTOl 

IF  THE  VALUE  OF  THE  REFLECTIVITY  ON  EITHER  SlOE  OF  THE  INTER- 
POLATION POINT  ARE  EQUAL « THEN  NO  INTERPOLATION  IS  NECESSARY. 

IF<ZW< I I .IMI *1 AOOR) <EQ.ZW< I I. IPI I 160  TO  16 

DELR=DELR— XCNT 

IFC DELR<LE<0<01 ) GOTOl 7 

THE  OIAGONAL  ELEMENTS  OF  THE  MATRIX  USED  FOR  THE  LAGRANGIAN 
INTERPOLATION  ARE  COMPUTED<  THE  OFF— O I AGONAL  ELEMENTS  ARE 
ALWAYS  THE  SAME  FOR  A PARTICULAR  RANGE  INTERVAL  < ANO  ARE 
ALWAYS  SET  AT  THE  BEGINNING  OF  THE  PROGRAM< 

ADDR=FLO AT  1 1 AOOR ) 

A(  1 < I ) =1  < 0*AOOR«-OELR 
A<  2 • 2 ) — DEL  R 
Af  3*  3) =— 1 < 0* AOOR* DEL R 
A ( 4 , 4 ) =— 2 • 0*  AOOR  4DELR 
00  3 K— 1,4 
PR00=1 <0 
00  4 J=I .4 
PRO0=PRO04A<  K.  J» 

RMOO-PROO 
SUM-0 <0 
K 1 = 0 

00  S K — I Ml .1 P2  < I AOOR 
Kl=Kl+  t 

ZPRM  = ZW(  I I ,K  l/RMl KI I 

SUM=SUM<ZPRM 

CONTINUE 

OIAG=I <0 

DO  6 K=1 ,4 

DI AG=0 I AG4 A<  K,K) 

CONTINUE 

. . ZInTrp=OIAG*SUM 
GOTO  10 

6 ZINTRP=ZWI II .IPI I 
GOTO 10 

7 ZINTRP=ZW<  I I . I Ml ♦ I AOOR) 

0 CONTINUE 

IF<ZINTRP<LE<  I. 0>ZINTRP=1.0 

WHEN  ICOUNT  IS  I.  WE  ARE  IN  THE  BASE  OR  ZERO  LEVEL  ANO  CANNOT 
INTEGRATE. 

IF<IC0UNT<LT<1 IGOTOIOO 

' THE  integration  of  liouio  water  content  is  accomplished. 

VIZI  II .1 )«VIZ< 1 1. I 14<  3.44E— 34  CpZf II, I 144.571 »ZlNTRP44. 571 1/2. 14DH 
100  CONTINUE 

PZ< II . I 1 =Z INTRP 
IF<ZINTRP<LE<  I <01 GOTOl 

IHT  IS  TESTED  TO  OETERMI NE  INTO  WHICH  LAVER  THE  INTERPOLATED 
REFLECTIVITY  <ZINTRPI  GOES. 
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c 


c 

55 


C 

56 


1 


C 

C 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


IF( I HT .GT. IS.ANO. IHT .LE.35 JGOTOS5 
I F< I HT . CT . 35 1 G0T056 

PVIZK  I I.  I )=PVIZI<  II.  I )FZI NTRP/FLOATI  IDIVKII  ) 
COTOI 

CONTINUE 

JFUDIV2II  >.£0.0)  GOTOI 

PVIZ2I l I . I I— PV I Z2( 1 1 .1 I+ZINTRP/FLOATI 1DIV2II 11 
GOTOI 

CONTINUE 

IF ( IHT.GT.50IGOTOI 
IFCI01V31I  I .LE.OIGOTOI 

PV  l Z3<  tl.t  l = PVIZ3(  II  .1  14-ZI  NTR  P/FLO  AT  IIOIV3IIII 

CONTINUE 

RETURN 

END 


SUBROUTINE  OD2  ( ILEFT. JOOwN.KI I .KKI I I 

COMMON/ BLOCK A/ ZW( IOO.tOO>,PVIZl<tOO.lOO).PVIZ2<IOO.IOOI.PVIZ3(IOO< 
I 1001 .VI Zt I 00. I 001 .KRANGE, Al  4.AI.IOIVII150I . I 01 V2<  1501 .101 V3C 1501 . 
2ZPI  100.200 I.PZI I00.100I.ZI51.5I  ) . CORECT ( 1 50 J . I COUNT, HT It  1001 . 

3HT2I 1001 • I AOO  * lAOOR 
INTEGER  AZONE.AZTWO. azthre. azmth 

This  SUOROUTINE  steps  through  the  rectangular  grid  coordinates 
and  computes  a value  for  each  using  an  interpolation  scheme  which 
selects  data  from  the  9 nearest  cylindrical  data  points  contained 

IN  THE  ZP  ARRAY.  the  OUTPUT  ARRAY,  Z.  IS  RETURNED  TO  THE  MAIN 
PROGRAM  FOR  PLOTTING  OR  FURTHER  USE.  SEE  VOGEL'S  I 19731  THESIS 
FOR  EXPLANATION  OF  THE  INTERPOLATION  SCHEME. 


DEFINITIONS 


Float i - distance  west  <♦>  or  east  i-i  of  tamu  in  km  to  the  x 

GRIO  COORDINATE  It) 

FLOATJ  - DISTANCE  NORTH  <♦)  OR  SOUTH  (-1  OF  TAMU  IN  KM  TO  THE  T 
GRID  COORDINATE  IJI 

R - THE  DIRECT  DISTANCE  FROM  TAMU  TO  THE  GRID  POINT  IN  KM 

THETA  - ANGLE.  IN  OEGREES  FROM  TRUE  NORTH,  OF  THE  GRIO  COORDINATE 

AZMTH  - THE  NEAREST  OATA  AZIMUTH  TO  THETA 

AZONE  - THE  NEAREST  AZIMUTH  TO  THETA  IN  CODED  FORM 

REAL I - The  NEAREST  AZ I MUTM  IN  OEGREES 

X - PERCENTAGE  OF  DISTANCE  FROM  NEAREST  ANGLE  TO  ANGLE  OF 
THE  GRIO  COORDINATE 

l«!  - THE  NEAREST  RANGE  GATE  TO  THE  GRIO  COORDINATE 
REAL J - THE  RANGE  OF  THE  RANGE  GATE  IN  KM 

V - PERCENTAGE  OF  DISTANCE  FROM  NEAREST  RANGE  GATE  TO  GRID  POINT 
AAA  THROUGH  III  - THE  9 DATA  POINTS  TO  BE  USED  IN  THE  QUADRATIC 
INTERPOLATION 

8 THROUGH  F AND  FONE  - INTERPOLATED  VALUES  OF  REFLECTIVITY 
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c 


DO  260  1*1. 51 
FLOATI=FLOATI2*I-ILEFTI 
00  260  J*l .51 
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Z tI,J)=0.0 
FLOAT J*FLO AT < 2* J- JOOWN  > 

R*S0RTtFL0ATI**24FL0ATJ**2l 
IFtR. LT. 2.0. OR. R.GT. 200.0)  GO  TO  260 
THETA=t ATAN2tFLOATI . FLOAT J I 1*57.2958 
IFtTHETA.LT. 0.0)  THET A=THETA* 360 
AZMTH= IF  IX ITHETAF0.5 I 

DIFFERENT  increments  MUST  qe  specified  for  either 
THE  2 DEG  OR  2 KM  MODE. 

IF«I ADD.EQ. 21 AZMTH= IF IXt THETA 1/2*2* 1 

THIS  INSURES  THAT  AZMTH  IS  ODD  IN  THE  2 DEGREE  MOOE 

IFtKI I .GT. 180. AND.AZMTH.LT. iao>GO  TO  213 
AZONE=AZMTH~KI I 
G0T021 4 

213  AZONE=AZMTH«-KKIl 

214  AZTWQ-AZONE+ I ADO 
AZTHRE=AZONE-IADD 
REAL I=FLOAT< AZMTH) 

X=THET A-REAL I 
IRI=»IF  IXtR+O.S) 

I Ft  I ADDR.EQ. 2) lRl=t IF  I X( R*. 999 )♦ I ) /2* 2 

THIS  INSURES  THAT  I R1  IS  EVEN  IN  THE  2 KM  MODE 

IR2=IRI+IA0DR 
JR3=IRl-IA00R 
REALJ^FLOATt IR1 ) 

V=R«REALJ 
C 

AAA=  ZP f AZONE . IR1 ) 

BOB— ZP  t AZONE . I R2 ) 

CCC*ZPt  AZONE . IR3) 

OOD»ZPt AZTWO. IRI ) 

EEEsZPt AZTW0.IR2) 

FFF*ZPtAZT*(0.IR3l 
• GGG=ZP  t AZTHRE. IRI) 

MHH=ZPt AZTHRE. IR2) 

III«ZPtAZTHRE.IR3) 

C 

IFC  AAA.NE .0. 0. OR .BBB .NE .0.0 .OR.CCC.NE .0.0 .OR .DUO.NE. 0. 0. OR.EEE.NE. 
CO.O.OR.FFF.NE.O.O.OR.GGG.NE.O.O.OR.HHH.NE.O.O.OR. I 1 1 .NE. 0.0)  GO 
CTO  236 
GO  TO  260 
C 

236  B*f D00-GGG)/2.0 
C»tBUS«CCC>/2.0 
0«DD0-AAA-B  ■ 

E*CCC-AAA*C 
IFtX.GE.O.O)  GO  TO  24B 

IFt V.GE.O*  0)  GO  TO  24S  .. 

F*I I I-AAA*B*C-D-E 

FONEMIt 

GO  TO  2S4 

c 

249  F* AAA-S»C»0»E«HHM 

FONEsHHH 
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n 


c 

lit 

121 


C 

I 12 
122 
103 


CO  TO  254 

IF1Y.GE.0.0)  GO  TO  252 
FxAAAFB-CFOfE-FFF 
FONE=FFF 
CO  TO  254 

f=eee-aaa-b-c-d-e 

F0NE=EEE 

ZZ  = AAAKO*XIM  C*Y)F|0*X4*2)f!E*Y**2)f!F*X«Y) 

IF  THE  INTERPOLATED  VALUE  OF  THE  REFLECTIVITY  IS  NEGATIVE*  THEN 
DO  A LINEAR  INTERPOLATION  USING  THE  4 DATA  POINTS  SURROUNDING 
THE  GRIO  POINT. 

IFtZZ.LT. 0.0 IGOT0500 

ZZxAMAXlIZZ. AH IN 11  A A A . DDO . GGG . 08B. CCC .FONE ) > 

IF( ZZ.LT.O .0 ) GO  TO  260 

Z !I.J)xZZ 

COTO260 

CONTINUE 

IRl  x IFIXl  R | 

IF! 1AODR.EO. 2) IRlxl IF  I X! Rf. 999 )F I ) /24 2 

I R2*I  R I FI AOOR 

DR  X R - FLOAT!  IRl  » 

DELT  x THETA  - FLOAT!  AZMTH  I 
IF10ELTI 100. I 01  . I 02 
OTx(I .OFABS! 1 .OFOELTI » 

ifikii.gt. tao. and. azmth. lt. iao>Go  to  no 

AZTWOxAZMTH-KI  I 
GOTO  t 2 0 

AZTWO* AZMTHFKK 1 1 
AZONE= AZTwO-I AOD 
COTO  103 

CONTINUE 

IF1KI I *GT .160. AND. AZMTH. LT. 1 00  I GO  TO  III 
A20NE=AZMTH-KI | 

GOTO  121 

AZONEeAXMTHFKKII 
AA  x ZPI AZONE. IR2) 

BB*ZP ! AZONE. IRl) 

COTO 104 

DTxDELT 

IFIAZMTH.LT. 180)  COTO 112 
AZONExAZMTH-KI I 
COTO  122 

AZONExAZMTHFKKII 
AZ TWO3 AZONE F I ADO 

AA  3 ZPI AZONE. IR2)  F ! ZPI AZTW 0. I R2 ) - ZP I AZONE. IR2) I • OT 
BB  x ZP! AZONE. IRl > F I ZP! AZT VO . I R 1 > - ZP1AZ0NE. IRl ) )•  DT 
Z II. J)  x BB  F ! AA  - BB)  * OR 


260  CONTINUE 
RETURN 
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END 


SUBROUTINE  PLOTZ  ( NO. NT.L VL . IL . JO, MM) 

DIMENSION  6UF( 1261 .OICI I0).NT(2) 

COMMON/BLOCK A/ZWI  100. I 00  > .P VI Z I I 1 00 . I 00 ) . PV I Z2(  1 00 . I 0 0 I . P V I Z3 ( I 00  . 

I tool .VIZI  1 00.100)  .KRANGE .A< 4,41.(01 VI ( 150) ,XOI V2< ISO) .101 V3I 150)# 
2ZPI 100.200) .PZ I 1 00. I00)*Z(51 ,51 ) .CORECTI  ISO) . 1COUNT.HTI ( 1 00) . 

3HT2I  100) • I ADO. IAOOR 
OATA  BLNK/IH  /.PLUS/IH4/ 

DATA  DIO/ 1 HO .IHt.IH2.lH3.lH4, 1HS. )H6. 1H7. IMS . 1H9/ 

This  subroutine  takes  a so  x so  rectangular  grid  array  and  plots 

FIRST  THE  LEFT  HALF  AND  THEN  THE  RIGHT  HALF  OF  The  ARRAY.  THE 
RESULTANT  PLOT  IS  25  INCHES  SQUARE  WHEN  BOTH  HALVES  ARE  PIECED 
TOGETHER.  AT  A 2 KM  GRIO  INTERVAL.  THE  MAP  COVERS  A 100  KM  SQUARE. 

starting  with  the  upper  left  hano  corner  of  the  rectangular  grid. 

DBZ  VALUES  FROM  THE  Z ARRAY  ARE  LOAOEO  A LINE  AT  A TIME  INTO  THE 
BUF  ARRAY.  FIRST.  THE  126  CHARACTERS  IN  BUF  ARE  SET  EQUAL  TO  A 
BLANK.  THEN  EVERY  25  CHARACTERS  ARE  SET  EQUAL  TO  A PLUS  (4).  THESE 
GRID  MARKS  REPRESENT  A DISTANCE  OF  10  KM. 


DEFINITIONS 

ARRAYS 

BUF  - STORAGE  LOCATION  OF  DATA  TO  BE  PRINTED  ON  EACH  LINE. 

126  CHARACTERS  LONG 

DIG  - CONTAINS  CHARACTER  REPRESENTATION  FOR  INTEGER  VALUES  0-9 
OTHERS 

BLNK  - A HOLLERITH  BLANK  CHARACTER 
PLUS  - A HOLLERITH  PLUS  SIGN  (♦) 

NN  - COUNTER  WHICH  STEPS  ACROSS  THE  PAGE  IN  INCREMENTS  OF  FIVE 
PRINT  LOCATIONS  REPRESENTING  2 KM  INTERVALS 
NUM  - INTEGER  VALUE  OF  DBZ  FOR  GIVEN  GRID  POINT 

NO 1 — COOEO  VALUE  REPRESENTING  100* S OIGIT  OF  DBZ  VALUE  I SHOULD 
BE  ZERO) 

N02  - COOED  VALUE  REPRESENTING  10* S OIGIT  OF  OBZ  VALUE 

N03  - CODEO  VALUE  REPRESENTING  UNITS  OIGIT  OF  DBZ  VALUE 

IW  - DISTANCE  IN  KM  EAST  or  WEST  TO  LOWER  LEFT  CORNER  OF  GRIO 

NS  -DISTANCE  IN  KM  NORTH  OR  SOUTH  TO  LOWER  LEFT  CORNER  OF  GRIO 

INN  - COUNTER  USEO  TO  COMPUTE  NN  FOR  RIGHT  HALF  OF  MAP 

WR1 TEI6.400) 

DO  330  JJ^l.51 
J=52-JJ 


DO  310  1 = 1 .126 
BUF< I)=OLNK 
HO  CONTINUE 

IF( MOO ( ( J— 1 ) .5 ) • NE.O ) GO  TO  315 

DO  314  1 = 1 .126.25 
BUFI IlsPLUS 
314  CONTINUE 

OBZ  DATA  FROM  THE  Z ARRAY  ARE  ROUNDEO  TO  THE  NEAREST  INTEGER 
ANO  LOAOEO  INTO  THE  PROPER  LOCATION  IN  THE  BUF  ARRAY.  OECIMAL 
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C NUMBERS  ARE  BROKEN  INTO  THE  TENS  AND  UNITS  OICITS.  WHICH  ARE 

C LOAOEO  AS  CHARACTERS  INTO  THE  BUF  ARRAY  AT  LOCATIONS  DETERMINED 

C BY  THE  COUNTER.  NN. 

C 

315  DO  327  1=2.26 
NN=5*( I-l I ♦! 

NUN* IF IXIZI I . J 1*0.5 > 

IFlNUM.LE.O)  00  TO  327 
IFINUM.LT. 10 1 GO  TO  325 
IFINUM.LT. 1001  GO  TO  323 
ND1*NUM/I 00* 1 
8UF{NN-2»=0I ClNDI ) 

323  ND2=M0D!NUM, 100)/I0frl 
BUF( NN- I | =Dl GI ND2 I 
325  NO 3= MODI NUM. 10171 
BUF!NN)=01G<  N03I 
327  CONTINUE 
C 

C AFTER  THE  ENTIRE  LINE  OF  DATA  IS  LOADED  I IT  COVERS  50  KM  EAST- 

C WEST  AT  A GIVEN  OISTANCE  NORTH!.  THE  LINE  IS  PRINTED  BY 

C PRINTING  THE  BUF  ARRAY. 

C 

WRITE (6. 401)  I BUF I 1 1 I I .1 11*1.126) 

330  CONTINUE 
C 

C THIS  SECTION  PRINTS  THE  CAPTIONS  FOR  THE  MAPS.  CAPTIONS  DEPENO 

c on  the  ileft  and  joown  entered. 

c 

IWEST*IL-2 
I W*I ABSI 1 VEST ) 

J SOUTH* JO- 2 
J S*IABS< JSOUTH) 

IFC IWEST.GE.OI  GO  TO  336 
IF( JSOUTH. GE.O)  GO  TO  340 
CO  TO  341 
C 

336  I F I JSOUTH. GE .0 ) GO  TO  33B 
GO  TO  339 

330  WRITEI6.402)  IV. JS 
60  TO  342 

339  WRI TE I 6*40  3)  IW.JS 
CO  TO  342 

340  WRI TEI 6.404)  IW.JS 
GO  TO  342 

341  WRITE! 6*405)  IW.JS 

342  GO  TO  <343.345. 347. 349*351 I.MM 

343  WRI TE! 6.406 ) ND.NT.LVL 
GO  TO  352 

345  WRITE! 6.407)  NO. NT 
GO  TO  352 

347  WRITE! 6.408)  NO. NT 
GO  TO  352 

349  WRITE! 6.409)  NO. NT 
CO  TO  352 

351  WRI TE! 6.41 0 ) NO. NT 

352  WRITE! 6.400) 

C 

C THE  RIGHT  HALF  OF  THE  MAP  IS  LOAOEO  AND  PRINTED  LIKE  THE 

C LEFT  HALF  WAS. 

C 

OO  377  JJ= 1.51 
J=S2-JJ 


uvuu 


00  3S7  1=1.126 
BUFI I )=BLNK 
397  CONTINUE 

c 

IFIMOOI! J-l) .Sl.NE.O)  CO  TO  362 

C 

OO  361  1*1 .1 26*25 
BUF(  I ) =PLUS 

361  CONTINUE 

C 

362  DO  375  1*27.51 
INN* 1*25 
NN*5*« INN-1 )*1 
NU<4<IFIX< 2 ( I . Jl+O.S) 

IFINUM.LE.O)  GO  TO  375 
IFINuM.LT. 10)  GO  TO  373 
IFINUM.LT. 100)  GO  TO  371 
NOKNUM/IOO+1 
BUF(NN-2)=0IG<NDI ) 

371  N02=M00(NUM, 1001/10.1 

BUFiNN-l ) = 01 G<  NO 2 ) 

373  NO3=M0O<  NUM. 10)4  1 
BUF<NN1<01G<  N03) 

375  CONTINUE 

C 

WRI TE<  6. 40 1 ) < BUF  < II  I ).I 11  = 1.126) 
377  CONTINUE 
C 

WRI TCI  6*4 1 1) 

WRI TE  <6.412)  NO.NT.LVL 


400 

401 

FORMAT! IHI I 

F0RMATI1H  • 1 26 A 1 • /// ) 

402 

FORMAT!*  * . • LE  * . / //* 
C*  KM  WEST  ANO  * . 1 3. • 

•••THE  LOWER  LEFT 
KM  SOUTH  OF  CLL • ) 

HA  NO 

GRID 

HARK 

IS 

• .13. 

403 

FORMAT!*  • ,*  LE • .///• 
C*  KM  WEST  AND  * . 13. • 

• , • T HE  LOWER  LEFT 
KM  NORTH  OF  CLL*) 

HANO 

GRIO 

MARK 

IS 

••13, 

404 

FORMAT! • • .*LE* .///• 

C*  KM  EAST  ANO  • • I 3. • 

* , • THE  LOWER  LEFT 
KM  SOUTH  OF  CLL*) 

HA  NO 

GRIO 

HARK 

IS 

• .13. 

405 

FORMAT!*  • ••LE* .///• 
C*  KM  EAST  ANO  * . 13. • 

• .'THE  LOWER  LEFT 
KM  NORTH  OF  CLL* ) 

HANO 

GRID 

MARK 

IS 

* . I 3. 

406  FORMAT I * • •• DATE* . A3. • T I ME • . II 6.4X . • REFLECT  1 V I TV  VALUES.  Z = 10** 
CIMAP  VALUE/10).  AT  THE*. 16.*  FT  LEVEL . • . 3X. • 2KM  X 2KM  GRID 
CINTERVAL.* ) 

407  FORMAT!*  • .'DATE* .A3. • T I ME * . 1 1 6. 4X. • PARTI AL  VI.  FROM  SURFACE  TO 
C15000  FT  < IN  DBZ  > .* ) 

408  FORMAT!*  • .* DATE* .A3. • T IME * . 1 I 6.4X. • PART I AL  VIZ  FROM  15000  TO 
C35000  FT  ! IN  DBZ).* I 

409  FORMAT!*  * . • DATE * . A3, • T I ME* . 1 1 6. 4X. ' PART I AL  VIZ  FROM  35000  TO 
C50000  FT  I IN  OBZ ) . • ) 

410  FORMAT!*  • .* DATE* .A3. • T IME* . 1 1 6.4X. 'TOTAL  VlL  FROM  SURFACE  TO 
C50000  FT  !1N  GM/M442 ) . • ) 

411  FORMAT! IH  .124X.2HRE) 

412  FORMAT!*  • .* OATE* .A3. • T IME * . I I 6.4X. 'LEVEL* . I 6. • FT. • » 

RETURN 

END 


SUBROUTINE  loaoicalibi 
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C THIS  SUBROUTINE  RECONSTRUCTS  THE  RADAR  CALIBRATION  CURVE  AND 

C CONVERTS  THE  256  DIGITAL  INTEGERS  TO  EQUIVALENT  DBM  VALUES, 

C IT  ALSO  PRINTS  OUT  THE  CAL  IB  ARRAY. 

C 

DIMENSION  OIGITSIBI, CALI B( 2 55  I 
READ! 5. I 1 00(01  GITS. DBM 
1100  FORMAT  !8F4.0«F6.2) 

IOIG  = IF  1 X(  DIGITS!  1 > ) 

c 

00  I I * I. IOIG 
I CAL  IB!  I I s DOM 

C 

DBM  » DBM  ♦ 10.0 

c 

DO  2 J * 2.8 
JJ  « J-l 

DBM  m DBM  - 1 0.0 
OELO  > OIGITS(J)  - DIGITSIJJI 
OELCAL  = 10. O/DELO 
IDELD  a IF  lx! OELDI 
C 

DO  3 K > 111 OELO 
KK  c DIGITS! JJ  > ♦ K 

3 CAL IB! KK ) = CAL  IB! KK* I I - OELCAL 


l CONTINUE 

• 

WR1TEI6. 20001  CAL  1 8 
2000  FORMATI'O* .20F6.1  » 
RETURN 
END 


FOLLOWING  THE  OECK  COMES  MORE  JCL  AND  THE  DATA  CAROS. 


X/GO.SYSIN  DO  • 

C • 

C DON'T  CHANGE  THIS  JCL  CARD. 

c 

C NEXT  COME  THE  DATA  CAROS.  FIRST  IS  THE  RADAR  CALIBRATION  DATA 

C card,  then  the  cards  with  the  initializing  values  for  the  scan 

C SEQUENCE.  ONE  SUCH  CARD  IS  REQUIRED  FOR  EACH  SCAN  SEQUENCE 

C THAT  IS  TO  BE  PROCESSED.  EXAMPLES  FOLLOW. 

C 

32  3T  60  98  139  171  201  233  107. S 

C 

C RAO AR  CALIBRATION  DATA 


0021 020691 81 002020 02201 0082 00230 
162-3826935100402002321008200238 

c 

C SCAN  SEQUENCE  CARDS  - 1 FOR  EACH  SCAN  TO  BE  PROCESSED 

c 

C the  FOLLOWING  JCL  CARDS  OEFINE  AND  CONTROL  THE  INPUT  AND 

C OUTPUT  TAPES. 

c 

X/FT08F001  DO  UN  1T-TAPE9. VOL-SER* ZZS6 IB.D I SP*! OLD. PASS! . 

X/  LABEL  — I t.NL*. IN  I , OCB-! RECFM=FB.LRECL=884 .8LKSI ZE-884 I 


136 


c 

C FT08  IS  THE  UNIT  NUMBER  { I RR I . FOO  I INDICATES  FOR  WHICH  RUN 

C THROUGH  THE  PROGRAM  THIS  INPUT  IS  TO  BE  USED.  ANO  ZZ5618  IS 

C THE  TAPE  NUMBER.  THE  SECOND  CARD  IS  A CONTINUATION  OF  THE 

C FIRST  AND  MUST  AL WATS  ACCOMPANY  IT.  L ADELsf 1 .... • INDICATES 

C WHICH  FILE  ON  THE  TAPE  YOU  WANT  TO  READ.  IN  THIS  CASE.  IT  IS 

C THE  FIRST  FILE.  FOR  THE  SECOND  SCAN  IBELOWI  IT  IS  THE  FIFTH 

C FILE.  NOTHING  ELSE  ON  THESE  CARDS  CHANGES. 

C 

X/FT08F002  DO  UN IT  = T APE9 . VOLsSER= ZZ56 1 8 .0  I SP=( OLD. PASS) . 

X/  LABELS ( 5. NL. . IN  I , DCB=( RECFM=FB ,LRECL=884 ,8LK SIZE  =884 > 

C 

C THE  FOLLOWING  PAIRS  OF  CARDS  ARE  FOR  THE  OUTPUT  TAPE. 

C 

X/FTIOFOOI  00  UN IT=T APE9. VOL=SER=  ZZ4 971 .LABEL=*<  1 5. NL. .OUT I . 

X/  DISP=( NEW. PASS) .DCB=<  RECFMsFB. LRECL— l 445. OLK SIZE* 1 3005 ) 

C 

C ONCE  AGAIN.  FTIO  IS  THE  UNIT  NUMBER  (IRIT).  AND  FOOI  IS  THE 

C SAME  AS  ABOVE.  OUTPUT  TAPE  NUMBER  IS  ZZ497I.  L A8EL=< l 5 . . . . . 

C MEANS  THAT  THE  OUTPUT  OF  THE  FIRST  PROGRAM  RUN  (5  MAPS ) WILL 

C BE  FILE  415  ON  THIS  TAPE.  LIKEWISE.  THE  OUTPUT  FROM  THE 

C SECOND  RUN  WILL  BE  FILE  *16.  NOTHING  ELSE  ON  THESE  CAROS 

C CHANGES. 

C 

X/FTIOF002  DO  UNIT=TAPE9. VOL=SER= ZZ497I .LABEL*! 16. NL. .OUTI . 

X/  DISP=<NEW.PASS> . DCB=<  RECFMsFB. LRECL-3  1445. GLKS I Z£* 13005) 

C 

X*ENO 

C THIS  IS  ALWAYS  THE  LAST  CARO  IN  THE  DECK. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


* RORMAP • - A CONREC  PLOTTING  AND  CONTOURING  PROGRAM 

THIS  PROGRAM  TAKES  THE  100  X 100  KM  OUTPUT  BOXES  <51  X 51 
ARRAYS)  OF  REFLECTIVITY  AND  VIL  FROM  THE  • RADAR*  PROGRAM  <APP. 
A)  AND  PLOTS  AN  AREA  OF  INTEREST  SPECIFIED  BY  THE  PROGRAMMER 
ON  THE  CALCOMP  PLOTTER.  THE  OUTPUT  IS  IN  THE  FORM  OF  A CON- 
TOURED MAP  OF  REFLECTIVITY  OR  VIL.  THE  VARIOUS  PARAMETERS  CAN 
ALL  BE  SPECIFIED  BY  THE  PROGRAMMER.  AND  THE  PROGRAM  IS  SUFFIC- 
IENTLY GENERAL  THAT  IT  CAN  BE  USED  FOR  ANY  GENERL  A I ZED  PLOT 
ROUTINE.  NOT  JUST  THE  •RADAR*  PROGRAM  USED  HERE. 

SOURCES  FOR  THIS  PROGRAM  ARE: 

(II  RE I O - RM  6 1 7 A OEM  BLDG 

<21  NCAR  LIBRARY  ROUTINE  TECHNICAL  MANUAL 

( 3 I OS/360  LOCAL  APPENDIX 

( 4 ) CONREC  APPENDIX 


THE  MECHANICS  OF  RUNNING  THE  PROGRAM  ARE  THE  SAME  AS  FOR  THE 
•RADAR*  PROGRAM  WITH  THE  FOLLOWING  ADDITIONAL  REQUIREMENT. 

normally  you  will  only  provide  the  input  tape  ano  the  dpc  will 

USE  ONE  OF  THEIR  TAPES  FOR  THE  OUTPUT  TAPE.  WHICH  IS  CALLFO 
•CALCOMP*.  ACCORDINGLY.  WHEN  YOU  FILL  OUT  THE  JOB  TICKET, 

ENTER  'CALCOMP*  AS  THE  OUTPUT  Tape.  YOU  MUST  ALSO  FILL  OUT 
A 'CALCOMP  JOB  TICKET*.  WHICH  IS  AVAILABLE  AT  THE  WINDOW  OF 
THE  MACHINE  ROOM  WHERE  YOU  TURN  THE  TAPES  IN.  THE  SPACE  MARkEO 
•TAPE  VOLUME*  ON  THIS  TICKET  IS  FOR  THC  OUTPUT  TAPE  AND  SHOULD 
BE  MARKED  'CALCOMP*.  UNLESS  YOU  PROVIDE  THEM  WITH  YOUR  Own 
OUTPUT  TAPE.  THE  FINISHED  PLOTS  MAY  BE  PICKED  UP  IN  THE  I/O 
ROOM.  OR  THEY  WILL  BE  DELIVERED  TO  TOUR  BOX  ADDRESS  WHEN  THEY 
ARE  COMPLETED. 


*4**44*4*44*4  JCL  - JOB  CONTROL  LANGUAGE  *4*4***4444*** 

WITH  THIS  JCL.  THE  PROGRAM  RUNS  IN  'FORTRAN  IV  G LEVEL  21*. 

THE  JCL  WILL  REQUIRE  SOME  MINOR  CHANGES  FOR  EACH  PROGRAM  RUN. 
NOTE  THAT  THE  FIRST  CHARACTER  ON  ALL  JCL  CARDS  <COL.  1)  SHOULD 
BE  A •/*.  REPRODUCED  HERE  AS  AN  *X*. 


X/RORMA P JOB  < X999.0060.002. 030 .XX) ,*YOUR  NAME 


c 

C 002  IS  TIME  IN  MINUTES.  DON'T  CHANGE. 

C 

x*p assword  xxxx 

X*SETUP 


c 

C THIS  TELLS  THE  MACHINE  NOT  TO  RUN  UNTIL  THE  TAPES 

C HAVE  BEEN  MOUNTED. 

C 

X4LEVEL  X 

X*J08PARM  K=0.R=256.T=2.TAPE9=2 
X/  EXEC  FTNGPLOT  .REG I 0n*=256K 
X/F ORT .SYS IN  OD  * 

c 

C DON* T CHANGE  THESE  LAST  THREE  CARDS. 

C 

c 

c 

C •*•«•*»*•*••*  THE  MAIN  PROGRAM  OECK  BEGINS  HERE. 

C 


O O U U uuou  oouuu  uuuuu  uuo  uou  o uuuuu  uu 


ONES  USED  TO  PROOUCE  THE 
HAVE  BEEN  REDUCED  TO 


C ******  NOTE  ♦*••*** 

C THE  PARAMETERS  SPECIFIED  HERE  WERE  THE 

C • NEAR  RANGE*  MAPS  IN  THIS  PAPER.  WHICH 

C 65%  OF  THE  SIZE  OF  THE  ORIGINALS* 

C A DETAILED  EXPLANATION  OF  THE  PARTS  OF  THIS  PROGRAM  CAN  BE 

C FOUND  IN  CHAPTER  tit* 

c 

c 

c 

DIMENSION  ZARRAYI51 .51 ). ZI26.26) 

C 

C ZARRAY  IS  THE  OUTPUT  ARRAY  FROM  THE  ‘RADAR*  PROGRAM,  AND  IS 

C IOO  X 100  KM  ON  A 2 KM  GRID.  Z IS  THE  PORTION  OF  ZARRAY  THAT 

C IS  TO  BE  PLOTTED  ANO  CONTOURED  OY  THIS  PROGRAM.  ITS  SIZE  ANO 

C POSITION  WITHIN  ZARRAY  ARE  DETERM I NOD  BY  THE  PROGRAMMER. 

C ZI26.26I  SPECIFIES  AN  AREA  SO  X 50  KM. 

C 

COMMON/CONREI/SIZEt-.SIZEM.SlZEP.NLA.NLM.XLT, YB T , SIDE. NREP.NCRT. ILA 
♦ 8.NULBLL. IOFFO.EXT. IOFFP.SPVAL. IOFFM. I SOL  10 

COMMON  STATEMENT  FOR  CONREC.  00  NOT  CHANGE.  FOR  DETAILS  ON 
THE  MEANING  OF  EACH  ARGUMENT.  SEE  CONREC  APPENDIX. 

SIZEL=»1.97 

SIZEL  - SIZE  OF  LINE  LABELS.  1.31=2MM.  OEFAULT  OPTION  (NO 
CARO ) IS  1.51 

NULBLL=I 

NULBUL  IS  THE  NUMBER  OF  UNLABELEO  LINES  BETWEEN  LABELED  LINES. 
DEFAULT=3,  OR  EVERY  FOURTH  LINE  IS  LABELED. 


INOEX-O 

Index  is  the  count  of  how  many  files  (scan  sequences i of  data 

ARE  READ  FROM  THE  RAOAR  TAPE.  REMEMBER.  THERE  ARE  FIVE  MAPS 
PER  SCAN  SEQUENCE. 

DO  99  ((>1.3 

THIS  RUN  PROCESSED  THREE  FILES  (SCAN  SEQUENCES!  OF  DATA. 
IC0UNT=0 

ICOUNT  IS  THE  MAP  NUMBER.  THERE  ARE  FIVE  MAPS  PER  FILE. 

I NO£X= I NDE  X+l 

ICOR=0 

JC0R»25 

ICOR  ANO  JCOR  ARE  USED  TO  SHIFT  THE  STARTING  POINT  OF  Z (LOWER 
LEFT  CORNER)  TO  A PARTICULAR  POINT  IS  ZARRAY.  FOR  EXAMPLE.  IF 
ICORsO  AND  JC0R>25.  THEN  Z( I • I )= ZARRAY! 1.26).  SEE  NOTE  ABOVE* 

15  CONTINUE 

READ! 8 • 6006. END*99 I ZARRAY 
06  FORMAT (200FS.I. 89F5. I I 

THIS  IS  THE  SAME  FORMAT  USED  ON  TAPE  IN  THE  RAOAR  PROGRAM. 


T 
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I COUNT  * I COUNT ♦ 1 

00  I 1*1 *26 
00  1 J*1.26 

ic*i*icor 

JC*J*JC0R 

ZU • J)*ZARRAYt IC.JCI 
CONTINUE 

THIS  TRANFERS  THE  OATA  TO  BE  CONTOURED  INTO  THE  Z ARRAY. 

CALL  SET! I 10.980. 1 1 0. 980 • 0. .25. .0. • 25. . I ) 

THIS  SUBROUTINE  PHYSICALLY  SETS  UP  THE  GRAPH.  SEE  CHAP.  III. 
CALL  PER IMi25«t,2S«l I 

this  subroutine  draws  the  oox  around  the  area  and  places  tick 

MARKS  AR0UN0  THE  PERIMETER  OF  THE  BOX.  SEE  CHAP.  III. 

CALL  PWRXI40.54S.28H*PRU*K1L0METERS  SOUTH  OF  CLL. 2B .I..1.1) 

CALL  PWRXI545.50.27H*PRU*KILOMETERS  EAST  OF  CLL . 27. 1 * . 0. 1 ) 

THIS  SUBROUTINE  PRINTS  THE  LABELS  ON  THE  MAP  AND  MUST  BE  CALLED 
ONCE  FOR  EACH  LABEL.  SEE  CHAP.  III. 

CALL  PWRXI 25. .90.7H*  I RU*  50 . 7.  I . , 0 . I ) 

CALL  PWRXI20..90.7H* I RU *40. 7. I. ,0.1) 

CALL  PWRXI I5..90.7H* IRU* 30.7.1. .0.1) 

CALL  PWRXI  10. .90. 7H* IRU*  20. 7. 1 . .0. I ) 

CALL  PWRXI  5. .90. 7H* IRU* 10. 7. I . .0. 1 i 
CALL  PWRXI 0. .90.6H* IRU*0.7. 1.  ,0.1  I 

THIS  PUTS  THE  NUMBERS  ALONG  THE  X AXIS. 

CALL  PWRXI  100.  0. . 7H • I RU *50.7, I..0.31 
CALL  PWRXI 100.  5..7H* IRU*40.7.1..0.3> 

CALL  PWRXI 10  0. 10. ,7H*  IRU*  30.7, 1 . ,0.3) 

CALL  PWRXI 100. 15. • 7 H • IRU1 20.7. 1..0.3I 
CALL  PWRXI 100.20. ,7H* IRU* 10.7.1.. 0.3) 

CALL  P WRXI100. 25. ,6H*IRU *0.7.1. .0.3) 

THIS  PUTS  THE  NUMBERS  ALONG  THE  Y AXIS. 


IF I I COUNT • GT .4  I GO  TO  5 

THIS  SELECTS  A DIFFERENT  SET  OF  CONTOUR  INTERVALS  FOR  THE  VI L 
MAP. 


C 

C 


CALL  CONRECt  Z<  2 . 2 1 . 26 . 24 . 24 . 1 0. .60. .5..1.-1.0) 

# 

THIS  SUBROUTINE  ACTUALLY  DOES  THE  PLOTTING  AND  CONTOURING  OF 
THE  DATA.  AND  IN  THIS  CASE,  CONTOURS  THE  REFLECTIVITY  MAPS 
FROM  10  TO  60  IN  5 UNIT  INCREMENTS.  SEE  CHAP.  III. 

GO  TO  10 

CONTINUE 

CALL  CONRECI ZI 2. 2 1 . 26 . 24 . 24. 1 . . 20. * I. . I .-I .01 


— - 


141 


c 

c 

c 

c 


THIS  TIME  THE  VIL  MAP  IS  DRAWN.  CONTOURED  FROM  1 TO  20  IN  I 
UNIT  INCREMENTS. 


10  CONTINUE 


I i 


C 

c 

c 

c 

c 


c 

c 

c 

c 


CALL  PLOT! 12.. 0. . -3) 

THIS  MOVES  THE  PEN  TO  A NEW  LOCATION  ANO  REINITIALIZES  THE 
COORDINATE  SYSTEM  IN  PREPARATION  FOR  QEGINNING  A NEW  MAP.  SEE 
CHAP.  III. 

CO  TO  IS 

99  CONTINUE 

CALL  PLOT<0..0..9«> 

THIS  EMPTIES  THE  OUTPUT  BUFFER  AnO  IS  NECESSARY  TO  END  THE 

PROGRAM. 

STOP 

■NO 


••••••••••••  THE  REST  OF  THE  JCL  FOLLOWS  THE  DECK. 


C 
C 

c 
c 
c 
c 
c 
c 

x/lkeo.ooname  oo  oisp*shr, 

*/  OSN4ME-USER.OCN.GUINASSO. J0HLI8 
A/LKEO.SYStN  OO  • 

INCLUOE  OONAME I CNRCSMTH. PWRX I 

c 

c DON'T  CHANCE  THESE  LAST  FOUR  CAROS. 

c 

C NEXT  COME  THE  INPUT  CARDS  - ONE  SET  OF  THESE  CAROS  FOR  EACH 

C FILE  THAT  IS  TO  BE  PROCESSED. 

c 

X/cb.FToaFOOl  00  UNI T-TAPE9. vOl=SER=ZZ497 I . 

X/  LABEL*! 03.NL. .INI .DISP=!  OLO .PASS) . 

X/  OCB*!ReCFM*F8.LRECL*)445.0LKS!ZE=13005> 

c 

C FTOB  IS  THE  UNIT  NUMBER  THAT  REFERS  TO  THE  INPUT  TAPE. 

C F00I  IS  The  FIRST  FILE  TO  BE  PROCESSEO. 

C LABEL* I INDICATES  THAT  THIS  FILE  IS  THE  THIRD  FILE  ON 

C the  input  tape. 

C ZZ497I  IS  THE  10  NUMBER  OF  THE  INPUT  TAPE. 

C 

X/GO.FT 08F002  DO  UNI T*T AP£9. V0L=SER*ZZ497 1 . 

X/  LABEL-! 05 .NL. . I N> .0  I SP= i OLO .PASS  I • 

X/  OCB*i RECFM-FB .LRECL—1 445. BLkSI ZE* 13005) 

c 

C the  SECOND  FILE  ISCAN  SEQUENCE)  TO  BE  PLOTTED  IS  THE  FIFTH 

C FILE  ON  THE  TAPE. 

C 

X/GO.FT08F003  00  UNI T*T APE9 . V0L-SER-ZZ497 I . 

X/  LABEL*! 06.NL..IN) .0  I SP  = ! OLD .PASS ) . 

X/  OCO* 1 RECFM-FO .LRECL* 1445. BLKS I ZE* 13005) 

c 
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